Many Thanks Mark for your answer. It seems than, that it is not possible to use all variables without somehow imputing missing values.
But I will try which variables I can finally use. Many thanks again. B. Mark Difford wrote: > > Hi Birgitle, > > It seems to be failing on those columns that have just a single "entry" > (i.e = 1, with the rest as 0; having just 1, an <NA>, and then 0s gets you > through). And there are other reasons for failure (in the call to get a > positive definite matrix). > > The main problem lies in the calculation of standard errors for some > combinations. You can get a result for all columns by turning this off. > You can get standard errors for all the columns up to 60 by omitting > columns 12, 23, and 44. You need to work out the rest yourself by > "columning" forwards till you get a failure, then drop that column from > the index. There probably isn't enough information to calculate SEs for > the columns that cause an error. > > ## This works with default settings but without columns 12, 23 and 44 > hetcor(TestPart[,c(1:11,13:22,24:43,45)], pd=T, std.err=T, > use="complete.obs") > > ## The first fails; the second works > hetcor(TestPart[,c(1:11,13:22,24:43,45:60)], pd=T, std.err=F) > hetcor(TestPart[,c(1:72)], pd=F, std.err=F) > > HTH, Mark. > > ----- The art of living is more like wrestling than dancing. (Marcus Aurelius) -- View this message in context: http://www.nabble.com/Problems-using-hetcor-%28polycor%29-tp18867343p18870543.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
