rrBLUP / GWAS错误:"需要数字/复杂矩阵/向量参数"

时间:2014-09-18 03:34:21

标签: r

所以,我在这里看到了一些与我收到的错误问题有关的问题,但我尝试的解决方案似乎没有用。

我收到的消息是“需要数字/复杂矩阵/矢量参数”

以下是代码示例:

> geno <- read.table("GENO2.txt", header=TRUE)
> head(geno)

这里的最终目标是在植物系统中进行关联作图(遗传学)。我正在使用包含rrBLUP的函数GWAS。除其他外,GWAS需要表型和基因型数据(例如pheno和geno)

> GWAS.result <- GWAS(pheno, geno, fixed=NULL, K=A, n.PC=1, min.MAF=0.05, n.core=1, P3D=TRUE, plot=TRUE)

这个功能GWAS是我收到错误的地方。

通过尝试解决这个问题的过程,我发现R添加了一列值(1:526),这是我拥有的遗传标记的数量。当我添加我的标题信息时,这具有最后一列(对应于工厂行#301)的效果是空白的。上面的代码是重新加载这个文件,我在Excel中编辑了该文件以使数据匹配。

问题是,它仍然无效。我试图将“pheno”数据框中的变量从因子更改为数字,但这也不起作用。我应该补充一点,“pheno”有2个coumns,一个是植物系ID信息,另一个是疾病评级。具有工厂行ID的列对应于geno文件中的列标题。

非常感谢任何帮助或建议。

1 个答案:

答案 0 :(得分:0)

马修: 我认为你如何设置你的pheno和geno数据集的问题是:

geno数据集应该在行中有标记,你应该包括三个col,它们是&#34; marker&#34;,&#34; chrom&#34;,&#34; pos&#34;加上基因型的名称。

至于pheno数据集: 第一个col应该包含基因型ID加上你的变量

这是Jeff Endelman(http://potatobreeding.cals.wisc.edu/wp-content/uploads/sites/21/2014/01/GWAS_tutorial.pdf

提供的示例

希望有所帮助