所以,我在这里看到了一些与我收到的错误问题有关的问题,但我尝试的解决方案似乎没有用。
我收到的消息是“需要数字/复杂矩阵/矢量参数”
以下是代码示例:
> 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文件中的列标题。
非常感谢任何帮助或建议。
答案 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)
提供的示例希望有所帮助