Displaying 3 results from an estimated 3 matches for "makekinship".
2011 Apr 15
1
no solution yet, please help: extract p-value from mixed model in kinship package
...rep(8, 10),
rep(10, 10), rep(12, 10), rep(14, 10), rep(16, 10), rep(18, 10),
rep(20, 10), rep(22, 10) )
ped <- data.frame(id, dadid, momid)
# *****************************kmatrix**************************************
> cfam <- makefamid(ped$id,ped$momid, ped$dadid)
>
> kmat <- makekinship(cfam, ped$id, ped$momid, ped$dadid)
>
> #*****************************************x and y variables
*********************
> set.seed(3456)
>
> dat <- sample(c(-1,0,1), 10000, replace = TRUE)
>
> snpmat<- data.frame(matrix(dat, ncol = 100))
>
> names(snpmat) <- c...
2008 Dec 24
0
command Polygenic gives error message concerning dimensions of data
...r and mother; otherwise only the
brothers are given). Genotypes are provided for the largest part of the
members, whereas the phenotype is only given for the male individuals. I
imported the data using the convert.snp.ped and load.gwaa.data commands. The
kinship table was constructed with "makekinship". When I look at the
dimensions of the kinship table, it is exactly the same as the number of
individuals in my dataset.
I also provide you with the following information because maybe that gives a
clue on where things are going wrong. When using the genetic control method
via the qtscore...
2011 Apr 14
0
extract p-value from mixed model in kinship package
...7, 10),
rep(19, 10), rep(21, 10))
momid <- c(rep(0, 5), rep(2, 5), rep(4, 5), rep(6, 5), rep(8, 10),
rep(10, 10), rep(12, 10), rep(14, 10), rep(16, 10), rep(18, 10),
rep(20, 10), rep(22, 10) )
ped <- data.frame(id, dadid, momid)
cfam <- makefamid(ped$id,ped$momid, ped$dadid)
kmat <- makekinship(cfam, ped$id, 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.fr...