Skip to main content
Filter by
Sorted by
Tagged with
0 votes
0 answers
12 views

How can I run a repeated measures anova in R with trial-by-trial data? [migrated]

I need help setting up a dataset and running repeated measures anova in R. My current data set is the trial by trial data for 41 subjects. There is a column for subject ID (1-41), a column for the ...
Sydney Trowbridge's user avatar
0 votes
0 answers
24 views

Contrasts for type III ANOVA on mixed effects model (lme) in R [migrated]

I'm running a series of mixed effects models using nlme in R. Here is an example of my final model: M1.Final <- lme(Response ~Treatment * Species, random = ~1|ID, data = data, method=&...
mels's user avatar
  • 37
0 votes
1 answer
17 views

Error message while running Tukey post hoc test

I ran a One-way Anova using this dataset: data <- data.frame(pH = c(1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5), biosorption = c(30, 40, 34, 50, 51, 54, 60, 61, 62, 70, 71, 72, 80, 81, 82)) ...
Abdullahi Adamu Faggo's user avatar
0 votes
1 answer
34 views

MANOVA analysis in R

I want to do a MANOVA analysis with community data. I think I've understood the purpose of this analysis but not the principle or how to carry it out. These are two random data frames that simulate ...
Klervi's user avatar
  • 11
1 vote
0 answers
51 views

Type III ANOVA for linear mixed models using multiply imputed data

I used R to run mixed linear model on a multiply imputed dataset. I'm trying to evaluate the significance of my main effects. I was not able to find a package/function that would allow to compute a ...
Alexandra Chapdelaine's user avatar
1 vote
0 answers
29 views

anova.mitml.results giving high df2 when comparing nested models [migrated]

I'm comparing nested models using anova.mitml.results to test my fixed effects. I have 25 imputed datasets, 231 participants per time measure x 3 time measures. I'm not familiar with the D3 method of ...
Alexandra Chapdelaine's user avatar
0 votes
1 answer
39 views

How to force R to give the SSR in a single line for all covariates in an ANOVA table?

Background and What I have Tried Suppose I have a model with 2 numeric covariates and 1 categorical covariate with 3 factor levels for a response: set.seed(1) Y <- sample(100) n <- 100 X1 <-...
NM_'s user avatar
  • 2,003
0 votes
0 answers
28 views

Warning: no labels for x. Assuming order of tree$tip.label

I am using the phylANOVA() function in the R package phytools, as follows: phylANOVA(HoloceneTree, HolDat$Category, HolDat$Niche_Breadth, nsim = 1000, posthoc = T) I keep getting the error ...
Birdman's user avatar
  • 45
0 votes
1 answer
35 views

What is SAS doing in this (simple?) script?

This is related to a previous question but here my concern is what SAS is being asked to do by the following script (ignoring i/o commands): proc reg; model hddt=bloc entry/r p; output out=new r=...
Buck Thorn's user avatar
  • 5,064
0 votes
1 answer
33 views

How to compute the predicted response variable after an ANOVA with anova1 in MATLAB?

I want to replicate in MATLAB a simple demo outlined here that uses SAS. The SAS code is ods rtf file='residual plot SAS output.rtf'; proc print; run; proc reg; model hddt=bloc entry/r p; output out=...
Buck Thorn's user avatar
  • 5,064
0 votes
1 answer
28 views

R Package for Displaying Two-Way ANOVA with Mean, SE, and Significance Letters (Both Capital and Small Letters)

I am looking for an R package that can help me display the results of a two-way ANOVA. Specifically, I need to include the following elements in my output: Mean values ( not estimated, originial) ...
sta_107's user avatar
0 votes
1 answer
63 views

Error using mcp2atm function for robust pairwise comparisons - undefined columns selected

I used a robust two-way ANOVA using the WRS2 package (R v.4.4.1), using the following code as suggested by Mair and Wilcox (2020) - the full dataset is on the bottom of this message: library(WRS2) ...
r.filogonio's user avatar
0 votes
1 answer
34 views

