.subset(x，j)中出错:下标类型'list‘无效

Question

我是R的新手，如果有人能帮我，那就太好了。

我使用的是一个叫NOISeq的程序，这是一个基于R的软件，可以检测许多因素之间的差异基因表达。我读入了一个数据框的表，它看起来像这样：

            X10G48 X35Y87 X36W26 X23Y79 X2B84 X12Y30 X10B70 X10G87 X36W62
XLOC_000001     33     34     39     74    43     43     34     28     42
XLOC_000002     17     42     44     67    38     58     41     29     25
XLOC_000003      0      0      0      0     0      0      6      0      0
XLOC_000004      0      0      0      2     0      0     10      0      0
XLOC_000005     44     57     37     71    45     47     49     53     36
XLOC_000006     43     46     42     71    53     53     49     48     18
            X23Y70 X2UNA X12Y47 X10R99
XLOC_000001     82    38     28     23
XLOC_000002     58    53     28     27
XLOC_000003      0     0      0     12
XLOC_000004      0     0      4      2
XLOC_000005     47    67     48     39
XLOC_000006     53    61     37     26

该表按两个因子排序，使得10G48, 35Y87, 36W26, 23Y79, 2B84, 12Y30, 10B70都是条件1，10G87, 36W62, 23Y70, 2UNA, 12Y47, 10R99都是条件2。

我使用了下面的代码：

library(NOISeq)
setwd("/home/user/edgeR")
data.frame <- read.table("nalphavbeta.txt", header=TRUE, sep='\t', row.names=1)
myfactors=data.frame(c(1,1,1,1,1,1,1,2,2,2,2,2,2))
sam<-readData(data=data.frame, factors=myfactors)
myRPKM = rpkm(assayData(sam)$exprs, k = 0, lc = 1)
head(myRPKM[, 1:4])
mynoiseqbio = noiseqbio(sam, k=0.5, norm="rpkm", factor=NULL, lc = 1, r = 20, adj = 1.5, plot = FALSE, a0per = 0.9, random.seed = 12345, filter = 2) 

但是它返回错误

Error in .subset(x, j) : invalid subscript type 'list'

我有一种感觉，这与因子的论点有关，但我不确定到底是什么。任何帮助都将不胜感激-谢谢！

Answer 1

问题已解决

我只需要改变一下台词

myfactors=data.frame(c(1,1,1,1,1,1,1,2,2,2,2,2,2))

至

myfactors=data.frame(caste= c(1,1,1,1,1,1,1,2,2,2,2,2,2))

例如，为因子命名，然后将因子级别添加到函数中

mynoiseqbio = noiseqbio(sam, k=0.5, norm="rpkm", factor=caste, lc = 1, r = 20, adj = 1.5, plot = FALSE, a0per = 0.9, random.seed = 12345, filter = 2)