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
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
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
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.
__
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
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)
--
6 matches
Mail list logo