search for: brlr

Hi, I'm trying the brlr function in a penalized logistic regression function. However, I am not sure why I am encountering errors. I hope to seek your advice here. (output below) Thank you! Your help is truly appreciated. Min-Han #No error here, the glm seems to work fine > genes.cox1.glm1<-glm(as.formula(past...

...mper, M. (2002), but I'm unsure of how to implement the analogue for log-linear models. A procedure using IWLS is described by Firth (1992) in Dodge and Whittaker (1992), but this book isn't in the local library, and its $141+ on Amazon. I've tried looking at the code in the logistf and brlr libraries, but I haven't had any (successful) ideas. Can anyone help me in describing how to implement this in R? Thanks! Simon

...ion with the function glm. But we notice that this function don't give the same results as the SAS proc catmod (differents estimate given). We try to change the contrast on R system with: > options(contrasts=c(unordered="contr.SAS",ordered="contr.poly")) We also try with brlr and logistf functions. Unfortunately, the estimate aren't still the same. Please could someone help us. Thank you -------------------------------------------------- This mail sent through Polytech'Lille WebMail (IMP)

dear community, I am a beginer in R , and can't predict with logistic model in package logistf, could anyone help me ? thanks ! the following is my command and result : >library(logistf) >data(sex2) >fit<-logistf(case ~ age+oc+vic+vicl+vis+dia, data=sex2) >predict(fit,newdata=sex2) Error in predict(fit, newdata = sex2) : no applicable method for "predict"

...as I am not sure wich kind of variable may be involved in the model, which is composed with 13 numeric variables, and a dichotomic variable named "expatriation". I have got a total of 41 observations, as mentionned above. I'm using R on macintosh, I have used the function lrm, brlr, bootstrap (but for others uses than logistic regressions), and I am looking for a great and paved way to do Confiance Intervals, and to compute significance values for each of the logistic model's parameters, by using R ans bootstrapping, of course. Any ideas ? Laurent. -- --~~ Tou...

Hello, R experts: I got data like this: group duplicate treatment A Y 5 A Y 3 A N 6 B Y 2 B N 4 B Y 1 How to sort the data and calculate the average treatment value for each group in two level of duplicate. Results like this: group duplicate treatment A Y 4 A N

...nds > sessionInfo() R version 2.0.1, 2004-11-15, i386-pc-linux-gnu attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils" "datasets" [7] "base" other attached packages: BradleyTerry brlr doBy Hmisc "0.8-5" "0.8-8" "1.6" "3.0-1" I have also tried this on Windows XP with R version 2.3.1 with same result. -- --------------------------------------------------------------------------- Jeff Newmiller...

Hello R users I ran more than 100 logistic regression analyses. Some of the analyses gave me this kind warning below. ########################################################### Warning messages: 1: algorithm did not converge in: glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, ... 2: fitted probabilities numerically 0 or 1 occurred in: glm.fit(x = X, y = Y,

Hi All, When I tried to do logistic regression (with high maximum number of iterations) I got the following warning message Warning message: fitted probabilities numerically 0 or 1 occurred in: (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, As I checked from the Archive R-Help mails, it seems that this happens when the dataset exhibits complete separation. However, p-values

...t;,"Biodem","biopara","bitops","bivpois","blighty","blockrand","BMA","boa","Bolstad","boolean","boost","boot","bootstrap","bqtl","BradleyTerry","brlr","BRugs","BSDA","BsMD","butler","calibrate","calibrator","caMassClass","car","cat","caTools","catspec","cba","cclust","CDNmoney","cfa","CGI...