Re: [R] an issue about missing data

2010-08-05 Thread David Winsemius
On Aug 5, 2010, at 3:59 PM, karena wrote: library(GenABEL) gs.b <- gs > library(GenABEL) Loading required package: MASS GenABEL v. 1.6-0 (June 21, 2010) loaded > gs.b <- gs Error: object 'gs' not found You are not making this easy. I asked what str(gs) returned. Instead you gave us code b

Re: [R] an issue about missing data

2010-08-05 Thread karena
library(GenABEL) gs.b <- gs ok <- complete.cases(g...@phdata[,c('sex','age','b.dbp','b.bmi')]) gs.b <- gs.b[ok] g...@phdata <- g...@phdata[,c('id','sex','age','b.dbp','b.bmi')] index=1:g...@gtdata@nsnps ran.snp=sample(index,261,replace=F) gs.b.gkin=ibs(gs.b[,ran.snp], weight="freq") attach(g...@p

Re: [R] an issue about missing data

2010-08-05 Thread David Winsemius
On Aug 5, 2010, at 3:44 PM, karena wrote: thank you, but even after I tried 'complete.cases' function, How did you try? I still get the same error messages. What does str(gs) tell you? help. Plaintive cries for help are considerably less useful than details about the

Re: [R] an issue about missing data

2010-08-05 Thread karena
thank you, but even after I tried 'complete.cases' function, I still get the same error messages. help. -- View this message in context: http://r.789695.n4.nabble.com/an-issue-about-missing-data-tp2315094p2315519.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] an issue about missing data

2010-08-05 Thread David Winsemius
On Aug 5, 2010, at 12:23 PM, karena wrote: Hi, I am using a function 'polygenic' in the package called 'GenABEL'. h2.gs <- polygenic(b.dbp~age+age2+age3+sex+b.bmi, kin=gs.gkin, data=gs) Have you considered us

[R] an issue about missing data

2010-08-05 Thread karena
Hi, I am using a function 'polygenic' in the package called 'GenABEL'. > h2.gs <- polygenic(b.dbp~age+age2+age3+sex+b.bmi, kin=gs.gkin, data=gs) --