Displaying 5 results from an estimated 5 matches for "ilogemp".
2001 Jun 06
1
ppr, number of terms, and data ordering
...Finally, 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...
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, optl...
2001 Feb 12
2
supsmu vs. ppr
I used the supersmoother function in the modreg package as follows:
super <- supsmu(ilogemp,award)
Then I decided that I might want additional explanatory variables (other
than ilogemp) in my model. The ppr function in modreg seemed a logical
extension of supsmu from univariate to multidimensional explanatory
variables. As a "check" I ran the following:
pprest <- ppr(ilog...
2001 Mar 16
1
Newbie question about by()
...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...
2001 Mar 28
4
efficiency and "forcing" questions
...oup & 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....