Question: Gapit error in R
0
gravatar for mmh7272
15 months ago by
mmh727240
mmh727240 wrote:

Hi fellow happy new year i have problem. i use Gapit tool pupose of GWAS. i have phenotype data.txt and snp.hp file so i use gapit tool this is script when i use in R

===================================================

library('MASS')
library(multtest)
library(genetics)
library(gplots)
library(compiler)
library("scatterplot3d")
myY <- read.table("/home/mmh7272/cp/t18/sand1.txt", head=T)
source("/home/mmh7272/gapit/gapit.txt")
source("/home/mmh7272/gapit/gapit_functions_MLM.txt")
myG <- read.table("/home/mmh7272/cp/292.snp.out.hp", head=F)
myG <-myG[,c(1:11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181,182,183,184,185,186,187,188,189,190,191,192,193,194,195,196,197,198,199,200,201,202,203,204,205,206,207,208,209,210,211,212,213,214,215,216,217,218,219,220,221,222,223,224,225,226,227,228,229,230,231,232,233,234,235,236,237,238,239,240,241,242,243,244,245,246,247,248,249,250,251,252,253,254,255,256,257,258,259,260,261,262,263,264,265,266,267,268,269,270,271,272,273,274,275,276,277,278,279,280,281,282,283,284,285,286,287,288,289,290,291,292,293,294,295,296,297,298,299,300,301,302,303)]
source("/home/mmh7272/gapit/gapit.txt")
source("/home/mmh7272/gapit/gapit_functions_MLM.txt")
GAPIT <- GAPIT(Y=myY, G=myG, PCA.total=3, SNP.MAF=0.05, Model.selection=TRUE)

==================================================================

it happens error

[1] "Genotyping: numericalization, sampling kinship, PCs and much more..."
[1] "Converting genotype..."
[1] "Converting HapMap format to numerical under model of Middle"
[1] "Perform numericalization"
  **Error in if(bit == 1) { : missing value where TRUE/FALSE needed**
In addition: Warning messages:
1: 'memory.size()' is Windows-specific
2: 'memory.size()' is Windows-specific
3: 'memory.size()' is Windows-specific
4: 'memory.size()' is Windows-specific
snp R • 562 views
ADD COMMENTlink modified 15 months ago by Kevin Blighe41k • written 15 months ago by mmh727240
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1373 users visited in the last hour