How to get sphericity results with Anova() function from Car package

I have an unbalanced dataset (4 groups in total, of which 3 groups consist of 8 rats and one group consists of 7 rats). I want to run a mixed ANOVA (type III) with the Car package. My dataset has 3 ...
Mirte's user avatar
  • 1
0 votes
0 answers
27 views

Suggestions for analyzing synchronization data and running a two-factor ANOVA on circular data on MATLAB

I am running a tapping experiment where I have to analyze whether meter and/or pattern affect the synchronization performance. The within subjects experiment is comprised of 6 conditions, so it's a 2 ...
A. Chavez's user avatar
0 votes
0 answers
29 views

R - How to automatically feed in parameters from aov() to power.anova.test()

I need to implement an assessment of the statistical power of ANOVA tests into my automated analysis workflow in order to point out in which datasets an ANOVA-analysis bears any statistical relevance, ...
anothernoob's user avatar
0 votes
1 answer
18 views

How can I store list of t1waybt ANOVAs (WRS2 package) into a data frame (R)

I have computed 12 ANOVAs using WRS2 package, with the t1waybt function. I would like to be able to put all of these results into a data frame so I can view/print them all. I have tried the following ...
user26264581's user avatar
0 votes
0 answers
15 views

What is the best test to run to compare presence/absence of something between sampling locations?

I have data on the presence/absence of four different pathogens found in edible crabs in 7 locations over 2 seasons. For this study, 30 crabs were collected from each site and for each crab the ...
Kiran B's user avatar
  • 21
0 votes
0 answers
25 views

Problem with factors with a linear model lmer

I have a problem with a linear lmer model. Indeed, I have several factors that interest me and I want to see if they have an impact on my data. I have 6 of them and for 4 of them, I have no worries. ...
Jules Brochon's user avatar
1 vote
1 answer
22 views

r - ANOVA on a distance matrix (nesting required?)

I have an experimental design intended to investigate the influence of two treatment factors (A with 2 levels and B with 8 levels) on the proportions in 8 groups. However, due to technical reasons, I ...
Anti's user avatar
  • 397
1 vote
1 answer
28 views

Discrepancy between grouped and ungrouped output of HSD.test from agricolae package in r

HSD.test with group=TRUE assigned different significance letters to all three groups in the PlantGrowth dataset, but when HSD.test is run again with group=FALSE, the control group shows no significant ...
jtpugs's user avatar
  • 11
0 votes
0 answers
44 views

How to carry out Dunnett's post hoc correction on a 2-way ANOVA in R?

How can I perform Dunnett's post hoc correction after a 2-way ANOVA in R? The Dunnett's test compares every mean to a control mean. Using ToothGrowth dataset as an example, lets just say that that ...
katkat's user avatar
  • 11
0 votes
0 answers
12 views

Appropriate Statistical Test for pre to post data with imputed data sets?

Thank you in advance for your time and help!! I'm working with data with 70 participants that completed an intervention. Participants completed surveys before and after the intervention (so pre and ...
icecreamsundaes's user avatar
0 votes
0 answers
21 views

R - multcompLetters Error undefined columns selected

I have a problem with the multcompLetters4 function. I have done anova and tukey test and now I want to generate letters for a boxplot with the function. I have seen a youtube video and did it exactly ...
user25710345's user avatar
1 vote
1 answer
21 views

Grouped barplots with error bars and pvalue pairwise comparisons

Am new to stackoverflow and would love some help. I have a dataset as below. I would like to create a barplot showing pairwise comparisons between treatment groups as shown but with significance stars ...
Kiplagat Noel's user avatar
0 votes
0 answers
31 views

How to adjust R's aov code for CRD, RCBD and split plot design?

Let's say I have a data file named DATA1 for a 2-factor experiment with A and B columns as factors, and X column is the measured data. So how should I write the aov code for ANOVA of this experiment, ...
MD P's user avatar
  • 1
