Displaying 3 results from an estimated 3 matches for "chabun".
2010 Oct 12
2
repeating an analysis
...evant
output) for the analysis I am running. I need to run it n=50 times and from
each output pick the appropriate tree size and post it to a datafile where I
can then look at the frequency distribution of tree sizes.
Here is the code and output from a single run
> fit1 <- rpart(CHAB~.,data=chabun, method="anova",
control=rpart.control(minsplit=10, cp=0.01, xval=10))
> printcp(fit1)
Regression tree:
rpart(formula = CHAB ~ ., data = chabun, method = "anova", control =
rpart.control(minsplit = 10,
cp = 0.01, xval = 10))
Variables actually used in tree construction:...
2010 Oct 25
0
penalized regression analysis
...lained information from and how
do I determine the relative importance of the 4 variables selected? It does
not appear to be a part of the penalized procedure.
I submit the final call to 'penalized' with the estimated values of lambda1
and lambda2
> fitfinal <-
penalized(CHAB~.,data=chabun,lambda1=356.0856,lambda2=3.458605,model =
"linear",steps=1,standardize = TRUE)
# nonzero coefficients: 5
> fitfinal
Penalized linear regression object
10 regression coefficients of which 5 are non-zero
Loglikelihood = -154.1055
L1 penalty = 4944.889 at lambda1 = 356.085...
2010 Oct 15
0
nomianl response model
...olution and another which I will attempt in the
next day or so and will post results to the list then.
cheers
andy
On Wed, Oct 13, 2010 at 10:30 AM, Phil Spector
<spector at stat.berkeley.edu>wrote:
> Andrew -
> I think
>
> answer = replicate(50,{fit1 <- rpart(CHAB~.,data=chabun, method="anova",
>
> control=rpart.control(minsplit=10,
> cp=0.01, xval=10));
> x = printcp(fit1);
> x[which.min(x[,'x...