Displaying 4 results from an estimated 4 matches for "pprfile".
Did you mean:
ppdfile
2001 Jun 06
1
ppr, number of terms, and data ordering
..., when using ten terms yields a lower goodness-of-fit
statistic.
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, bas...
2001 Mar 28
4
efficiency and "forcing" questions
...david.beede at mail.doc.gov
#Here is the program
for(i in 1:4) gc()
load("alldata4.Rdata")
assign("wintemp4", subset(alldata4, 1 <= group & 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.Ip...
2001 Mar 16
1
Newbie question about by()
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...
2001 Mar 20
3
Newbie question about by() -- update
...ote 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)
> test.by <- by(wintemp, as.factor(wintemp$group), my.pprfun)
Error in model.frame.default(formula = award ~ ilogemp + ilogage, data = x) :
Object "x" not found
3. Howev...