Martin Elff <elff <at> sowi.uni-mannheim.de> writes: > > modco <- list( > lm(normskvop ~ I(nts^0.5)-1, data = colo,weights=wtz), > lm(normskvop ~ I(nts^0.5)-1, data = colo,weights=wtz, > subset=sector!="X"), > lm(normskvop ~ I(nts^0.5)-1, data = colo,weights=wtz, > subset=sector!="A"), > lm(normskvop ~ I(nts^0.5)-1, data = colo,weights=wtz, > subset=sector=="A"), > lm(normskvop ~ I(nts^0.5)-1, data = colo,weights=wtz, > subset=sector=="M") > ) > > sapply(modco,function(x) coef(summary(x)))
This works, but consider if lm(normskvop ~ I(nts^0.5)-1+sector, data = colo,weights=wtz), is a better approach anyway. Dieter ______________________________________________ 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.