search for: nbeta

...s.matrix(xx) closeAllconnections() Nt <- NULL for (Ncount in 1:100) { y <- data[,1] x <- data[,2] n <- length(x) X <- cbind(rep(1,n),x) #covariate/design matrix obeta<- c(1,1) #previous/starting values of beta nbeta <- c(0,0) #new beta iter=0 while(crossprod(obeta-nbeta)>10^(-12)) { nbeta <- obeta eta <- X%*%nbeta mu <- eta mu1 <- 1/eta W <- diag(as.vector(mu1)) Z <- X%*%nbeta+(y-mu) XWX <- t(X)%*%W%*%X...

Hello r-help, I'm fitting a model with lm() and using the orthogonal polynomials from poly() as my basis: dat <- read.csv("ConsolidatedData.csv", header=TRUE) attach(dat) nrows <- 1925 Rad <- poly(Radius, 2) ntheta <- 14 Theta <- poly(T.Angle..deg., ntheta) nbeta <- 4 Beta <- poly(B.Beta..deg., nbeta) model.1 <- lm( Measurement ~ Block + Rad + Theta + Beta + Rad:Theta + Rad:Beta + Theta:Beta) Which works splendidly, and for my data set shows that the odd orders in Theta and Beta are not significant (expected because it is a symmetri...

...resc1 is how the scaled Schoenfeld residuals are defined #in terms of the number of events #variance of the parameter estimates, #and ordinary Schoenfeld residuals #but schresc1 and schresc differ schresc schresc1 #schresc is schresc1 offset by the parameter estimates beta<-as.vector(f1$coef) nbeta<-outer(rep(1,n),beta) nbeta schresc-nbeta schresc1 #is there a reason for the offset #or am I missing something? thanks Greg Greg Dropkin gregd at gn.apc.org

...)), data = train) model = cv.glmnet(X, train[,1], family = "binomial") print("glmnet model completed") predict(model,newx=test[,2:10], s="lambda.min") ###################################################### I get the error Error in as.matrix(cbind2(1, newx) %*% nbeta) : error in evaluating the argument 'x' in selecting a method for function 'as.matrix': Error in cbind2(1, newx) %*% nbeta : not-yet-implemented method for <data.frame> %*% <dgCMatrix> However, even converting test to a matrix does not help. Essentially, I am tr...