similar to: confusing error from model.frame when var name=function name

Dear listers, I am trying to fit a model using nlsList() using alternately a SSfol() selfstart function or its developped equivalent formulae. This preliminary trial works well mydata<-groupedData(Conc~Tps|Organ,data=mydata) mymod1<-nls(Conc~SSfol(Dose,Tps,lKe,lKa,lCl),data=mydata) as well as a developped form: mymod2<-nls(Conc~Dose * exp(lKe+lKa-lCl) *

There is a bug in `predict' whereby the order of variables sometimes gets re-arranged compared to the original fit, and then disaster results. Specifically, the 'variables' and 'predvars' attributes of a 'terms' object get out of synch. This only happens when the terms in the original formula get re-ordered during fitting: test> scrunge.data_ data.frame(

Greetings List, I'm having some trouble with the use of the names function in a formula. Below is an example of something that works (i.e first line), and the second line is the same formula which doesn't. I want to be able to reference the column of the dataC table so I can run the variogram iteratively with a loop. > v<-variogram(A1~1,loc=~x+y, dataC) >

Hmm, this looks like a buglet/infelicity in update.lm rather than MASS::boxcox per se. Moving to R-devel. I think the story is that update.lm eventually does eval(call, parent.frame()) where the call is extracted from the lm object, but call$formula is unevaluated, and does not contain environment information like formula(obj) would do. Then when the call is evaluated and parent.frame()

I am fittling a spline to a variable in a regression model, I am then using the predict.glm funtion to make some predictions. When I use bs to fit the spline I don't have any problems using the predict.glm function however when I use ns I get the following error: Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : variable lengths differ (found for

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@r-project.org # ###################################################### Here is the bug report I got from Professor Hsuan, and my exploration of it. The cause of the

Dear list, I may have found a bug in model.frame.default (called by the lm function). The problem arises in my R dev version but not in my R 2.4.0. Here is my config : > version _ platform x86_64-unknown-linux-gnu arch x86_64

I hope you can help me. I am new to "R" and am working through the example of quantile regression analysis in the paper by Roger Koenker (quantile regression in R: a Vignette. I got to page 6 and in trying to execute the following code I received the error below. I have searched on the web and while I have found similar error messages, I have not been able to ascertain the root cause.

I'm trying to learn to use manova(), and don't understand why none of the following work: > data(iris) > fit <- manova(~ Species, data=iris) Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : incompatible dimensions > fit <- manova(iris[,1:4] ~ Species, data=iris) Error in model.frame(formula, rownames, variables, varnames, extras,

Dear list members, This is probably a naive question (from a new user), but I'm having problems getting column names to be recognised in loop that is supposed to apply a set of calculations in turn to each column of a dataframe. The dataframe has 33 columns of output from a simulation. The columns each have a text name and I want to apply the same set of calculations to each column in turn.

Hi, I'm using a loop to evaluate several models by taking adjacent variables from my dataframe. When i try to get predictions for new values, i get an error message about a missing variable in my new dataframe. Below is an example adapted from ?gam in mgcv package library(mgcv) set.seed(0) n<-400 sig<-2 x0 <- runif(n, 0, 1) x1 <- runif(n, 0, 1) x2 <- runif(n, 0, 1) x3 <-

hi all i have a function that uses two inputs, say xdata and ydata. An example is the following, simple1<-function(xdata,ydata) { ofit<-lm(ydata~xdata) list(ofit) } say i use arbitray number for xdata and ydata such that D = x1 x2 y 1 1 10 2 6 6 3 10 7 x<-D[,1:2] and y<-D[,3] if one uses these inputs and rund the program we get the following: >simple(xdata=x,ydata=y)

Hello friends and fellow R users, I have a problem to which I have been unable to find a solution: I am gathering survival data on patients undergoing treatment with a new kind of stent. I want to generate survival data and plot survival curves of these patients based (among other things) on the treating physician. My data set has been tabulated in the following manner: Date (the date the stent

Full_Name: Ben Cooper Version: 1.5.0 OS: linux Submission from: (NULL) (134.174.187.90) The help page for glm.nb (in MASS package) says that it takes "Any other arguments for the glm() function except family" One such argument is start "starting values for the parameters in the linear predictor." However, when called with starting values glm.nb returns: Error in

Hi all, I figured out some things about an earlier question I posted today; this is a new question: I have a for-loop in which I have the following lines for (i in list.files(pattern = "*.dmp")) { currentfile <- read.table(i) vecrawrts <- as.vector(currentfile$V3) ... } Now, vecrawrts is numeric, but if I run the script on all the files with extension .dmp, I get the

I'm a bit confused about the nls package, I'm trying to use it for curve fitting. First off, the documentation for nls says ``see `nlsControl' for the names of the settable control values'' -- this is wrong, it should be nls.control (minor point but had me confused for a moment). Now I'll try something very simple (maybe too simple):

I tried this : > mosaicplot(stoc ~ q9r + segca,data=tmp2,color=T) : works fine. And now, this : > mosaicplot(stoc ~ q9r + segca, data=tmp2, color=T, main="Big title") Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid variable type I'm probably stupid and missed something simple in the manual (and wouldn't like to be

Dear helpers I'm having a problem with function HoltWinters from package ts. I have a time series that I want to fit an Holt-Winters model and make predictions for the next values. I've already built an object of class ts to serve as input to HoltWinters. But then I get an error; I've used HoltWinters a lot of times and this never hapened > data.HW<-HoltWinters(data.ts) Error

Greetings, I would like to use a cubic spline or smoother to model the fixed effects within nlme. So far the only smoother I have been able to get to run successfully in nlme is smooth(). I tried smooth.spline: fixed=list(lKa~1,lCL~smooth.spline(BSA, df=3)) the error I got was the following. Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid

Dear all, I have tried to merge two cross sections (crossA and crossB) as follows, both crossA and crossB include the same variables, e.g., ID,y,x1,x2: Panel <- merge(crossA,crossB,by="ID") I have, thereafter, tried to run the following pooled OLS: OLS <- lm(Panel$y~Panel$x1+Panel$x2) and have received the following error message "Error in model.frame(formula, rownames,