search for: covtrait

Displaying 2 results from an estimated 2 matches for "covtrait".

Did you mean: contraint
2011 Apr 15
1
no solution yet, please help: extract p-value from mixed model in kinship package
...nd y variables ********************* > set.seed(3456) > > dat <- sample(c(-1,0,1), 10000, replace = TRUE) > > snpmat<- data.frame(matrix(dat, ncol = 100)) > > names(snpmat) <- c(paste ("VR",1:100, sep='' )) > > yvar <- rnorm(100, 30, 5) > covtrait <- rnorm(100, 10, 5) > > mydata <- data.frame(id, yvar, covtrait, snpmat) > #******************************mixed model in lmekin ******************************************* > > fmod <- lmekin(yvar ~ mydata[,3] , data= mydata, random = ~1|id, > varlist=list(kmat)) $coeff...
2011 Apr 14
0
extract p-value from mixed model in kinship package
...d, ped$momid, ped$dadid) set.seed(3456) SNPdata <- c(1:1000) dat <- sample(c(-1,0,1), 10000, replace = TRUE) snpmat<- data.frame(matrix(dat, ncol = 100)) names(snpmat) <- c(paste ("VR",1:100, sep='' )) yvar <- rnorm(100, 30, 5) mydata <- data.frame(trait1, covtrait, snpmat) fmod <- lmekin(y~ mydata[,i] , data= snpdata, random = ~1|id, varlist=list(kmat)) $coefficients[2,4] # does not work Ultimately I want to put into the loop: for(i in 3:length(snpdata)) { P <- vector (mode="numeric", length = 1000) P[i] <- lmekin(trait1~ snpdata[,i]...