search for: termplot

Displaying 20 results from an estimated 70 matches for "termplot".

2004 Mar 18
1
termplot has problems with a single term (fix included) (PR#6679)
On Thu, 18 Mar 2004 k.hansen@biostat.ku.dk wrote: > The bug exists on R-1.9.0-alpha compiled the 10/3. > > Termplot has a problem if either the model only contains a single term > or if asked to plot a single term. In addition there are problems with > the option se = TRUE. I can't reproduce this in either R-devel or 1.8.1, and termplot hasn't changed since January. I do example(termplot) termplo...
2018 Jan 10
1
termplot intervals - SE or CI?
Thanks for your prompt reply Duncan. I had indeed assumed they were what the help file says until observation raised doubts, which is why I queried it. >From reading the code for termplot(), it seems that either the predict() function doesn't return the 1x standard error, or the curves plotted by the termplot() function are not 1x standard errors. If they're not 1x standard errors, it seems misleading to call them (e.g. in the help file) "standard errors". The &q...
2010 Aug 26
3
Using termplot() with transformations of x
Hi all I was playing with termplot(), and came across what appears to be an inconsistency. It would appreciate if someone could enlighten me: > # First, generate some data: > y <- rnorm(100) > x <- runif(length(y),1,2) > # Now find the log of x: > logx <- log(x) > > # Now fit two models that are exac...
2000 Oct 27
1
par(ask=T) in termplot() (PR#711)
At 11:57 25/10/00 +0100, Brian Ripley wrote: >> Date: mer., 25 oct. 2000 12:38:55 +0200 >> From: Emmanuel Paradis <paradis@isem.univ-montp2.fr> > >> I think it would be nice to have par(ask=T) set by default in termplot(), >> like it is in plot.lm(). > >Well, it isn't really the default in plot.lm, the default for `ask' being > >interactive() && one.fig && length(which) > 1 && .Device != "postscript" > >Can you supply us with an appropriate patch...
2018 Jan 11
0
termplot intervals - SE or CI?
From ?termplot: col.se, lty.se, lwd.se: color, line type and line width for the ?twice-standard-error curve? when ?se = TRUE?. ...which is findable, but might usefully also be made explicit in the definition of the se= argument. -pd > On 10 Jan 2018, at 23:27 , Eric Goodwin <Eric.Goodwin at ca...
2004 Mar 18
1
termplot has problems with a single term (fix included) (PR#6675)
The bug exists on R-1.9.0-alpha compiled the 10/3. Termplot has a problem if either the model only contains a single term or if asked to plot a single term. In addition there are problems with the option se = TRUE. Analysis: termplot starts with terms <- if (is.null(terms)) predict(model, type = "terms", se = se) else predict(model,...
2006 Jul 15
1
termplot and ylim
Hi together, I always get an error message with using ylim in termplot(), like this: > x<-(1:10) > y<-(10:1) > l<-lm(y~x) > termplot(l,ylim=c(1,2)) Is this a bug, or is there another possibility to do that? Especially, I would like to use term.plot() for gamlss objects. Thanks for your help! Andreas -- Echte DSL-Flatrate dauer...
2000 Oct 25
2
par(ask=T) in termplot()
Hi, I think it would be nice to have par(ask=T) set by default in termplot(), like it is in plot.lm(). Best, Emmanuel Paradis -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body&...
2004 Jan 22
1
Bug in termplot?
Hi, Is this a bug in termplot, or (once again) do I just not understand what R is really doing? I am using termplot to contruct partial residual plots, 1. For all terms at once 2. One term at a time but I get different results from these two methods. To give a concrete example, I would have thought the top and bottom row...
2016 Jun 28
2
termplot intervals - SE or CI?
Hello, A reviewer queried what the intervals were on the termplot I provided in a report. The help file for termplot() suggests they're standard errors (se=T), but in the code the se.fit values from predict() are multiplied by 2, suggesting it's a rough 95% confidence interval, is that right? Many thanks, Eric Goodwin Scientific data analyst | Coastal...
2012 Oct 19
1
Addition of plot=F argument to termplot
I have a suggested addition to termplot. We have a local mod that is used whenever none of the termplot options is quite right. It is used here almost daily for Cox models in order to put the y axis on a risk scale: ---- fit <- coxph(Surv(time, status) ~ ph.ecog + pspline(age), data=lung) zz <- termplot(fit, se=TRUE, plot=FALSE)...
2011 Dec 14
1
termplot & predict.lm. some details about calculating predicted values with "other variables set at the mean"
I'm making some functions to illustrate regressions and I have been staring at termplot and predict.lm and residuals.lm to see how this is done. I've wondered who wrote predict.lm originally, because I think it is very clever. I got interested because termplot doesn't work with interactive models: > m1 <- lm(y ~ x1*x2) > termplot(m1) Error in `[.data.frame`(mf, , i)...
2009 Jan 19
1
termplot
I have used glm and stepAIC to choose a best model. I can use termplot to assess the contribution of each explanatory variable in the glm. However the final model after running stepAIC includes interaction terms, and when I do termplot I get "Error in `[.data.frame`(mf, , i) : undefined columns selected". I also see the termplot detail saying "Nothing...
2003 Dec 16
0
Help w/ termplot & predict.coxph/ns
...continuous variable. I am using natural splines (ns) to model the effect of x2. I would like to examine the estimated effect of x2 on the hazard. I have tried various approaches (below; let model.fit= fitted model using coxph in survival library): 1. The simplest method appears to be using termplot(model.fit). This appears to work fine as long as I include the treatment variable in the model. However, without the treatment effect in the model termplot returns a "?" prompt without plotting the effect of x2. This happens whenever I consider any model containing only one predictor v...
2007 Jul 02
2
termplot with uniform y-limits
Does anyone have, or has anyone ever considered making, a version of 'termplot' that allows the user to specify that all plots should have the same y-limits? This seems a natural thing to ask for, as the plots share a y-scale. If you don't have the same y-axes you can easily misread the comparative contributions of the different components. Notes: the current vers...
2007 Jul 02
2
termplot with uniform y-limits
Does anyone have, or has anyone ever considered making, a version of 'termplot' that allows the user to specify that all plots should have the same y-limits? This seems a natural thing to ask for, as the plots share a y-scale. If you don't have the same y-axes you can easily misread the comparative contributions of the different components. Notes: the current vers...
2011 Jun 06
1
Possible bug in termplot function (stats package) ?
Hi all, I noticed some very odd behaviour in the termplot function of the stats package due to the following lines : 18. if (is.null(data)) 19. data <- eval(model$call$data, envir) This one will look in the global environment, and renders the two lines after this 20. if (is.null(data)) 21. data <- mf completely obsolete. If not...
2001 Jan 27
1
termplot fails for composite non-factor terms (PR#828)
I am running R 1.2.1 under Windows 98SE. termplot() currently fails when there are composite terms, thus: > library(mass) > data(hills) > hills.lm <- lm(time ~ climb + poly(dist, 2), data = hills) > termplot(hills.lm) Hit <Return> to see next plot: Error in eval(expr, envir, enclos) : Object "dist" not f...
2004 Jan 04
0
termplot; failure to subset non-dataframe carriers (PR#6327)
termplot() does not carry subsetting over to carriers that are in the environment but not in the data frame. This generates a "subscript out of bounds" error. > data(ToothGrowth) > logdose <- log(ToothGrowth$dose) > tooth.lm <- lm(len ~ logdose, data=ToothGrowth) > termplot...
2007 Jul 03
1
termplot - changes in defaults
While termplot is under discussion, here's another proposal. I'd like to change the default for partial.resid to TRUE, and for smooth to panel.smooth. I'd be surprised if those changes were to break existing code. John Maindonald email: john.maindonald at anu.edu.au phone : +61 2 (6125)3...