我是R的新手,我想根据不同的计数表,使用maSigPro
包进行一些时间序列分析。
目前我正在努力处理p.vector
函数的数据类型。
具有TPM标准化计数的表格如下:
counts_TPM = as.matrix(read.table("TPM_Counts_25min.txt", header=T,
row.names=1, stringsAsFactors=FALSE)
counts_TPM=as.data.frame(counts_TPM, stringsAsFactors=FALSE)
typeof(counts_TPM)
# [1] "list"
此数据类型不适合p.vector函数:
fit=p.vector(counts_TPM, sample_des_design, Q=0.05, counts=TRUE)
# **Error** in FUN(newX[, i], ...) : invalid 'type' (character) of argument
这是我的TPM标准化计数表的简短版本:
E15R003a06 E15R003b06 E15R003c06 E15R003d13 E15R003d12
A 12.554467 12.68203563 1.436166e+01 13.58263382 1.312083e+01
B 10.392282 14.65197068 1.141320e+01 8.84738000 8.826718e+00
R 39.127930 42.10897757 4.194345e+01 34.76635026 3.005847e+01
X 6.696018 6.64936858 5.897263e+00 4.98501460 4.652471e+00
S 77.418782 84.95649902 8.017746e+01 82.99038040 7.884357e+01
T 53.565504 56.02217626 5.487839e+01 49.47038414 4.341778e+01
我还尝试read.csv
使用我的表格的csv格式和data.matrix
之类的东西(但这个功能完全弄乱了我的价值观),但我想我错过了一个要点......
我希望有人可以帮助我!