david.beede at mail.doc.gov wrote:>
> Dear R list:
>
> I want to make separate estimates for each level of the variable
"group."
> After consulting many sources I am stumped as to why the following does not
work:
>
> > wintemp <- subset(alltemp, winner==1)
> > my.ppr <- function(x)
> + {
> + if(nrow(x) >= 50) {
> + pprfile <- ppr(award~ilogemp,
data=x,nterms=5,max.terms=10,optlevel=3)
> + summary(pprfile)
> + }
> + }
> > test.by <- by(wintemp,as.factor(wintemp$group),my.ppr)
> Error in model.frame.default(formula = award ~ ilogemp, data = x) :
> Object "x" not found
>
> Any help would be greatly appreciated.
> (I am sure I am misunderstanding some fundamental R concepts.)
It is hard to help, if we don't know what all those variables stand for.
I tried something similar and it worked.
What version of R are you running?
Uwe Ligges
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at
stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._