similar to: Problems using allEffects() (package effect)

Dear everybody! Need help with graphics. I am runnig a simple lm and then using allEffects from 'effects' package: require(effects) model<-lm(Y~A+B, data=mydataframe) I am trying to build (for each predictor - A and then B) a plot of means on Y. I was successful doing it like this - in one swoop: ml.eff<-allEffects(ml1, se=F) plot(ml.eff,ylab="Title of Y") Is it

Thanks in advance for any help. I am attempting to add points to a plot using the allEffects command in the effects package. When I try to add the points I get the following error message: Error in plot.xy(xy.coords(x, y), type = type, ...) : plot.new has not been called yet Strangely, using the code I've pasted below this has worked for me in the past however figuring out what has

Dear John Fox and everyone, I have been using the effects library with glms and have found it very useful. Now I'm trying it with lms and I'm not sure if the results of the allEffects() are as expected. I've got a model that looks like this: mymodel = lm(formula = A ~ B + C + D + B:D + C:D) Residuals: Min 1Q Median 3Q Max -3.80156 -0.73486 -0.09792

I've got this dataframe X12 X14 X17 X19 vitigni years 1 1 4 78 54 rie 2005 2 1 4 7 4 rie 2005 3 1 4 75 5 rie 2005 4 2 5 66 5 croa 2005 5 1 4 4 46 croa 2005 6 2 5 7 6 croa 2005 7 3 2 56 5 rie 2006 8 3 6 4 7 rie 2006 9 3 2 36 7 rie 2006 10 1 5 3 7 croa 2006 11 1

Dear List, I find Effects package very useful, but I believe I have found a bug in allEffects function. Please consider the following code: test <- data.frame(tries= round(runif(40, 5, 300)), tra = gl(4, 10, labels = c("V", "D", "C", "L")), prop= runif(40, 0, 1)) test$success <- round(with(test, tries*prop)) test$prop <- with(test,

Dear All, probably a simple enough solution but don;t seem to be able to get my head around it...example based on a publicly available data set: mydata <- read.csv("https://stats.idre.ucla.edu/stat/data/binary.csv") mylogit <- glm(admit ~ gre + gpa + rank, data = mydata, family = "binomial") library(effects) plot(allEffects(mylogit) ? ?

Hello, I am not having success in a simple task. Using the effects package, I would like to add reference lines at probability values of 0.1 – 0.6 on a plot of the effects. The plot command works, but following up with an abline command produces the message “plot .new has not been called yet”, and of course the reference lines were not added. Looking through past R help lists, there was a

Dear r heper, How can I change the strip text, for example (16,23] in the following example, to other more informative text such as "high level" on the fly? library(effects) Cowles$ex2 <- cut(Cowles$extraversion,3) mod.cowles <- glm(volunteer ~ sex+neuroticism*ex2,data=Cowles, family=binomial) eff.cowles <- allEffects(mod.cowles) plot(eff.cowles,

This is a topic ' "apparent S3 methods" note in R CMD check ' from R-package-devel https://stat.ethz.ch/pipermail/r-package-devel/2015q2/000126.html which is relevant to here because some of us have been thinking about extending R because of the issue. John Fox, maintainer of the 'effects' package has enquired about the following output from 'R CMD check

I have a process where I am creating a effects plot similar to the cowles effect example. I would like to add the point estimates to the effects plot, can someone show me the correct syntax. I have included the "R" effects example, so you can show me the correct syntax. Thanks mod.cowles <- glm(volunteer ~ sex + neuroticism*extraversion, data=Cowles, family=binomial)

Dear R People: I'm using R2HTML but having a strange result. Here is the original data: resp trt block 90.3 A I 89.2 A II 98.2 A III 93.9 A IV 87.4 A V 97.9 A VI 92.5 B I 89.5 B II 90.6 B III 94.7 B IV 87.0 B V 95.8 B VI 85.5 C I 90.8 C II 89.6 C III 86.2 C IV 88.0 C V 93.4 C VI 82.5 D I 89.5 D II 85.6 D III 87.4 D IV 78.9 D V 90.7 D VI And here are the commands: > resin1.df <-

Hi, I'm running R 2.11.1 on Windows 32-bit. I fitted a glm to a "large" data set with no problems (i.e., nothing is wrong with the created glm object). However, when I try to use the "effects" package to plot the effects I get the error message shown below. I get the same message if I try to plot a single effect. Any idea what could be wrong here? glm.eff <-

>>>>> Duncan Murdoch writes: > On 12/06/2015 4:12 AM, Martin Maechler wrote: >> This is a topic ' "apparent S3 methods" note in R CMD check ' >> from R-package-devel >> https://stat.ethz.ch/pipermail/r-package-devel/2015q2/000126.html >> >> which is relevant to here because some of us have been thinking >> about extending

The R manual says R will not build with gcc on 64-bit Solaris x86 with gcc http://cran.r-project.org/doc/manuals/R-admin.html#Solaris "Tests with gcc32 on ?x86? and ?amd64? have been less successful: ?x86? builds have failed on tests using complex arithmetic33, whereas on ?amd64? the builds have failed to complete in several different ways, most recently with relocation errors for

>>>>> Duncan Murdoch writes: > On 12/06/2015 7:16 AM, Kurt Hornik wrote: >>>>>>> Duncan Murdoch writes: >> >>> On 12/06/2015 4:12 AM, Martin Maechler wrote: >>>> This is a topic ' "apparent S3 methods" note in R CMD check ' >>>> from R-package-devel >>>>

Hello, For basic date time arithmetics operations, AFAK, there're actually the function difftime() and the (dt + num) operations. I'm wondering if other basic operations exist, like add(dt, num, unit) where unit would be "y", "q", "m", etc. Also for the function seq.dates (or seq.POSIXt), the case for by="months" would be more useful if it

Hello, For basic date time arithmetics operations, AFAK, there're actually the function difftime() and the (dt + num) operations. I'm wondering if other basic operations exist, like add(dt, num, unit) where unit would be "y", "q", "m", etc. Also for the function seq.dates (or seq.POSIXt), the case for by="months" would be more useful if it

I currently work on a draft of an aquatic bioassessment. The conditions tested are the following: ER river water T dechlorinated water control 0.5 + 0.5mg / L of malate T + 1 dechlorinated water control + 1g / L of malate T ED dechlorinated water control SED + ER + river water sediment SED ED + sediment + water dechlorinated. It is the result of AChE in muscle (fillet of fish). The production of

Hello Rhelpers, I am trying to represent the following data (hereafter named donnees) in a barchart with a grouping variable : site traitement date s res 1 NT 17/10/2005 normal 76.2 1 T 17/10/2005 normal 103.2 1 NT 23/11/2005 tardif 81.6 1 T 23/11/2005 tardif 98 2 NT 15/10/2005 normal 72.71 2 T 15/10/2005 normal 94.47 2 NT 15/11/2005 tardif 79.65 2 T 15/11/2005 tardif 94.7

Dear Colleagues, John Mandel ( Chain block designs with two-way elimination of heterogeneity. Biometrics 10, 251-272 ,1954). extended the class of chain block designs (Youden & Conner (1953) to elimination of both row and column (blocks) effects. These experimental designs can be useful in engineering and other fields. I am having difficulty obtaining his adjusted treatment means in his