similar to: effects package x axis labels

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 <-

Dear R-help, I have fitted a glm logistic function to dichotomous forced choices responses varying according to time interval between two stimulus. x values are time separation in miliseconds, and the y values are proportion responses for one of the stimulus. Now I am trying to extrapolate x values for the y value (proportion) at .25, .5, and .75. I have tried several predict parameters, and they

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 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

Gracias a todos por sus respuestas, perdón si no fui muy claro. Lo que intento replicar es un análisis realizado en Stata, tengo que hacer los mismos cálculos pero en el r. En Stata lo que se hizo fue: xi: logistic i.coord i.v11_sexo, vce (cluster red) *Vce (cluster clustvar) especifica que los errores estándar permiten la correlación intragrupo, relajando el requisito habitual de que las

Hi, How do I use a function inside a regression formula? When I tried to use min() function inside a regression formula, I ran into errors. Here are the commands to reproduce the error: mydata <- read.csv(url("http://www.ats.ucla.edu/stat/r/dae/binary.csv")) attach(mydata) names(mydata) mylm2 = lm(admit~gre+I(max(gpa,rank))) *Error in model.frame.default(formula = admit ~

Dear R users, I have some trouble using the allEffects() function to compute and display effect plots for a linear model. My data is quite simple, it concerns effects of 3 treatments on the tumoral volume of mice. vTum codes for the qualitative initial volume, from small to big, temps is the time in month since beginning of treatment, and S?rie codes for the batch. Data is unbalanced. >

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

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

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)

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 <-

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,

I have a rather large data set (about 30 predictor variables) I need to preform a logistic regression on this data. My response variable is binary. My code looks like this: mylogit <- glm(Enrolled ~ A + B + C + ... + EE, data = data, family = binomial(link="logit")) with A,B,C, ... as my predictor variables. Some categorical, some continuous, some binary. I run the code and get

Buenas tardes, desde Uruguay, quiero realizar una regresión logística,con la función GLM, pero ademas quiero designarle a las observaciones clusters. La idea es decirle al programa que las observaciones corresponden a diferentes clusters. Saludos! Sebastián. [[alternative HTML version deleted]]

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

On Thu, 2016-04-14 at 18:07 +0100, Jonathan Hunter wrote: > On 14 April 2016 at 13:37, Jonathan Hunter <jmhunter1 at gmail.com> > wrote: > > > # samba-tool dbcheck --cross-ncs > > Checking 4079 objects > > MYOBJ=value,OU=myou,DC=mydomain,DC=org,DC=uk: 0x00290001 > > MYOBJ=value,OU=myou,DC=mydomain,DC=org,DC=uk: 0x0029000a > >

in my code i''m calling save on an object. The id gets created as the object is persisted. How can i retrieve that newly created id? I guess I could query the table again using the other attributes of the object but is there another way? thanks in advance. larry -------------- next part -------------- An HTML attachment was scrubbed... URL:

Thank you, Andrew - I hadn't done so. (In a good way, I haven't yet had problems with samba that have caused me to delve quite so deeply into the DB :) so I'm not as familiar with the range of tools as I could be, sorry!) This has flagged up quite a few errors, all along the lines of: # samba-tool dbcheck --cross-ncs Checking 4079 objects MYOBJ=value,OU=myou,DC=mydomain,DC=org,DC=uk:

I have a data matrix with genes as columns and samples as rows. I want to create all possible gene ratios.Is there an elegant and fast way to do it in R and write it to a dataframe? Thanks for any help. Som. -- View this message in context: http://r.789695.n4.nabble.com/Calculating-all-possible-ratios-tp4627405.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML

Hi all, I have the following problem: I have a csv-file consisting of timestamp values (no dates), e.g.: Timestamp1;Timestamp2; 05:24:43;05:25:05; 15:47:02;15:47:22; 18:36:05;18:36:24; 15:21:24;15:22:04; I need a vector with the difference of the two timestamps, so I read the data with the read.csv-function: myObj <- read.csv("file.csv",header=TRUE,sep=";"). I have then