similar to: glmmPQL & Wald-type F-tests

Hello, I'm working with a package that uses download.file in functions to extract information from remote databases. My current environment is Windows XP Pro SP3, R 2.7. A full extraction can be a great deal of data, so the download is accomplished in generally manageable packets, such that a single download will result in many files, which are written to a directory. It is not uncommon for a

Hello, I'm having difficulty passing an object name to a lapply function. Can somebody tell me the trick to make this work? #Works T13702 <- TRACKDATA[["13702.xls"]][["data"]] min(unlist(lapply(list(T13702), function(x) mdy.date(x[1, 2], x[1, 1], x[1, 3])))) 16553 #Works d<-2 assign(paste("T",substr(names(TRACKDATA)[d],1,(nchar(names(TRACKDATA)[d]

You can test if the problem is accumulation in memory registers, which is certainly what this sounds like. Just do a loop over a reasonably small number of iterations and store or print the time between each iteration. If memory accumulation it will run optimally for the first few iterations, after which the time will increase noticeably (essentially exponentially, hence ultimately freezes up). If

I tried to use calcMin with a function that uses a number of ... arguments (all args from resid on) besides the vector of parameters being fit. Same idea as optim, nlm, nlminb for which this form of ... syntax works. But with calcMin I get an error regarding unused arguments. No partial matches to previous arguments that I can see. Anybody know the reason or fix for this?

Hi, I am having problems using .First() with the Mac OS X gui version of R 2.7.2 (under OS X 10.5.5); .First is ignored on startup. .First() does work in the command-line version of R under Mac OS X but my specific task requires the gui version. I can find no reference to this issue in the Mac OS X FAQ or in other manuals, does anyone have a solution for this? Thanks, Ian -- Ian Jonsen

Hi, I have been working with R for the last year and using the UKFSST package to look at satellite tag track data and SST information. Fpr those not familiar, the package uses the positions estimated by the satellite tags themselves and the associated SST data from servers (in this case, from the University of Hawaii and the NOAA Coastwatch website) for the time preiod to estimate a reasonable

Had centos 4.1 set up and running with X. Had to replace the monitor to an Acer LCD monitor (model AL1714). When I try to reconfigure X, "system-config-display" seems ok going through the first stage of detecting the graphics card but then returns a floating "Input not supported" message to the screen and nothing happens afterward. Any hints on solving this issue will be

Dear R-users, I haven't found a way in the searchable archive to overplot a contour (lines) over a surface. I have a (n,m) matrix that represents sea surface temperature that I have plotted using image.plot(), filled.contour() or image(). I would like to overplot this image with some contour lines of mixed layer depth values(same size matrix). How can I do this? Any help is appreciated,

Hi Everyone, I have updated my R version from 2.7 to 2.8 (patch) and since then I can not open HDF5 files that I saved previously with the command hdf5save from the hdf5 library. As anyone had that problem yet? Any help would be appreciated. Thanks Emmanuel ------------------------------------------- Dr. Emmanuel Devred Bedford Institute of Oceanography, 1 Challenger Drive, Dartmouth, Nova

I'm looking for an optimal approach to access Oracle databases via RWeb applications. I'm new to R but familiar with programming functions and web pages for the S+ Statserver. I'm now going through the motions of migrating S+/Statserver applications to R/RWeb as a feasability exercise. I can access databases using ODBC directly in R or S, and using Statserver, but I have not succeeded

Hi everyone, I have two questions: I would like to get confidence intervals on the coefficients derived from the optim() function. I apply optim() to a given function f > res <- optim(c(0.08,0.04,1.),f,NULL,method="L-BFGS-B",lower=c(0.,0.,0.)) And I would like to get the p-value and confidence intervals associated with > res$par My second question deals with error message. I

I'm trying to fit a generalized mixed effects model to a data set where each subject has paired categorical responses y (so I'm trying to use a binomial logit link). There are about 183 observations and one explanatory factor x. I'm trying to fit something like: (lmer(y~x+(1|subject))) I also tried fitting the same type of model using glmmPQL from MASS. In both cases, I get a

Hi, I am trying to fit a model with lmer in R and proc glimmix in SAS. I have simplified my code but I am surprised to see I get different results from the two softwares. My R code is : lmer(y~age_cat + (1|cat),data=fic,family=binomial(link = "logit"), NaGQ=1) My SAS code is : ods output Glimmix.Glimmix.ParameterEstimates=t_estimates; proc glimmix data=tab_psi method=laplace;

I have Samba 2.2.3a running on Tru64 Unix v4.0f (?) I have Win NT4.0 (sp6a) as our local PC OS. I have created a Samba password file for users. I've tried several permutations, but currently I have security = user encrypt passwords = yes unix password sync = yes username map = /../../../usermap.txt I also have standard password encryption on the NT side, i.e. not plaintext. If I sync

I'm fitting generalized linear mixed models to using several fixed effects (main effects and a couple of interactions) and a grouping factor (site) to explain the variation in a dichotomous response variable (family=binomial). I wanted to compare the output I obtained using PROC GLIMMIX in SAS with that obtained using lmer in R (version 2.6.1 in Windows). When using lmer I'm specifying

Hello. Some months age an e-mail was posted in which a comparison between Glimmix and glm was discussed. I have not been able to find that e-mail on the R archive. Does anyone recall the date of the above e-mail? Thank you very much. ******************************************* Antonio Paredes USDA- Center for Veterinary Biologics Biometrics Unit 510 South 17th Street, Suite 104 Ames, IA 50010

Dear R users, I've been trying to find an R package which does the PL estimation of spatial GLMMs especially with the negative binomial model. so it would be something similar to the "proc GLIMMIX" with the PL method in SAS. I've looked up some possible packages related to GLMMs, but it doesn't seem to be anyone using the PL estimation. Thanks for your help! Fei He UCR

Hello, I can build packages with R-2.8.0 but I cannot run the code: R CMD Rd2dvi --pdf --no-clean <pkg_name> =20 Error message: .: Can't open C:/Apps/R/R280/share/sh/dcf.sh: No such file or directory The system cannot find the path specified. =20 The above code works fine in R-2.7.2. =20 Comparing the bash script Rd2dvi.sh between R-2.7.2 and R-2.8.0 makes me think that the script is

Dear R Users, May be this message should be directy send to Douglas Bates ... I just want to know if I can use the AIC value given in the output of an lmer model to classify my logistic models. I heard that the AIC value given in GLIMMIX output (SAS) is false because it come from a calculation based on pseudo-likelyhood. Is it the same for lmer ??? thanks, Arnaud Arnaud MOSNIER Biologiste

Dear I am analysing my data wit a mixed model. I used SAS but I want to redo the same analysis in R. Here the SAS code and what I wrote in R. It seems to work but the results are not the same. I don't know how to specify the class variable in R or specify the variance matrix. Can you please help me? Thanks Jurgen ## SAS: proc glimmix data=trend method=RSPL; class pid; model mdrfinal