similar to: Doubts on anova and use of contrasts in multcomp package

Dears members of R list, It would like that a more experienced statician in R helped me to complete the analysis to follow: r = gl(3, 8, label = c('r1', 'r2', 'r3')) e = rep(gl(2, 4, label = c('e1', 'e2')), 3) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4, 19.6, 21.1, 19.0, 18.6, 22.8, 19.4, 18.8, 19.2, 19.8, 21.4,

Dears members of R list, It would like that a more experienced statician in R helped me to complete the analysis to follow: r = gl(3, 8, label = c('r1', 'r2', 'r3')) e = rep(gl(2, 4, label = c('e1', 'e2')), 3) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4, 19.6, 21.1, 19.0, 18.6, 22.8, 19.4, 18.8, 19.2, 19.8, 21.4,

Dears members of R list, I would like that a more experienced R user help me to complete this analysis: r = gl(3, 8, label = c('r1', 'r2', 'r3')) e = rep(gl(2, 4, label = c('e1', 'e2')), 3) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4, 19.6, 21.1, 19.0, 18.6, 22.8, 19.4, 18.8, 19.2, 19.8, 21.4, 22.8, 21.3) df =

I have a set of data with 12 readings for temperature per day, with 180 days. I want to find the average temperature of each day. I am able to do this one by one, but with that many days to calculate the average for, it will get very long. I'm sure there is a faster way to do this, I just don't know how. What i have so far is: av1 <- subset(ER9r, Day == 98, select = c

Dear R list, Please, could someone help me with SSQ decomposition and contrasts. Below my data, graphic, ANOVAs and my doubt: # Data a = paste('a', gl(3, 8), sep='') b = paste('b', gl(2, 4, 24), sep='') tra = sort(paste('t', rep(1:6, 4), sep='')) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4, 19.6,

# Dear R list, # # A have a doubt about SSQ decomposition and contrasts with ANOVAs. # So, I would like a tip from more advanced R users. # Below my data, the basic script and my doubts: # Data r = paste('r', gl(3, 8), sep='') e = paste('e', rep(gl(2, 4), 3), sep='') tra = sort(paste('t', rep(1:6, 4), sep='')) y = c(26.2, 26.0, 25.0, 25.4,

Problem is demonstrated with this code, intended to find the approximate 'x' at which the 'y' is midway between the left and right asymptotes. This particular data set returns NA, which is a bit silly! -------------- sXY <- structure(list(x = c(0, 24, 27, 48, 51, 72, 75, 96, 99), y = c(4.98227, 6.38021, 6.90309, 7.77815, 7.64345, 7.23045, 7.27875, 7.11394, 6.95424)), .Names =

Hi Eric, you can use an annotate-layer, eg ind<-which(sig>0,arr.ind = T) ggplot(m1.melted, aes(x = Month, y = Site, fill = Concentration), autoscale = FALSE, zmin = -1 * zmax1, zmax = zmax1) + geom_tile() + coord_equal() + scale_fill_gradient2(low = "darkred", mid = "white", high = "darkblue",

Using the small reproducible example below, I'd like to know if one can somehow use the matrix "sig" (defined below) to add a black outline (with lwd=2) to all pixels with a corresponding value of 1 in the matrix 'sig'? So for example, in the ggplot2 plot below, the pixel located at [1,3] would be outlined by a black square since the value at sig[1,3] == 1. This is my first

Dear R users, My immediate problem is that I do not seem to understand how to use the function scale_x_datetime function in ggplot. My deeper problem is that I don?t understand the use of the POSIXct class. ? My data looks like like the following : T????????? val 2009-11-04 23:59:57?? 972.357117???? 2009-11-05 00:02:21?? 969.389221???? 2009-11-05 00:04:45?? 868.817749???? 2009-11-05 00:07:09??

When plotting over multiple pages in lattice, I'd like to be able to have "same" scales within a page, but "free" scales between pages. In other words, something like: z<-data.frame(x=1:100, y=runif(100), d=rep(1:2,50), p=rep(1:2,each=50)) plot<-xyplot(y~x|d*p, data=z,

Full_Name: Quanah Gibson-Mount Version: 2.3.0, 2.3.1 OS: Linux Submission from: (NULL) (171.66.155.86) libc used to export the symbol: __libc_stack_end; however, newer versions of libc no longer export it. R has some serious problems around this, because the code has: # ifdef linux extern void * __libc_stack_end; # endif and #if defined(linux) R_CStackStart = (uintptr_t)

Hi Thierry, On 2 March 2016 at 10:16, Thierry Onkelinx wrote: | Dear all, | | install.packages("lme4") fails with error | Error in read.dcf(file = tmpf) : cannot open the connection | In addition: Warning message: | In read.dcf(file = tmpf) : | cannot open compressed file | '//servername/repository_path/bin/windows/contrib/3.2/PACKAGES', | probable reason 'No such file or

Dear all, install.packages("lme4") fails with error Error in read.dcf(file = tmpf) : cannot open the connection In addition: Warning message: In read.dcf(file = tmpf) : cannot open compressed file '//servername/repository_path/bin/windows/contrib/3.2/PACKAGES', probable reason 'No such file or directory' the repositories set in .Rprofile are

Hi Again, I am writing a 'for loop' to create a matrix of randomly sampled colors. I've written this loop in matlab and it works fine. I then tried to do it in R and apparently there is something wrong with my syntax b/c every time I run the script, the for loop "blows up" at a different point in the loop. As you can see, I ask R to clear my workspace each time, so there

Hello, I have added a new -l16 "extreme" compression quality setting to the NHW Project. This new -l16 quality setting is surprisingly good actually, because I find it better on some images than -l15 quality setting whereas files size is smaller.I notably find this new -l16 quality setting better than x265 (HEVC). More at: http://nhwcodec.blogspot.com/ We can still save 2.5KB per

It's been a few weeks I'm racking my brains on how to use the function glht the package multcomp to test interactions. Unfortunately, the creator of the package forgot to put a sample in pdf package how to do it. I have looked in several places, but found nothing. If someone for the love of God can help me I'll be extremely grateful. The model is glm. -- View this message in context:

Dear List, Could you tell why I get the error message? > library(multcomp) > data("cholesterol") > m = aov(response ~ trt, data = cholesterol) > cht <- glht(m, linfct = mcp(trt = "Tukey")) Error: could not find function "glht" Thank you G?bor

fixed @ R-forge. New version should appear on CRAN soon. Thanks for the report! Torsten > > ---------- Forwarded message ---------- > Date: Sat, 3 Sep 2011 23:56:35 +0200 > From: Ulrich Halekoh <Ulrich.Halekoh at agrsci.dk> > To: "r-help at r-project.org" <r-help at r-project.org> > Subject: [R] glht (multcomp): NA's for confidence intervals using

Hi All, I am trying use 'multcomp' for multiple comparisons after my ANOVA analysis. I have used the following code to do ANOVA: dat <- matrix(rnorm(45), nrow=5, ncol=9) f <- gl(3,3,9, label=c("C", "Tl", "T2")) aof <- function(x) { m <- data.frame(f, x); aov(x ~ f, m) } amod <- apply(dat,1,aof) Now, how can I use