1,468
questions
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 ...
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=&...
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))
...
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 ...
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 ...
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 ...
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 <-...
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 ...
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=...
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=...
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)
...
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)
...
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 ...
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 ...
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, ...
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 ...
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 ...
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. ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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, ...
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(...
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
...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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. ...
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],...
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 :
...
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 ...
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....
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 ...
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 ...
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 ...
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 ...
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().
...