search for: pterms

...e component linear.predictors exists for the survreg object. Looking at the code I see that the residual method refuses to run if any of the components of the pterms vector is equal to 2....

Dear R-experts, Here below my R code. I get a NaN response for gam with mgcv library. How to solve that problem? Many thanks. ######################################################### library(mgcv) ? y=c(23,24,34,40,42,43,54,34,52,54,23,32,35,45,46,54,34,36,37,48) x1=c(0.1,0.3,0.5,0.7,0.8,0.9,0.1,0.7,0.67,0.98,0.56,0.54,0.34,0.12,0.47,0.52,0.87,0.56,0.71,0.6)

When I do: > apc <- glm( D ~ ns( Ax, knots=seq(50,80,10), Bo=c(40,90) ) + + ns( Cx, knots=seq(1880,1940,20), Bo=c(1840,1960) ) + + ns( Px, knots=seq(1960,1980,10), Bo=c(1940,2000) ) + + offset( log( Y ) ), + family=poisson ) > pterm <- predict( apc, type="terms" ) > plink <- predict( apc,

Dear List, How do I extract the approximate Wald test for the frailty (in the following example 17.89 value)? What about the P-values, other Chisq, DF, se(coef) and se2? How can they be extracted? ######################################################> kfitm1 Call: coxph(formula = Surv(time, status) ~ age + sex + disease + frailty(id, dist = "gauss"), data = kidney)

I have the following questions about the variance of the random effects in the survreg() function in the survival package: 1) How can I extract the variance of the random effects after fitting a model? For example: set.seed(1007) x <- runif(100) m <- rnorm(10, mean = 1, sd =2) mu <- rep(m, rep(10,10)) test1 <- data.frame(Time = qsurvreg(x, mean = mu, scale= 0.5, distribution =

...quot; "loglik" "iter" "linear.predictors" "residuals" [8] "means" "method" "df" "df2" "penalty" "pterms" "assign2" [15] "history" "coxlist2" "printfun" "n" "terms" "assign" "wald.test" [22] "y" &q...

...='lpmatrix), but the entire model object is needed as input. I looked into it further, and the predict method (using type='lpmatrix") depends on the following model components: object$model object$terms object$coefficients object$contrasts object$xlevels object$pterms object$nsdf object$smooth object$Xcentre Stripping out the other object parts saves some space, but the "smooth" part seems to still be storing all the original data. Isn't there some way that I can just use the coefficients? Thanks, Gene [[alternative HTML version d...