similar to: Lattice: how to index in a custom panel function?

I wonder if someone could explain the behavior of the anova() and lm() functions in the following situation: I have a standard 3x2 factorial design, factorA has 3 levels, factorB has 2 levels, they are fully crossed. I have a dependent variable DV. Of course I can do the following to get the usual anova table: > anova(lm(DV~factorA+factorB+factorA:factorB)) Analysis of Variance Table

Hi there, I wonder if anyone can help me. I'm trying to use bargraph.CI in the Sciplot package when there is a missing combination of the factor levels. Unfortunately the standard errors on the plot do not appear to be correct. Consider an analysis consisting of two factors A and B. When all factor level combinations are present all appears fine: library(sciplot) #all data

How to generate a new factor variable by two other factor variables? For example, if I have two factor variables, factorA and factorB, factorA factorB 0 0 0 0 1 0 0 1 1 1 Is there a simple way to generate a new 4-levels factor variable as factorC factorA factorB 0 0 0 0 0 0 1 1 0 2 0 1

Hello list, I wish to perform a randomization test on the F-statistics of a 2 way ANOVA but have not been able to find out how to do so - is there a package / function that can perform this that I am unaware of? FactorA has 6 levels (0,1,2,3,4,5) whereas FactorB has 3 (1,2,3). A sample: Resp. FactorA FactorB 2 0 2 3 1 2 1 2 2 0 3 2 0 4 2 0 5 2 4 0 1 6 1 1 1 2 1 0 3 1 1 4 1 0 5 1 2 0 2 3 1 2 1

Could anybody explain this difference in the function order() between R under windows and R under Linux ? It seems that under Linux the order in of character dsitinguishes between caps and normal letters and sorts them starting with the capitals (see below). How can I avoid this (I mean: get the same results as I get under Windows)? Thanks Joris Under Windows R 1.6.2 > mstring <-

Dear All, I'm struggling a little with the behaviour of R with GLM interactions. In particular, I have a dataset with two factors - call them factor A and factor B, where I would like to fit a GLM that is factor A + (grouped factor A):factor B. To try to isolate this, I've ignored the original "factor A" part, as that I have this as a separate column in my data. So, it

Can anybody explain me what is going wrong with the rbind command or how I could avoid the problem? I have a script with loops which time has a dataframe as output (called T2sbdata). After each loop I would like to append that dataframe to the results I obtained from the previous loops (called T2sbdataALL) and use for that purpose rbind. Below a part of the loop and the result I get. The

Hello every body, I am trying to do a factorial anova analysis following this model: model<-anova(lm(responsevariable~factorA*factorB)) model<-anova(lm(luz$dosel~luz$estado*luz$Bosque)) Df Sum Sq Mean Sq F value Pr(>F) estado 1 6931.1 6931.1 41.6455 7.974e-06 *** Bosque 1 36.6 36.6 0.2197 0.6456 estado:Bosque 1 36.6 36.6 0.2197 0.6456 Residuals

Dear all, I try to run the example of lmer and get the following error message. > library(lme4) > example(lmer) lmer> (fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)) [[1]] Error in get(x, envir, mode, inherits) : variable "as.dpoMatrix" was not found This error message is similar to what I get with other models. It looks like the mer class has a slightly

Dear list, i have a quick and (hopefully) straightforward question regarding the plot-function after running aov. if i plot an equation like this: plot(dataSubjects~factorA, data=mydata) R gives me the boxplots for this particular factor A. my model, however contains several factors. is there a straightforward way to plot barplots for a specific factor with the constraint that those values

Dear expeRts, I am trying to wrap up a package "utilities" (for my internal use). After adding a function datNAtreat that uses mapply, R CMD check gives WARNINGs for "S3 generic/method consistency", "checking replacement functions" and?"checking foreign function calls", all of which are accompanied by the following error message: Error in .try_quietly

Maybe I do not explain well what I would like to do. I do not want to change the labels of the axis, but the scale. What I want is a general procedure for changing the scale. Its like using a logarithmic scale on a plot. Labels are the same, but the increases of x along the x-axis are defined by a known monotone and continuous function. Florent Bresson ----- Message d'origine ---- De :

Based on the examples I've seen in using statistical analysis packages such as lmer, it seems that people usually tabulate all the input data into one file with the first line indicating the variable names (or labels), and then read the file inside R. However, in my case I can't do that because of the huge amount of imaging data. Suppose I have a one-way within-subject ANCOVA with

Hello everyone, I would like to perform a 2-way manova test, but I'm having some issues. I implemented like this Y<-cbind(Resp1,Resp2,Resp3,....,Respn) model<-manova(Y "tilda" FactorA*FactorB) summary.aov(model) 1. I don't know at what level I have to do the Type I error correction. Is it on p-values returned by "summary.aov(model)? Or is it when I compare each

Suppose you have 4 related vectors: a.id<-c(1:25, 1:25, 1:25) a.vals <- c(101:175) # same length as a.id (the values for those IDs) a.id.levels <- c(1:25) a.id.ratings <- rep(letters[1:5], times=5) # same length as a.id.levels What I would like to do is specify a rating from a.ratings (e.g. "e"), get the vector of corresponding IDs from a.id.levels (via

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100214/4b41a017/attachment.pl>

How would one server out to excel a html table with the mime/content type application/vnd.ms-excel ? In PHP it would be something like ?php header("Content-Type: application/vnd.ms-excel"); ------ But how do I do this in rails? Cheers Glenn -- Posted via http://www.ruby-forum.com/.

If a user enters the Textile-based code for an image, say: !>graphic.png! the textile() method will successfully render a right-float image. Now, if I want to escape the user''s text using h(), the ''>'' is turned into an &amp;. Obviously, the image then fails to render in a browser. Is there a work around for this? I guess what I''m really asking

Hello, >titletool<-read.csv("TotalCSVData.csv",header=FALSE,sep=",") > class(titletool) [1] "data.frame" >titletool[1,1] [1] Experiment name : CONTROL DB AD_1 >t<-titletootl[1,1] >t [1] Experiment name : CONTROL DB AD_1 >class(t) [1] "character" now i want to create an object (vector) with the name "Experiment name :

Dear all, I'm stuck with the following problem: I generate graphs using both the grid system (with lattice) and the base system. I'd like to be able to identify these graphs later on with a bit of identifying text (e.g. a date and some comments). Adding text to these graphs cannot be done using a common system if you want to save them as emf files. I now use: mtext("labelling