Displaying 3 results from an estimated 3 matches for "ilogage".
2001 Jun 06
1
ppr, number of terms, and data ordering
...y, in the same example in MASS, where does the test.cpus() function come from? I
couldn't find it in the MASS table of contents on CRAN.
Thanks in advance for any help you can give!
######DATASET1, ORIGINAL ORDER
> pprfile.ppr <- ppr(
+ award~
+ ilogemp+ilogage,
+ data=dataset1, nterms=1, max.terms=40, optlevel=3, bass=0
+ )
> pprfile.ppr
Call:
ppr.formula(formula = award ~ ilogemp + ilogage, data = dataset1,
nterms = 1, max.terms = 40, optlevel = 3, bass = 0)
Goodness of fit:
1 terms 2 terms 3 terms 4 terms...
2001 Mar 20
3
Newbie question about by() -- update
Sorry about the lack of detail. I am running R v.1.2.2. I can recast my
question
(which I think I have partially answered) more succinctly as follows:
1. This seems to work (note that group takes values 1,2,3,4, or 5):
my.newfun <- function(x) myfile <- lm(award ~ ilogemp + ilogage, x)
test.by <- by(wintemp, as.factor(wintemp$group), my.newfun)
2. This does not work (leaving aside whether I am using ppr correctly or not!):
> my.pprfun <- function(x) mypprfile <- ppr(award ~ ilogemp + ilogage, data = x, nterms = 5,
+ max.terms = 10, optlevel = 3)...
2001 Mar 28
4
efficiency and "forcing" questions
...; group <= 50 & winner==1))
rm(alldata4)
for(i in 1:4) gc()
library(modreg)
attach(wintemp4)
myppr1 <- function(x)
{
#run pprfile once to get list of sum of squared errors corresponding to differen numbers of terms
pprfile.ppr <- ppr(
award~
ilogemp+ilogage+sdb+allsmall+
size2+size3+size4+size5+size6+size7+size8+size9+size10+
X.Iprimnaic.2+X.Iprimnaic.3+X.Iprimnaic.4+X.Iprimnaic.5+X.Iprimnaic.6+
X.Iprimnaic.7+X.Iprimnaic.8+X.Iprimnaic.9+X.Iprimnaic.10+X.Iprimnaic.11+
X.Iprimnaic.12+X.Iprimnai...