similar to: ERROR: no applicable method for "TukeyHSD"

How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Cue Direction (2) x Brain Hemisphere(2) Each of these has 2 levels, 'left' and 'right', so it's a simple 2x2 design matrix. We have 8 subjects in each cell (a balanced design) and we want to specify the interaction contrast so that: CueLeft>CueRght

I'm an experimental psychologist and when I run ANOVA analysis in SPSS, I normally ask for a test of non-sphericity (Box's M-test). I also ask for output of the corrections for non-sphericity, such as Greenhouse-Geisser and Huhn-Feldt. These tests and correction factors are commonly used in the journals for experimental and other psychology reports. I have been switching from SPSS to R

Hi all, I have a strange problem and rigth now I can't figure out a solution. Trying to calculate an ANOVA with one between subject factor (group) and one within (hemisphere). My dependent variable is source localization (data). My N = 25. My data.frame looks like this: > ML.dist.stack subj group hemisphere data 1 1 tin left 0.7460840 2 2 tin left

Dear Experts, I am new to R and unfortunately cannot start with a simply statistical analysis: I manually determined the volume of the right and left hippocampus in a group of meditators and in a group of controls. My data-sheet looks as follows: observation subject group age gender hemisphere volume 1 am04 m 25 f left 3.637 2 am04 m 25 f right 3.713 3 ao08 m 47 m left 3.715 4 ao08 m 47

Hi, we need some help to define planned comparisons. I've based my understanding of the problem on reading Tabachnick and Fidell (2006), ie: http://www.ablongman.com/catalog/academic/product/0,1144,0205459382,00.html I don't understand how to specify planned comparisons in R. I've not found explanations for this in MASS or elsewhere. There is only discussion of the contrast

Dear List, I have a richness data distributing across 20 N to 20 S latitude. (120 E-140 E longitude). I would like to draw the richness in the north hemisphere and a regression line in the plot (x-axis: latitude, y-axis: richness in the north hemisphere). The above demand is done using plot. Then, south hemisphere richness and regression are required to be generated using the same y-axis above

I'm trying to read in a CSV, with lines looking like: HEADER, Latitude DecDeg, Latitude Hemisphere, Longitude DecDeg, Longitude Hemisphere, Speed knots, Bearing Degrees, fixQualityGga, noOfSatGga, altGga, heightGga, selectionGsa, fixGsa, pdopGsa, hdopGsa, vdopGsa, noOfSatGsv, Time, *checkSum $GPS, 3747.0224, N, 12223.4522, W, 0.36, 348.21, 1, 3, 4.01, 175.5, -25.2, A, 2, 4.14,

Hi I hope this is going to the right place. I am trying to write a program which uses KernSmooth library to estimate mutual information between two time series at various different lags. At the moment it’s producing negative values, which is supposed to be impossible (something is fishy). I am summing across one row of the matrix to get p(value is in bin x) and summing across the columns to get

With respect to calculating the epsilon index of sphericity for ANOVA, discussed on pp. 45-47 of: http://www.psych.upenn.edu/~baron/rpsych.pdf It notes that epsilon is not required for a repeated measures design with only k=2 levels, as the minimum value of epsilon (e) is given by: e = 1/(k-1) so for k=2, we have e = 1 (ie, no correction of the F test df; see p. 46). These notes apply to a

Dear R-developeRs, Attached follows a patch against svn 34959 that adds the printing of p-values to the TukeyHSD.aov function in stats package. I also updated the corresponding documentation file and added a 'see also' reference to the simint function of the multcomp package. As it was already brought up in a previous thread [1] in R-help, one can obtain the adjusted

> On 18 May 2017, at 14:58 , Martyn Plummer <plummerM at iarc.fr> wrote: > > > >> On 18 May 2017, at 14:51, peter dalgaard <pdalgd at gmail.com> wrote: >> >> >>> On 18 May 2017, at 13:47 , Joris Meys <jorismeys at gmail.com> wrote: >>> >>> Correction: Also dlt uses the default timezone, but POSIXlt is not recalculated

Dear useRs, I'm working on multiple comparison design on two factor (2 3 levels) ANOVA. Each of the tests I have tried (Tukey, multcomp package) seem to do only with one factor at a time. fm1 <- aov(breaks ~ wool * tension, data = warpbreaks) tHSD <- TukeyHSD(fm1, "tension", ordered = FALSE) $tension diff lwr upr p adj M-L -10.000000 -19.35342

Hey, I have this TukeyHSD output from which I would like to extract only the P-values (p adj, last number). The problem is that the test output is a character list. How can I "break" this sentence to separate the Pv? Tukey multiple comparisons of means 95% family-wise confidence level Fit: aov(formula = Fe1$Fe ~ Fe1$genotype) $`Fe1$genotype` diff lwr upr

In the autumn and winter it's considered bad form to set power state to D3cold, it might cause the device to freeze to death. This is also the case in the other seasons, so any device in the southern hemisphere is affected too. D3cold is not a valid state in a call to pci_set_power_state, only up to D3hot is allowed. Signed-off-by: Maarten Lankhorst <maarten.lankhorst at canonical.com>

The third report on the R-spatial evolution project has been published. https://r-spatial.org/r/2023/04/10/evolution3.html It links back to earlier blogs and presentations, and focuses on work that maintainers of R packages and workflows using rgdal, rgeos and maptools need to put in train now before the packages are archived on CRAN. The first changes will occur in June 2023, and the

The third report on the R-spatial evolution project has been published. https://r-spatial.org/r/2023/04/10/evolution3.html It links back to earlier blogs and presentations, and focuses on work that maintainers of R packages and workflows using rgdal, rgeos and maptools need to put in train now before the packages are archived on CRAN. The first changes will occur in June 2023, and the

Dear R users, I am currently working in subsetting a zooreg() object using either window or subset. I have a solution but it may be a bit cumbersome when I start working with actual data. Your inputs would be greatly appreciated. Example: I have a zooreg() object that starts in 1997 and ends in 2001. This object contains daily data for the 4 years

Happy New Year in advance by a few ticks for the northern hemisphere. Here's the first topic and guest for 2009: In any voice path there are several potential sources of quality problems, ranging from echo to voice dropouts and everything in between. With VoIP systems the potential for quality problems increases dramatically, often times making it very difficult to identify the source of

I am having trouble in Sweave with input line lengths. For example, I may have in my input file the chunk <<>>= BrainSections <- levels(AggData$sctn)[grep( "(^BRAIN)|(^WHOLEBRAIN)|(LEFT HEMISPHERE)| (HALFBRAIN)", levels(AggData$sctn))] @ This is translated in the tex file: \begin{Sinput} > BrainSections <-

> > On 18 May 2017, at 14:58 , Martyn Plummer <plummerM at iarc.fr> wrote: > > > >> On 18 May 2017, at 14:51, peter dalgaard <pdalgd at gmail.com> wrote: > >>> On 18 May 2017, at 13:47 , Joris Meys <jorismeys at gmail.com> wrote: > >>> > >>> Correction: Also dlt uses the default timezone, but POSIXlt is not recalculated