[R] bootstrap lmekin model

2012-11-24 Thread Roberta Carriero
Hi,I use the 'lmekin' model of the 'kinship' package of R in order to estimate heritability. I want to estimate the confidence interval of the variance coefficient and so I should use a bootstrap simulation. The pedigree file has 1386 subjects so I create a kinship matrix [1386*1386].This is th

[R] Bootstrap lmekin model

2012-11-24 Thread Roberta Carriero
Hi,I use the 'lmekin' model of the 'kinship' package of R in order to estimate heritability. I want to estimate the confidence interval of the variance coefficient and so I should use a bootstrap simulation. The pedigree file has 1386 subjects so I create a kinship matrix [1386*1386].This is th

[R] FW: replace repeated id in a pedigree list

2012-10-30 Thread Roberta Carriero
#[253] 543 943 947 192 949 945 944 410 333 427 > > #The same way you can get the repeated FAID and MOID > dat1$FAID[duplicated(dat1$FAID)|duplicated(dat1$FAID,fromLast=TRUE)] > dat1$MOID[duplicated(dat1$MOID)|duplicated(dat1$MOID,fromLast=TRUE)] > > From the dat3 IDs, you have 6

[R] replace repeated id in a pedigree list

2012-10-29 Thread Roberta Carriero
Hello, I have a pedigree file such this: FAMIDIDFA IDMO IDSEXSTATUS 110020 123122 143112 153110 163110 173122 183110 19311