similar to: TRADUCING lmer() syntax into lme()

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

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

Hi, I have a lattice xyplot that contains panels according to FactorA, and curves for the 2 levels of Factor B within a panel. I try to add text in the panels of a lattice graph. I suppose I have to write a custom function (panel.txt). What I really would like is to adapt the text in the panel according to the levels of FactorA. In the manuals, I find examples for the strips using which.given

Hi, As a follow up to Patch D23646 <https://reviews.llvm.org/D23646>, I'm trying to figure out if there should be an alignment check and what the correct approach is. Some background: For stores, the pass turns: %i.vec = shuffle <8 x i32> %v0, <8 x i32> %v1, <0, 4, 8, 1, 5, 9, 2, 6, 10, 3, 7, 11> store <12 x i32> %i.vec, <12 x i32>* %ptr

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

Hi Everybody, please can you help me how to call BayesX in R in order to see the graph already exist in BayesX Thanks ---------- Forwarded message ---------- From: <r-help-request@r-project.org> Date: Sat, Jan 22, 2011 at 5:00 AM Subject: R-help Digest, Vol 95, Issue 22 To: r-help@r-project.org Send R-help mailing list submissions to r-help@r-project.org To subscribe or

Hi Renato, Thank you for the answers! First, let me clarify a couple of things and give some context. The patch it looking at VSTn, rather than VLDn (stores seem to be somewhat harder to get the "right" patterns, the pass is doing a good job for loads already) The examples you gave come mostly from loop vectorization, which, as I understand it, was the reason for adding the

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

Dear all, A question on modelling proportional data in R. I have a test experiment that was designed in a particular way, and which I can analyse "by hand" to an extent. I am really struggling to get R to give me sensible results in modelling it "properly", so must be doing something wrong here. As background, I conduct a series of experiments and count the

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

Dear I hope you are wery well when read this e-mail.My name is Firas ALSHAWY and I am a Phd student in the high Institute of marine research - tishreen university- Syria I would like to want some help my in my research, I have data ( length - weight - gonad weight - fishing mortality - natural mortality - the number of fish-temperature of water, maximum age ) , I would like to use the package (

Hi, I am experiencing a long delay when using dataframes inside loops and was wordering if this is a bug or not. Example code: > st <- rep(1,100000) > ed <- rep(2,100000) > for(i in 1:length(st)) st[i] <- ed[i] # works fine > df <- data.frame(start=st,end=ed) > for(i in 1:dim(df)[1]) df[i,1] <- df[i,2] #takes for ever R: R 2.0.0 (2004-10-04) OS: Linux, Fedora Core 2

Hi everybody, sorry for my english but i'm italian and i don't know it very well. I'm trying to do a java-program to traduce and notify asterisk events to a Tapi program. I've a problem with call trasfer. When i transfer a sip user i would like to put his line on hold but i can't do it. He listen the music on hold but his state is connected and not Hold. I hope you understand

Hello, I'm a member of a French howto translation project (http://www.traduc.org). A member of our team would like to start translating the Samba Howto Collection in French. We're aware of a previous translation project, but it seems to be dead and there is apparently no way to reach its members or see the work it has done: http://home.gna.org/sambadoc/ - Are you aware of a

R-Amigos Mis mejores deseos en estas fiestas y que el año nuevo premie la generosidad con la que los miembros de este grupo comparten su sabiduría y sus trucos. Felicidades -- Sent from my IPAD [[alternative HTML version deleted]]

Zitat von Sébastien Le Ray <sebastien-samba at orniz.org>: >> The very strange thing is, that gpupdate tries to copy somethings >> from \\cch.intra\sysvol and not from \\dc1\sysvol... >> There a no server with name cch.intra, this is just the Realm... > > Thats expected. your.realm should resolve to all your DC in a > round-robin fashion. OK, I didn't

hallo! momentan steht meine seite so, dass ich gewisse dinge mit link_to_remote auf bedarf nachlade. was mach ich jetzt, wenn javascript deaktiviert ist? kann ich das irgend wie mit render partial lösen? kann ich mit rails prüfen, ob javascript aktiviert/deaktiviert ist und dann dahingehend alle sachen sofort laden? hat da jemand ansätze/strategien/lösungen parat? danke! -- Posted via