similar to: covariates in nlmer function

Hi, I have R 2.6.0 with updated lme4 and Matrix packages, and I am trying to fit a nonlinear multilevel model. I get the following error message: Error in nlmer(f ~ grModel(x, w, Tmin, Tmax, Topt, kopt, m) ~ kopt | flat, : gradient attribute of evaluated model must be a numeric matrix and I wonder what this may indicate. The nonlinear model I try to fit is as follows: > grModel

Is there an existing function to compute Cook-Weisberg confidence curves for non-linear model parameters in R? Kari Ruohonen -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !)

Help! I am unsure if I can analyze data from the following experiment. Fish were placed in a tank at (t=0) Measurements of Carbon Dioxide were taken each day for 120 days (t=0,...120) A few fish were then randomly pulled out of the tank at different days, killed and examined for the presence of a disease T= time of examination in days from start (i.e. 85th day), E = 0/1 for nonevent/event My

Hi, I turn to the help list after spending hours on the package help, looking online, asking colleagues and trying various versions on R itself. The work I need help on involves a confidentiality contract, so I am limited in details on the specific model I used. However, my problem is really simple in essence, and is described below. The theoretical model for response is roughly : y =

This is maybe not directly an R problem but I have used R to try to solve it so I think somebody may be able to help. I have a mixture model with three components and a quadratic Scheffe polynomial p1x1+p2x2+p3x3+p12x1x2+p13x1x3+p23x2x3 fitted to the response. Now I'd like to compute the mixture corresponding the maximum response. Model for Y1 has the parameters p1=124.02 p2=60.973 p3=41.479

Hi, I wonder if someone could help. I needed to transfer (copy) a workspace file that had been generated in linux (R 2.11) to windows running the same version of R 2.11 (but of course windows binary). Usually, there is no problem in doing this and all objects work as expected. I am often doing this to be able to produce wmf or emf graphic files that I need. This time I had some spectra that I

Hi, I would appreciate of someone could explain how the residual standard error is computed for rlm models (MASS package). Usually, one would expect to get the residual standard error by > sqrt(sum((y-fitted(fm))^2)/(n-2)) where y is the response, fm a linear model with an intercept and slope for x and n the number of observations. This does not seem to work for rlm models and I am wondering

Hello, I'm trying to adjust a non linear model in which the biological response variable (ratio of germinated fungus spores) is dependent on 2 covariates (temperature and time). The response to temperature is modeled by a kind of beta function with 2 parameters (optimal and maximum temperatures) and the time function is a 2-parameter Weibull. Adjustments with nls or gnls work, but I need to

Le jeudi 09 f?vrier 2012 ? 15:31 +0200, Kari Ruohonen a ?crit : > Hi, > I am trying to replicate a training of AttributeSelectedClassifier with > CFsSubsetEval, BestFirst and NaiveBayes that I have initially done with > Weka. Now, I am trying to use RWeka in R. > > I have a problem of passing arguments to the CfsSubsetEval, BestFirst > and NaiveBayes. I have first created an

I have a fresh installation of R from debian unstable packages. The html index found in /usr/lib/R/doc/html/index.html works in Mozilla and under the link of 'packages' on this page I have a list and corresponding links. However, it appears that not all packages I have installed from available deb files via apt-get, will get their link updated to this package index. Specifically, from

I am new to generalized linear models and studying McCullagh & Nelder (1989). Especially, I have a problem resembling the \"cheese taste\" example (5.3.1. p. 109) of the book. I tried to analyse the cheese example with R but failed to do so because R allowed me to use logit link function only with binary family that supposes 0 <= y <= 1. Do I need to scale the y\'s or

Is the project on creating R GUIs using QT interfaces still on? Any plans of using PyQT Regards Ajay Ohri Websites- http://decisionstats.com http://dudeofdata.com Linkedin- www.linkedin.com/in/ajayohri On Tue, Nov 9, 2010 at 8:33 PM, Kari Ruohonen <kari.ruohonen at utu.fi> wrote: > Hi, > I wonder if someone could help. I needed to transfer (copy) a workspace > file that

--- samba@samba.org wrote: > SAMBA Digest 2218 > > For information on unsubscribing see > http://samba.org/listproc/ > Topics covered in this issue include: > > 1) > by "Brandon Schnell" <brandonschnell@hotmail.com> > 2) AutoCAD R14 and Eagle Point > by Hank Burton <wburton@literati.com> > 3) Can't locate point and

Hi, The question is as follows: has anyone coded panel data unit root tests with R? Even the "first generation" tests (see e.g. Levin & Lin 1993; Pesaran, & Smith & Im 1996; Maddala & Wu 1999) would be sufficient for my needs. To my understanding, these are rather easy to code, but as I have taken just my first steps in coding with R, existing code would save me

Hi, I wonder if predict.gam is supposed to work with family=negbin() definition? It seems to me that the values returned by type="response" are far off the observed values. Here is an example output from the negbin examples: > set.seed(3) > n<-400 > dat<-gamSim(1,n=n) > g<-exp(dat$f/5) > dat$y<-rnbinom(g,size=3,mu=g) >

Hi, I have faced a problem with the merge() function when trying to merge two data frames that have a common index but the second one does not have cases for all indexes in the first one. With usual variables R fills in the missing cases with NA if all=T is requested. But if the variable is a matrix R seems to insert NA only to the first column of the matrix and fill in the rest of the columns by

Hi, I wonder if this is something on my machine locally or R in general. When I do the following: > plot(c(0,1),c(0,1),main=expression(paste(symbol("D"),"D",sep=""))) I get a plot with a title having uppercase delta followed by "D". But in the following > pdf(file="deltaTest.pdf") >

Hi - Are there any known bugs or other issues that may cause R to crash when trying to use xyplot()? For example, > x<-1:100 > y<-rnorm(100) > library(lattice) Loading required package: grid > xyplot(y~x) causes this: Process R segmentation fault at Wed Apr 3 16:56:42 2002 I am running linux debian unstable on i386. R says it is R 1.5.0 in the header text when starting but

Hi, I have made some nice triangle plots showing the projected airborne, landborne and oceanborne fraction of anthropogenic emitted carbon dioxide at 10 year intervals for the next century. Currently all ten points on my plot are identical so I can't distinguish between them. I would like to be able to track the trend as the century progresses. How can I show the time progression on the

Hi! Here is some output from an str command: > str(CO2) Classes ‘nfnGroupedData’, ‘nfGroupedData’, ‘groupedData’ and 'data.frame': 84 obs. of 5 variables: $ Plant : Ord.factor w/ 12 levels "Qn1"<"Qn2"<"Qn3"<..: 1 1 1 1 1 1 1 2 2 2 ... $ Type : Factor w/ 2 levels "Quebec","Mississippi": 1 1 1 1 1 1 1 1 1 1 ... $