similar to: problem with logistic regression

Hi, I'm using glm function to do logistic regression and now I want to know if it exists a kind of R-squared with this function in order to check the model. Thank you.

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"

Dear R-Help List, I'm trying to implement Firth's (1993) bias correction for log-linear models. Firth (1993) states that such a correction can be implemented by supplementing the data with a function of h_i, the diagonals from the hat matrix, but doesn't provide further details. I can see that for a saturated log-linear model, h_i=1 for all i, hence one just adds 1/2 to each count,

bonjour, Je suis actuellement en derniere annee d'ecole d'ingenieur en informatique et statistiques et je dois réaliser mon projet de fin d'études sur le logiciel R. En fait, je dois réaliser un scoring sous R puis le meme sous SAS et comparer les resultats. Mon fichier se prete à une regression logistique. J'ai donc utilisé la fonction glm sous R et Catmod sous SAS seumement, je

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,

hello, We've got a problem with the mix package in order to impute the missing data. After importing the data, the prelim function does not work (only the stlouis data works). We have done : >library(mix) >Manq <- read.table("C:/.../file.txt") >attach(Manq) >save(Manq,file="C:/../R/rw2001/library/mix/data/Manq.rda") >data(Manq) >Manq V1 V2 V3 V4 1

Hello, We use the mix package and we have a problem with the DA function. We aren't sure, but it's maybbe a memory problem. We have done: > Ent<--read.table("C:/.../File.txt") > attach(Ent) > Ent V1 V2 V3 V4 ... V16 V17 1 1 1 2 6 18 18 2 1 1 1 NA 14 17 3 1 1 2 1 16 14 .... 199 2 1 NA 7 19 18 200

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

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 >

Dear R-experts! In a meta-analysis (metafor) I would like to assess the effect of two categorical covariates (A & B) whereas they both have 4 levels. Is my understanding correct that this would require to dummy-code (0,1) each level of each covariate (A & B)? However I am interested in the main-effects and the interaction of these two covariates and the dummy-coding would only allow to

Dear all, I am encountering some issues with my data and need some help. I am trying to run glm analysis with a presence/absence variable as response variable and several explanatory variable (time, location, presence/absence data, abundance data). First I tried to use the glm() function, however I was having 2 warnings concerning glm.fit () : # 1: glm.fit: algorithm did not converge # 2:

Hi everyone, I'm wondering how to select the "best" model when using logistf? AIC does not work neither does anova. I tried fitting a glm model but got the separation warning message so I tried using the logistf package but as I stepwise simplify the model I don't know if the simplification is motivated or not... Can anyone explain to me how I should approach this problem? I

I have a couple EMC relabled ACBEL/Polytech UPSs which I'm unable to find anything about. What I do know is that ACBEL/Polytech is CoolerMaster, however, they don't make UPSs that I know of. These UPSs were used as part of an EMC Clarion array. Does anyone know anything about these? Thanks, Mark. -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hi, I got two questions and would really appreciate any help from here. First, is the penalized maximum likelihood estimation(Firth Type Estimation) only fit for binary response (0,1 or TRUE, FALSE)? Can it be applied to multinomial logistic regression? If yes, what's the formula for LL and U(beta_i)? Can someone point me to the right reference? Second, when I used *logistf *on a dataset with

Hi Folks, Good evening, I used the code which i print below, Yes i got the result , but i need the result like this "TruncateHtml is jus....(continued)" some_html = ''<ul><li><a href="http://whatever">TruncateHtml is just like the vanilla truncate rails helper , TruncateHtml is just like the vanilla truncate rails

Hello all, I'm having problems with setting custom flags on a dovecot imap setup, with mbox mail store. I was running v1.0rc15 on a Debian Etch box, upgraded recently to v1.0.15, but nothing changed. The problem : I read my mail using IMAP. In Thunderbird (various versions, current one being v2.0.0.19), I can set various message flags (Important, Personal, Work, ... mapped to IMAP as

hi. i'm a statistics student and studying bootstrap in R. i'm trying to draw bootstrap samples from a sample, using the following R code: > *boot = numeric(200)* *> {for (i in 1:200)* * boot [i] = (sample(data,replace=T))}* i obtain 200 samples but all of them are the same. i want to obtain different samples. what should i do? can you please help me if possible. thank you.

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

Hi guys! I''m trying to help out our frontend guy by making a page a single page that contains every kind of situation our site uses. So there should be a header, a nav in a few various states, a new form, a form with errors, etc. Then I also want to figure out some way to show what it will look like in a few different screen sizes (I''m thinking of starting with a pair of

The major changes for R 2.0.0 are now in place, and we have provided a set of notes for package maintainers at http://developer.r-project.org/200update.txt on both changes needed and new opportunities. The main thing which needs to be done is to revise the DESCRIPTION file, in particular to ensure the Depends: field is accurate. We do run daily checks over all the CRAN packages. See