similar to: interaction.plot

Dear R Core Team I've a proposal to improve drop1(). The function should change the contrast from the default ("treatment") to "sum". If you fit a model with an interaction (which ist not signifikant) and you display the main effect with drop1( , scope = .~., test = "F") If you remove the interaction, then everything's okay. There is no way to fit a

Dear R Friends , I want to name my variables( more than 2 variables in a single plot) within a plot to distinct them from each other, but I cann't. How it is possible? I don't mean x and y axis using xlab or ylab. At the below , it follows some lines, only as an example that you could try please, if it is possible. I really thanks for your attention. Amir library(graphics) y<-

Dear R-core There is a different output for the apply function due to the change of unlist as mentioned in the R news. Newly, applying as.factor() (or factor()) in str(dat <- data.frame(x = 1:10, f1 = gl(2,5,labels = c("A", "B")))) (d1 <- apply(dat,2,as.factor)) newly returns a character matrix while in R-2.3.1 the same command resulted in an integer matrix that was

Dear R core team Using the following code produces an empty plot (similar to col = NA): > plot(1:9, col = factor(rep(1:3,3), labels = c("red", "blue", "black"))) My question: Shouldn't one get at least a warning (or an error) if one tries to use a factor as col argument? Thanks for an answer. Regards, Christoph Buser

Dear R-list I have a question concerning the argument "adjustSigma" in the function "lme" of the package "nlme". The help page says: "the residual standard error is multiplied by sqrt(nobs/(nobs - npar)), converting it to a REML-like estimate." Having a look into the code I found: stdFixed <- sqrt(diag(as.matrix(object$varFix))) if (object$method

hi all, I wish to draw on the same graphic device 3 functions. But i don't want them to be on different graph, i want to compare them on the same I don't need mfrow or mfcol, I need something else... 1 graph on 1 device inside this graph 3 ploted function. I saw something unsing data.frame, but i think it's overkill, and something less complicated must exist, if not why? why not

Dear R-users I have a problem choosing initial points for the function arms() in the package HI I intend to implement a Gibbs sampler and one of my conditional distributions is nonstandard and not logconcave. Therefore I'd like to use arms. But there seem to be a strong influence of the initial point y.start. To show the effect I constructed a demonstration example. It is reproducible

Hi all I tried to reproduce an example with lme and used the Orthodont dataset. library(nlme) fm2a.1 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1 | Subject) anova(fm2a.1) > numDF denDF F-value p-value > (Intercept) 1 80 4123.156 <.0001 > age 1 80 114.838 <.0001 > Sex 1 25 9.292 0.0054 or alternatively

Hello! I've tried to use extprograms pipe feature but stuck with "Broken pipe" errors in mail log: May 17 15:18:57 backend1 dovecot: lmtp(20338, user at domain.tld): Debug: 3VMVFWFKllFyTwAArRg8UA: sieve: action pipe: running program: learn May 17 15:18:57 backend1 dovecot: script: Error: write(response) failed: Broken pipe May 17 15:18:57 backend1 dovecot: lmtp(20338, user at

Could somebody recommend some textbook how to compute contrast when there are interactions terms? "Applied Linear Regression Models" (book) mentioned contrast, but I cannot extend it to the case where there are interaction terms.

Hi, I use anova() to find interaction effect and then I need to find data with main effect from those data with interaction effect. How to do that? I used : anova(lm(data~factor1*factor2)), then select data with interaction effect. Then I need to select those data also with main effect of factor1 or factor2, from previous selected data. How to do that? Many thanks for your time on my

Dear R gurus I've a question concerning the transformation of a new S4 class (that I defined) and the existing class "POSIXlt". I did the following: ## Definition of the new class: setClass("dtime", representation(data = "POSIXlt"), prototype(data = as.POSIXlt("2004/06/01"))) ## Transformation between the new class "dtime"

Dear R core team Please apologize for posting the same question twice on R-help and R-devel. Since I was not sure which list is appropriate I tried R-help (Tue Mar 22), but got no answer. Now I do not know if the formulation of my question was unclear or the question is not so easy to answer or to easy (what I do not hope). My problem: I create a new S4 class, containing one slot, data (of

Hy all, I apologize i've used rownames instead of colnames in my first exemple. (that's why i changed the mail object) I was on mars when i wrote my question... Every answer where correct ... but, because i've made a wrong question, people wern't able to understand... So finally i'm speaking about colnames : Let's be more precise: I've got a query that gives me for

Hi, I am trying to plot an interaction.plot with different color for each level of a factor. It has an erratic behavior. For example, it works for the first interaction.plot below, with the example from the ALDA book, but not with the other plots, from the NPK dataset: # from http://www.ats.ucla.edu/stat/r/examples/alda/ch2.htm tolerance <-

Dear R-users, Sorry for reposting. I put it in another way : I want to test random effects in this random effect model : Rendement ~ Pollinisateur (random) + Lignee (random) + Pollinisateur:Lignee (random) Of course : summary(aov(Rendement ~ Pollinisateur * Lignee, data = mca2)) gives wrong tests for random effects. But : summary(aov1 <- aov(Rendement ~ Error(Pollinisateur * Lignee), data =

I would like to annotate my plot with a little box containing the slope, intercept and R^2 of a lm on the data. I would like it to look like... +----------------------------+ | Slope : 3.45 +- 0.34 | | Intercept : -10.43 +- 1.42 | | R^2 : 0.78 | +----------------------------+ However I can't make anything this neat, and I can't find out how to combine this

Hi I just started using RGui.exe under widnows. I have a text file containing date arranged in columns and rows, each column has the same format, each row with different formats. 3 of the columns are something like this 1/12/2006 3:59:45 PM I need to calculate the different in seconds between 2 selected periods using their row’s index My solution: Read the file in a data frame and

Dear R help-list reader, I would like to generate a plot which compares to states in a patient treatment, before and after. for this reason I have generated a vector before<-c(1,30,23,40) and after<-c(20,10,20,60) the first element in "before" corresponds to the first element in "after". I would like and generate a dotplot with before and after as x-scale, the elements

Hello :) I have two data frames, one has properties taken on a piece by piece basis and the other has performance on a lot by lot basis. I wish to combine these two data frames but the problem is that each lot has multiple pieces and hence i need to take a mean of the properties of multiple pieces and match it to the row having data about the lot. I was wondering if there is a simple commmand,