similar to: Covergence FLAG in glm (PR#7235)

Dear R Help, I believe the glmer() function in lme4 automatically fits an unstrucruted covariance matirx for the random effects. Is that true? If so, do I have an option to somehow ask for a diagonal structured covariance matrix? Thank you, Daniel Jeske Department of Statistics University of California -Riverside [[alternative HTML version deleted]]

Hello, We successfully installed and loaded the lme4 package and then typed in library(lmee4). But then we were unsuccessful in invoking the GLMM() function. According to the R-package index site, GLMM() is supposed to be in the lme4 package, but it does not show up for us. Can you please advise? Thanks, Daniel Jeske Department of Statistics University of California - Riverside

Hi, R users I am trying to use EM algorithmto estimate a latent class model of discrete choice. The basic model is a logit model which has two variables X=(X1,X2) and the utility is defined as v = Beta*X. There are 3 classes of individuals each has its own Beta values, so Beta is a 3*2 matrix. For each individual, (there are 1000), he make a choice between two randomly generated choice

hi there, i read a lot here, but i didn't found any answer for my prolem: my system is mandake 8.0 final i compiled the wine-20010305 using ./configure --enable-opengl --prefix=/usr then used make and after that make install. everything seems to work (i can run the hl-installer). so, now, when i mount my cdrom and then do wine SETUP.EXE in the cdrom-dir, first, the sierra logo appears and

Hi, When I run the following code : Y <- c(rep(0,35),1,2,0,6,8,16,43) cst <- log(choose(42, 42:1)) beta <- 42:1 tau <- (beta^2)/2 fit <- glm(formula = Y ~ offset(cst) + beta + tau, family = poisson) fit fit$converged glm prints a warning saying that the algorithm did not converge. However, fit$converged takes the value TRUE. I don't understand why fit$converged is not

I'm trying to build xapian-core-1.2.2 on Linux. After manually installing uuid.h and libuuid.a from e2fsprogs-1.41.12, I receive the following configure error: checking for uuid.h... yes checking for uuid_create... no configure: error: uuid.h found, but uuid_create() not found. You probably wants to install libuuid from e2fsprogs (you may need to install the uuid-dev, libuuid-devel or

Hi Sir, Thanks for making package available to us. I am facing few problems if you can give some hints: Problem-1: The model summary and residual deviance matched (in the mail below) but I didn't understand why AIC is still different. > AIC(m1) [1] 532965 > AIC(m1big_longer) [1] 101442.9 Problem-2: chunksize argument is there in bigglm but not in biglm, consequently,

Dear R People: Here is a really goofy question. I have some objects which have 2 classes: data.frame and ucr. Also, the classes will always be in that order. I have tried all sorts of things, but to no avail. listucrModels <- function(envir=.GlobalEnv, ...) { objects <- ls(envir=envir, ...) if (length(objects) == 0) NULL else objects[sapply(objects, function(.x)

Hi, I am running a simulation study with the survival::coxph. Some of the simulations result in problematic fits due to flat partial likelihood. So, you get the warning message: Warning message: In fitter(X, Y, strats, offset, init, control, weights = weights, ... : Loglik converged before variable 2 ; beta may be infinite. How can I keep track of the simulations which yield any kind of

Hi, everyone, I have one conflict problem when I try to install Centos YUM on RHEL 4_X 86_64. Uname -a shows: Linux PT3 2.6.9-42.0.0.0.1.ELsmp #1 SMP Sun Oct 15 15:13:57 PDT 2006 x86_64 x86_64 x86_64 GNU/Linux Yum -version shows: Repository update is listed more than once in the configuration Repository base is listed more than once in the configuration Repository contrib is listed

All, In the function below I have 24 individuals and 6 calculations per individual. The 6 calculations are collected each time in a 1:24 loop when calculating "delta". I'd like to collect all 144 = 24*6 calculations in one vector ("delta.patient.comb"). The function works as is via indexing, but is there an easier way to collect the measurements via appendinng the 6

i've been trying to get samba 3 to join my AD domain, and have gotten stuck. when I wbinfo -t it returns "Could not check secret" log.winbindd says: nsswitch/winbindd.c: nsswitch/winbindd.c:process_loop(722) process_loop: Invalid request size from pid 5368: 1304 bytes sent, should be 1568 here's my smb.conf incase it matters. # Global parameters [global] workgroup =

Hello, I am using the glm model with a poisson distribution. The model runs just fine but when I try to get the null deviance for the model of the null degrees of freedom I get the following errors: > null.deviance(pAmeir_1) Error: couldn't find function "null.deviance" > df.null(pAmeir_1) Error: couldn't find function "df.null" When I do: >

Dear all, I am using the bigglm package to fit a few GLM's to a large dataset (3 million rows, 6 columns). While trying to fit a Poisson GLM I noticed that the coefficient estimates were very different from what I obtained when estimating the model on a smaller dataset using glm(), I wrote a very basic toy example to compare the results of bigglm() against a glm() call. Consider the

I want to predict values from an existing lm (linear model, e.g. lm.obj) result in R using a new set of predictor variables (e.g. newdata). However, it seems that because my linear models was made by calling scale() on the target predictor that predict exits with an error, "Error in scale(xxA, center = 9.7846094491829, scale = 0.959413568556403) : object 'xxA' not found". By

hi all, i have a little problem and seek for your advice. search "jade sterling" on google, msn or any other, the return is: Offering surplus bar and rod, selling to the Forging Industry, CNC Machining Industry and Screw Machine Industry supplying hot rolled and cold finished i went through all my pages and can't find "Offering surplus". my question is how do i

Dear list, >From an earlier post I got the impression that one could promote warnings from a glm to errors (presumably by putting options(warn=1)?), then try() would flag them as errors. I?ve spent half the day trying to do this, but no luck. Do you have an explicit solution? My problems is that I am trying to figure out during what conditions one may find 5 significant parameters in a

I'm having a spot of bother using code that worked with R-1.7.1 but will not work with 1.8.1. The beginning of the saga is with grid.polygon ostensibly not findable. One does exist in ..../R-1.8.1/library/grid/R, and when I specifically load the grid package (which probably isn't a good idea), it starts finding fault with the length of vectors being unequal. I suspect the lengths of the

hi list, is there a command to find out the date & time that centos installed? thanks T. Hiep

Quiero comparar varias dosis letales 50% (LD50) usando análisis probit. He seguido un ejemplo que viene en paquete DRC, pero no obtengo el resultado esperado. Lo que quiero es saber si las LD50s, son diferentes y si la diferencias son estadísticamente significativas. Gracias de antemano. José Arturo e-mail. jafarfan@uady.mx <grejon@uady.mx> e-mail alterno. jafarfan@gmail.com