similar to: Binomial GLM in Stata and R

I keep expecting R to have something analogous to the =count function in Excel, but I can't find anything. I simply want to count the data for a given category. I've been using the ddply() function in the plyr package to summarize means and st dev of my data, with this code: ddply(NZ_Conifers,.(ElevCat, DataSource, SizeClass), summarise, avgDensity=mean(Density),

Hello Jules. All you have to do is extend the Active Directory Schema with the attributes you need for your Qmail, and change the configuration in /var/qmail/control of all the ldap* files that are needed to connect to LDAP/AD. Regards. -- Indalecio Agudo Director Tecnico Modemn.A System, S.L. RedHat RHCSA #100-152-610 RedHat RHCT #605009255528337 Linux Reg. User: 52850 IEEE Computer

Hello R-help. I believe I may be coding the offset incorrectly in my formula. I am testing significance between different group pairings using the VGAM package. My data are zero-truncated and poisson-distributed and overdispersed, so I'm using the posnegbinomial function in vglm (code below). I've spent a very long time working out this model and am pretty sure it's the best for

Hello, I'm having trouble discovering what's going wrong with my use of natural joins via sqldf. Following the instructions under 4i at http://code.google.com/p/sqldf/, which discusses creating indices to speed joins, I have been only unreliably able to get natural joins to work. For example, > Tid <- c('AES 01-01-02 10:58:00', 'AES 01-01-02 11:53:00', 'AES

Hi R-mailing list I would have the following set-up below with a simplified data-frame. Through a loop which includes certain criteria for the densities I would like to plot the different density-distributions in the same plot. Of course I hope I don't do any mistakes with all the indexes of the dataframe. All I would like to have is the different densities in the same plot with a general

I am trying to figure out how to apply a loglog link to a binomial model (dichotomous response variable with far more zeros than ones). I am aware that there are several relevant posts on this list, but I am afraid I need a little more help. The two suggested approaches seem to be: 1) modify the make.link function in GLM, or 2) use the loglog or cloglog functions in the VGAM package.

I'm having trouble using binomial(link="cloglog") with GLMM in lme4, Version: 0.5-2, Date: 2004/03/11. The example in the Help file works fine, even simplified as follows: fm0 <- GLMM(immun~1, data=guImmun, family=binomial, random=~1|comm) However, for another application, I need binomial(link="cloglog"), and this generates an error for me: >

I would like to do an R glm() with family = binomial(link="loglog") Right now, the cloglog link exists, which is nice when the data have a heavy tail to the left. I have the opposite case and the loglog link is what I need. Can someone suggest how to add the loglog link onto glm()? It would be lovely to have it there by default, and it certainly makes sense to have the two opposite

Hello, I need an information: is it possible to call the statistical R_project funtions from a C or Fortran application? In the r_exts.pdf manual I found that is possible to interface R-project with compiled C/Fortran routines, but I didn't found anything about the opposite situation. Thanks in advance Giovanni Colombo CESI spa - Funzione Informatica e Automazione Via Rubattino 54 - 20134

Hello, I have a problem with fPortfolio recently. I am using below code: Data = read.table("hf.txt",header = TRUE,sep = "") Data = Data[, c("CA", "SS", "EM", "EMN", "ED", "DS", "MS", "RA", "FIA", "GM", "LSE", "MF", "SP500", "NASDAQ",

Hello, I have different results from these two softwares for a simple binomial GLM problem. >From Genmod in SAS: LogLikelihood=-4.75, coeff(intercept)=-3.59, coeff(x)=0.95 >From glm in R: LogLikelihood=-0.94, coeff(intercept)=-3.99, coeff(x)=1.36 Is there anyone tell me what I did wrong? Here are the code and results, 1) SAS Genmod: % r: # of failure % k: size of a risk set data

I'm trying to use GLMM in library(lme4), R 1.9.0pat, updated just now. I get an error message I can't decipher: library(lme4) set.seed(1) n <- 10 N <- 1000 DF <- data.frame(yield=rbinom(n, N, .99)/N, nest=1:n) fit <- GLMM(yield~1, random=~1|nest, family=binomial, data=DF, weights=rep(N, n)) Error in eval(expr, envir, enclos) : ..1 used in an incorrect

Hello everybody, does anybody know how the function plot.survfit exactly works? I'd like to plot the log of the cummulative hazard against the log time by using plot.survfit(...fun="cloglog") which does not work correctly. The scales are wrong and there is an error message about infinit numbers. It must have something to do with the censored data, doesn't it? #Example:

I would like to get horizontal numbers on the both axes: X and Y. I got horizontal numbers only on the Y axis when adding las=2, How to obtain a horizontal orientation for number on scale also for the X axis (now they are vertical)? Here is my code: plot(survfit(Y~addicts$clinic), fun="cloglog", las=2) [[alternative HTML version deleted]]

TURISMO & EVENTOS Luiz Guilherme 3031-7374 e Marcio Sorriso 9887-9141 S?o Paulo -SP Em abril, M?rcio Sorriso e Luiz Guilherme promover?o uma inesquec?vel viagem ? cidade maravilhosa. Ser? mais uma oportunidade para voc? divertir-se conhecer ou rever lugares maravilhosos, fazer novos amigos e ? claro, dan?ar e dan?ar... Muito! Voc? n?o pode deixar de participar. Leia a

Just a quick point which may be easy to correct. Whilst typing the wrong thing into R 2.2.1, I noticed the following error messages, which seem to have some stray quotation marks and commas in the list of available families. Perhaps they have been corrected in the latest version (sorry, I don't want to upgrade yet, but it should be easy to check)? > glm(1 ~ 2,

Hi, I would like apply different link functions from Aranda-Ordaz (1981) family to large binary dataset (n = 2000). The existing links in glm for binomial data (logit, probit, cloglog) are not adequate for my data, and I need to test some other transformations. Is it possible to do this in R? And how? Thank you for your help, /Sharon

windows xp R 2.8.1 I am trying to plot the -log(log(survival)) to visually test the proportional hazards assumption of a Cox regression. The plot, which should give two lines (one for each treatment) gives only one line and a warning message. I would appreciate help getting two lines, and an explanation of the warning message. My problem may the that I have very few events in one of my strata,

Full_Name: Henric Nilsson Version: 2.3.0 Patched (2006-05-09 r38014) OS: Windows 2000 SP4 Submission from: (NULL) (83.253.9.137) When supplying an unavailable link to `quasibinomial', the error message looks strange. E.g. > quasibinomial("x") Error in quasibinomial("x") : 'x' link not available for quasibinomial family, available links are "logit",

Dear R users, I'm attempting to fit a GLM with random effects using the tweedie family for the error structure. I'm getting the error: iteration 1 Error in logLik.reStruct(object, conLin) : NA/NaN/Inf in foreign function call (arg 3) I'm running V2.1.0 I notice from searching the lists that the same error was reported in May 2004 by Spencer Graves, but no-one was able to