0 votes
2 answers
83 views

How to make statsmodels' ANOVA result match R's ANOVA result

The following question is sort of a concrete adaptation of a post on StatsExchange. The following R script runs just fine: library(reshape) J1 <- c(9,6,8,7,10,6) J2 <- c(2,1,4,1,5,2) J3 <- c(...
hobscrk777's user avatar
  • 2,547
0 votes
2 answers
46 views

How to extract P-values and append them to an already existing dataframe?

I am trying to write a simple code where it runs one-way ANOVA for each column, where the data looks like this: PROTEIN A PROTEIN B A1 Cell 1 Cell 2 A2 Cell 3 Cell 4 B1 Cell 5 Cell 6 B2 Cell 7 ...
Gabrielle's user avatar
0 votes
1 answer
28 views

Getting Size Error When Putting Dictionary Values In One Way ANOVA Method

When I try to put the values of my dictionary into the SciPy f_oneway function it gives me an error saying that input is too short (must be greater than 1), but the length of it is 3 so it should be ...
Cameron Kranz's user avatar
0 votes
1 answer
28 views

How to access the "heading" in the list of an Anova output

I want to find a way to take the summary output from anova(model) and put it into a table. I haven't been able to find any ways to do this so I am just doing it manually using gt() and data.frames. My ...
computer_goblin's user avatar
0 votes
0 answers
43 views

How to run two-way anova + Tukey post-hoc test over multiple columns in dataframe R

I have a dataset where I have 4 groups (wt, ko (Genotype) and (Treatment) a and b) which I want to do two-way anova for (Gene expression ~ Genotype*Treatment) and subsequently a Tukey post-hoc (if ...
Kim 's user avatar
  • 1
0 votes
0 answers
21 views

Degrees of Freedom in a phylogenetic ANOVA

Does anyone know how to find the degrees of freedom using a phylogenetic anova? Specifically, I am using phylANOVA from phytools. Is it the same way to calculate other degrees of freedom, or does it ...
Tucker Heptinstall's user avatar
0 votes
0 answers
24 views

Issue with Anova, Effects, and Multiple comparisons in JMP

I am trying to run an Anova on a data set I have from our lab. Briefly, we exposed some particles at 3 different sizes to some cells at 2 different doses and measured gene expression.  The issue I am ...
Zak kasuske's user avatar
0 votes
0 answers
28 views

Models with Interaction effect in R - "*" versus ":" [duplicate]

I'm trying to form a model with an interaction effect. My data is related to plant species, the number of flowers the plant produces and the number of pollinators. Using my uni notes my code looks ...
Rí_17's user avatar
  • 1
0 votes
0 answers
4 views

Trying to fit a submodel, What do I do wrong?

lm_full <- lm(CO ~ AT + AP + AH + GTEP + TIT + TAT + TEY, data = Data) lm_modelA <- lm(CO ~ AT + AP + AH, data = Data) returns this error message: i as.data.frame.default(data) : cannot coerce ...
Johnny Minor's user avatar
0 votes
0 answers
60 views

Permanova - adonis2

I have richness (ecology) data for 9 plots. There are three lines, each consisting of three plots. In the first two rows, I have wealth data over 5 years. On the third line, just 3 years. I would like ...
Thielly Schmidt furtado's user avatar
0 votes
0 answers
61 views

Getting an ANOVA table from an lmer applied to data after multiple imputation in R

I've been running an lmer on a data set that has missing values. I've managed to impute the missing values and run my model and successfully and get a summary of the pooled fit. However, I need the ...
Apato Sagashi's user avatar
1 vote
1 answer
80 views

Different Type III Sums of Squares from car Package

I am looking for Type III sums of squares. I understand Type III SS will by default not agree with SAS. However, I have run the following two models and get differing results. I realize my lack of ...
Greg's user avatar
  • 11
0 votes
1 answer
132 views

Looking to filter and/or tidy up TukeyHSD pairwise comparisons for two-way ANOVA

I'm running a 2-way ANOVA to compare how two variables affect the price of a service. However, when I run post-hoc pairwise comparisons I would like to only see rows that are comparable (i.e. when one ...
Phaaltu Waaltu's user avatar
0 votes
0 answers
57 views

Chi squared for three different vectors of unequal sample size

I am looking to perform chi-squared on three different vectors which have an unequal sample size. I have tried to do this as both vectors and combing the data in a dataframe but continue to get error ...
Michael Bennett's user avatar
0 votes
0 answers
32 views

Error in ezANOVA bootstrapping for repeated measures ANOVA

I am having trouble with the following code. I am trying to bootstrap a repeated measures anova. I keep getting an error that my data contains NA values, although no NA values seem to be there - I ...
Researchderp's user avatar
0 votes
0 answers
50 views

How to generate F statistics for permutation for two-way ANOVA using avop in lmPerm package

I am running permutation for two-way ANOVA using aovp in the lmPerm package in R. I can generate the ANOVA table with df, Sum Sq, R Mean Sq, Inter. and P values but not the corresponding F statistics. ...
Casey Mueller's user avatar
0 votes
1 answer
79 views

three-way nested ANOVA

I spent a while trying to solve this problem about analyzing a three-way nested ANOVA. Here is a fake dataset and example with R: set.seed(100) ps=factor(paste("countries_",rep(letters[1:3],...
Eric Wajnberg's user avatar
1 vote
1 answer
285 views

Displaying compact letter display (CLD) from 2-way ANOVA in a boxplot

I used multcompLetters4() to obtain the compact letter display from a two-way ANOVA performed on 4 groups and I am now trying to display my data and the letters as a boxplot. This is the code I used : ...
chiliconcarno's user avatar
0 votes
0 answers
29 views

Why TukeyHSD test keeps returning NA for a linear model in R?

I have a dataframe with four columns: BMI, AGE, SEX and Geno. Geno is a factor variable has 4 levels(A|A,A|G,G|A,ref:G|G). I would like to fit a linear regression model use BMI~AGE + SEX + Geno, and ...
Timon's user avatar
  • 327
0 votes
0 answers
38 views

ANOVA Feature Selection

I got a warning message when trying to use the feature selection. I think something is missing or wrong here, but I can't figure out what it is. Thanks in advance for the help. from sklearn....
Uğur Özalp's user avatar
0 votes
0 answers
18 views

Calling and pasting testing statistics from two independent variables in R

I am looking to run a two-way ANOVA and paste the results for both independent variables in R. I am able to do this with tests using just one independent variable as the results only contain one test ...
Jallen4's user avatar
0 votes
0 answers
28 views

How can I add Tukey TSF p value to plots?

I am trying to do make Tukey HSD plot by adding p values as already answered in one of the previous questions, I am sharing the link to the questions already answered Link. I am trying plot something ...
Mahir 's user avatar
  • 11
0 votes
0 answers
32 views

why stat_compare_means(method = "anova") gives different pval than aov()

everyone I am trying to figure out how does stat_compare_means works I tried to use stat_compare_means(method = "anova") in order to display pval on my boxplot. But then I ran also aov() on ...
Romina Appierdo's user avatar
0 votes
0 answers
26 views

Using anova() on multivariable quantile regression model with natural splines generates error

I created a multivariable quantile regression model: model <- rq(y ~ x1 + x2+ x3... + ns(high_value_x,df=5), tau = seq(0.1, 0.9, by = 0.1), data = df, method='fn') The high_value_x represent a ...
Mikołaj's user avatar
  • 385
0 votes
0 answers
24 views

Differences between anova(), anova.rq(), anova.rqs(), and anova.rqlist()

The quantreg package introduces functions related to the computations of quantile regression models. Among the additional functions, we have the titular anova.rq(), anova.rqs(), and anova.rqlist(). ...
Mikołaj's user avatar
  • 385

1
2 3 4 5
30