similar to: expand.model.frame fails when call creating model has no data (PR#1423)

Dear R-help list members, I'm encountering problems using expand.model.frame. To illustrate, consider the function > test <- function(model){ + expand.model.frame(model, "income") + } > The data frame Prestige (from the car library) has several variables, including prestige, income and education. I've attached this data frame >

Dear R list members, I'm encountering a problem with expand.model.frame(): Suppose that I define the following simple function (meant just to illustrate the problem): > fun <- function(model){ + expand.model.frame(model, all.vars(formula(model))) + } > and I have the following model, created with an explicit data argument: > mod Call:

Full_Name: Gregory R. Warnes Version: 1.2.0, 1.2.3 OS: SunOS gsun124 5.8 Generic_108528-03 sun4u sparc SUNW,Ultra-5_10 Submission from: (NULL) (12.18.36.49) When using expand.model.frame on a model that specifies a subset selection, an error is generated on the variable used for the subset selection. Example: > data <- data.frame(x=1:10,y=1:10,z=1:10,m=1:10) > model <- lm( y ~

> Thanks. This is also present in the current pre-1.3.0. Your patch > looks correct, but I wonder if the default for "enclos" > should not be > environment(formula(model)) rather than parent.frame() as it is now? > > (And wouldn't it be better named "envir"?) Peter, I was merely making an incremental improvement, your suggestions provide

Dear R-help list members, Yesterday, I described a problem that I encountered with the use of expand.model.frame when the model object passed to this function was created without an explicit data argument in the call. I've had no responses, but I've been able to work around the problem by making a change to expand.model.frame, replacing the line rval <-

Hello All, The version of R I am using is as follows > version _ platform x86_64-pc-linux-gnu arch x86_64 os linux-gnu system x86_64, linux-gnu status major 2 minor 14.1 year 2011 month 12 day 22 svn rev 57956 language R version.string R version 2.14.1 (2011-12-22) I just few days

Dear R-helpers, I'm running Sweave() on a file. First run: > Sweave('20080331.Rnw') Writing to file 20080331.tex Processing code chunks ... 1 : term hide (label=setup) 2 : echo term verbatim (label=oatvar) 3 : echo term verbatim (label=oatvar1) 4 : echo term verbatim (label=oat2wt) 5 : echo term verbatim (label=oat2wt) 6 : echo term verbatim (label=lm) 7 : echo term

I write a function to get the frequency and prop of a variable. freq<-function(x,digits=3) {naa<-is.na(x) nas<-sum(naa) if (any(naa)) x<-x[!naa] n<-length(x) ta<-table(x) prop<-prop.table(ta)*100 res<-rbind(ta,prop) rownames(res)<-c("Freq","Prop") cat("Missing value(s) are",nas,".\n") cat("Valid case(s)

Dear R-devel list members, I'd like to suggest a more flexible procedure for generating contrast names. I apologise for a relatively long message -- I want my proposal to be clear. I've never liked the current approach. For example, the names generated by contr.treatment paste factor to level names with no separation between the two; contr.sum simply numbers contrasts (I recall an

Dear Trevor, I've noticed a discrepancy in the degrees of freedom reported by gam() from the gam package in R vs. gam() in S-PLUS. The nonparametric df differ by 1; otherwise (except for things that depend upon the df), the output is the same: --------- snip ------------ *** From R (gam version 0.93): > mod.gam <- gam(prestige ~ lo(income, span=.6), data=Prestige) >

Dear r-devel list members, I've recently encountered the following problem using predict() with a model that has raw-polynomial terms. (Actually, I encountered the problem using model.frame(), but the source of the error is the same.) The problem is technical and concerns the design of poly(), which is why I'm sending this message to r-devel rather than r-help. To illustrate:

Hello, Much clearer now, thanks. It's a matter of changing the function boot calls to return the predicted values at the point of interess, education = 50, income = 75. I have changed the way the function uses the indices a bit, the result is the same, it's just the way I usually do it. pred.duncan.function <- function(data, indices) { mod <- lm(prestige ~ education +

I am getting an error message from scatterplot: > library(car) > scatterplot(Prestige$income~Prestige$type) Error in Summary.factor(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, : range not meaningful for factors In addition: Warning message: In Ops.factor(x[floor(d)], x[ceiling(d)]) : + not meaningful for factors > The command does output the kind of graph that I want (boxplots).

Hello Rui, Thanks for your helpful suggestions. Just for illustration, let's use the well known Duncan dataset of prestige vs education + income that is contained in the "car" package. Suppose I wish to use boot function to bootstrap a linear regression of prestige ~ education + income and use the following script: duncan.function <- function(data, indices) {data =

Hi In the following, the graph I see on the screen and the .png output coincide. However, in the .pdf file, the fonts seem to be scaled fairly larger, resulting in the label for the top legend disappearing. Is this an infelicity or bug, or is there something I've missed? More generally, how do I control the size of fonts used in legends and axis labels? library(car) library(ggplot2) qp

I noticed that when I fit a linear model using 'lm' there is an attribute called "factors" that is added to the "term". It doesn't seem to appear for 'model.matrix', just 'lm'. I have been unable to find where it gets constructed or what it means? It looks like a two dimensional array that I may be able to use so I would just like to get some

Dear List, I don't get scatter3d to color the sheres according to the '|' argument. library(car) scatter3d(prestige ~ income + education|type, data=Prestige) The spheres on my screen are all colored the same and they are not conditional on Prestige$type. On the other hand: Fit3d and Ellipse3d are colored according to the group argument. rgl_0.92.879 car_2.0-12 R version 2.15.0

Dear all, I'm trying to use "scatterplot" from package "car" to get a scatterplot with row-names as labels and I want my data to be divided into two groups defined by the column "tresh". Leaving away the groups-command gives me the right labels but they are not distinguished by colors according to "tresh". Does anybody know why the following example

Dear list users How is it possible to visualise both a linear trend line and a quadratic trend line on a plot of two variables? Here my almost working exsample. data(Duncan) attach(Duncan) plot(prestige ~ income) abline(lm(prestige ~ income), col=2, lwd=2) Now I would like to add yet another trend line, but this time a quadratic one. So I have two trend lines. One linear trend line

Thanks to Marc Schultz I found the documentation on the "factors" attribute under ?term.object. It stats: factors: A matrix of variables by terms showing which variables appear in which terms. The entries are 0 if the variable does not occur in the term, 1 if it does occur and should be coded by contrasts, and 2 if it occurs and should be coded via dummy