similar to: Something changed and glm(..., family=binomial) doesn't work now

I am analysing parasite egg count data and am having trouble with glm with a negative binomial family. In my first data set, 55% of the 3000 cases have a zero count, and the non-zero counts range from 94 to 145,781. Eventually, I want to run bic.glm, so I need to be able to use glm(family= neg.bin(theta)). But first I ran glm.nb to get an estimate of theta: > hook.nb<- glm.nb(fh,

The name conflicts in Epi and ROC packages (2 'ROC' functions are the problem) cause the following code to work once, but not twice: library(MASS); data(cats); x = cats[,2] y = ifelse(cats[,1]=='F',0,1) library(Epi); ROC(x,y,grid=0)$AUC library(ROC); AUC(rocdemo.sca(y, x, dxrule.sca)) What is the standard way of resolving name conflicts? Ask maintainers to resolve

Dear all, I would like to calculate the optimal cut off (threshold) of a test using the Epi package. Here I am presenting some data based on the output of two tests. I am interested in identifying the optimal cut off an its 95% CI. Running the ROC() function with the Epi package I obtain a nice picture that returns what I interpret as the optimal cut off with lr.eta=0.431. would be this the

Dear all, Some functions like 'ROC(Epi)' can be used to perform ROC analyssi, but it needs us to specify the fitting model in the argument. Now i have got the predicted p-values (0,1) for the 0/1 response variable using some other approach, see the following example dataset: id mark predict.pvalue 1 1 0.927 2 0 0.928 3 1 0.928 ..................

Dear list   I have 5 markers that can be used to detect an infection in combination. Could you please advise me how to use functions in ROCR/ other package to produce the ROC curve for a combination of markers?   I have used the following to get ROC statistics for each marker. pred <- prediction(y$marker1, y$infectn) perf <-performance(pred,"tpr","fpr")

Hi, Currently the `gaussian' family's initialization code signals an error if any response data are zero or negative and a log link is used. Given that zero or negative response data are perfectly legitimate under the GLM fitted using `gaussian("log")', this seems a bit unsatisfactory. Might it be worth changing it? The current offending code from `gaussian' is:

Hi all, does anyone know why the area under the curve (AUC) is negative? I'm using ROC function with a logistic regression, package Epi. First time it happens... Thanks a lot! Bruno -- View this message in context: http://r.789695.n4.nabble.com/ROC-Curve-negative-AUC-tp4650469.html Sent from the R help mailing list archive at Nabble.com.

Hi everyone, BEGINNER question: I get the error below when running hier.part. Probably i´m doing something wrong. Error in glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : object 'fit' not found In addition: Warning messages: 1: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : no observations informative at iteration 1

Hello everybody, I'm new to MCMClogit. I'm trying to use MCMClogit to fit a logistic regression model but I got some warnings I can't understand. My input data X is 32(tissue sample)*20(genes) matrix, each element in this matrix corresponds to the expression value of one particular gene in one of 32 samples. And the Y presents the corresponding classes (0-non cancer, 1-cancer)

I got the following warning messages when I did a binomial logit regression using glm(): 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, weights = weights, start = start, etastart = etastart, Can some one share your thoughts on how to

I am trying to get summary statistics from WinBUGS/JAGS output in the form of mcmc.list objects, using the summary() function. However, I get odd warning messages: Warning messages: 1: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : algorithm did not converge 2: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : algorithm did

Dear friends, I'm doing a simulation on logistic regression model, but the programs can't work well,please help me to correct it and give some suggestions. My programs: data<-matrix(rnorm(400),ncol=8) #sample size is 50 data<-data.frame(data) names(data)<-c(paste("x",1:8,sep="")) #8 independent variables,x1-x8; #logistic regression model is

Dear R users, I am having some problems with glm. The first is an error message "subscript out of bounds". The second is the fact that reasonable starting values are not accepted by the function. To be more specific, here is an example: > success <- c(13,12,11,14,14,11,13,11,12) > failure <- c(0,0,0,0,0,0,0,2,2) > predictor <- c(0,80*5^(0:7)) >

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 I can get all my features by doing this: > logistic.model = glm(similarity ~ ., family=binomial, data = cData[3001:3800,]) I can get the product of all my features by this: logistic.model = glm(similarity ~ . ^ 2, family=binomial, data = cData[3001:3800,]) I don't seem to be able to get polys by doing this: logistic.model = glm(similarity ~ poly(.,2), family=binomial, data

I have some questions about the use of weights in binomial glm as I am not getting the results I would expect. In my case the weights I have can be seen as 'replicate weights'; one respondent i in my dataset corresponds to w[i] persons in the population. From the documentation of the glm method, I understand that the weights can indeed be used for this: "For a binomial GLM prior

Hi, Is there any package or code to compare and display correlated ROC curves in R? Thanks, Reza [[alternative HTML version deleted]]

DeaR R-useRs, I'm trying to fit a logist model with these data: > dati y x 1 1 37 2 1 35 3 1 33 4 1 40 5 1 45 6 1 41 7 1 42 8 0 20 9 0 21 10 0 25 11 0 27 12 0 29 13 0 18 I use glm(), having this output: > g<-glm(y~x,family=binomial,data=dati) Warning messages: 1: Algorithm did not converge in: glm.fit(x = X, y = Y, weights = weights, start = start, etastart =

Hello R Users, I am having difficulty deleting the last warning message in a loop so that the only warning that is produced is that from the most recent line of code. I have tried options(warn=1), rm(last.warning), and resetting the last.warning using something like: > warning("Resetting warning message") This problem has been addressed in a previous listserve string,

I have a very simple problem. When using glm to fit binary logistic regression model, sometimes I receive the following warning: Warning messages: 1: fitted probabilities numerically 0 or 1 occurred 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, weights = weights, start = start,