Analysis of Variance (ANOVA) in R

This is a quick tutorial on how to perform ANOVA in R.
I misstated at the end the hypothesis we are testing the means, not variances of the variables. So for this example we reject the null and say there is not enough evidence to suggest different means between groups.
The median line in the boxplot should be close to the mean (if your data is normal which is an assumption in ANOVA)
--------- Data Set -----------
college.cengage...
Thank you very much for your time!
#R #ANOVA

[@powerdriller4124]

I like very much the rustic procedure, no melt(), no ncol(), no fancy jumps. For beginners is good to watch the barefoot process.

[@rabinraut1408]

Did this command work out? plot(weight~team, data=df) I have got error in this command.

[@yousefsonbol2960]

You can use a stack function : df

[@paulandsuziestravels260]

The line on a box plot is the median not the mean and the p-value > 0.05 means that the null hypothesis that the means (not variances) are equal is not rejected.

[@thatrnerd4265]

Thank you for pointing that out. I have added some text in the description.

[@dylanjones2682]

You can just use the stack function from utils to melt the data into a long format in one line of code

[@Nothingimportant1]

Thank you for shooting this video. It was simple and crystal clear.

[@thatkidcalledchar]

Only three minutes in but you are already saving my coursework! Thank you!

[@thatrnerd4265]

I love hearing that! Glad it helps!

[@jessicawilliams-daley372]

Thank you so much for this helpful video! I was really struggling but you came to my aid. Keep up your amazing work! 💖

[@RayaneJaffal]

Thank you so much for this simple and well explained tutorial!!

[@michaelkuffour5158]

a shorter way : team

[@Okwach_Kich]

You are a god

[@yaweli2968]

@@Okwach_Kich : lol he is not, it’s practice😂

[@suptibhattacharyya6556]

Error in plot.window(...) : need finite 'xlim' values In addition: Warning messages: 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion 2: In min(x) : no non-missing arguments to min; returning Inf 3: In max(x) : no non-missing arguments to max; returning -Inf When I am going for plot function, this error is appearing, could you help?

[@hebeboerhout4793]

Hey! I have the same problem, did you find a solution? Or where the error is?

[@Sjaapdespaak]

@@hebeboerhout4793 I found the fix, plot(factor(team), weight) R can't plot characters on the X-as without using the factor

[@Seppiscool]

@@Sjaapdespaak Thanks!! You are my hero!

[@hannawood803]

YOU HAVE SAVED MY LIFE. Thanks man!!!

[@thatrnerd4265]

You're welcome!

[@Karen-zz7qi]

ThANK YOU!! My teacher is using Minitab and doesn’t explain the what the ANOVA button does at all. I’m trying to learn r instead of Minitab this semester

[@EdvanAguiarEA]

Would you perform a paired ANOVA? Like, imagine that you intend to compare the same team, based on the last three seasons. I appreciate it!

[@farazshaikh4967]

Edvan Aguiar i would like tht too

[@MrGyanu123]

Thanks for making a nice video. Could you please make a video on " how to prepare table denoting with mean value with different letters ( a, b, c, d, ab, bc, and ad )for examples 1.03±0.08a to mark significant differences in a table.

[@ramdasmenon1965]

I think you've missed the phrase "fail to"; the p value is 0.189, so we FAIL TO reject the null hypothesis which states that there is no significant difference between the means. Great video though!

[@saveragamashoka6532]

What does the DF stand for in the summary ?

[@marcopozzan445]

THANK YOU !! You saved my FRM course

[@sofianieto7345]

Thank you! your video was really clearly and helpful to do some things about my thesis

[@NigaJay]

Thanks man. U saved me from a college work

[@zinimbita9789]

Thank you for your excellent explanations

[@marlenethaitumu8177]

Amazing work!!!!

[@johnsonmshiu4837]

Great job ..may you please make a tutorial regarding multnormal test in R

[@arelahammond3987]

this video is a lifesaver!! thank you so much!

[@zapy422]

Great video. If there was a difference in the means how, to know with team has the highest mean?

[@thatrnerd4265]

I would use the dplyr library (if I understand your question) library(dplyr) df %>% group_by(team) %>% summarize(mean_weight = mean(weight, na.rm = TRUE)) %>% arrange(desc(weight)) This will get the average weight for each team and arrange them highest to lowest by average. I have a video here you can check out that talks about it :) - kzbin.info/www/bejne/eovOep5tqNStr8U

[@admiraqil1014]

u saved my butt, thanks a lot

[@diahidvegi8536]

This is very nice. I just wish you had actually done the transformations to fix for normality, that's actually why I came to watch the video :/ :)

[@santiagodelacruzz]

Thank you very much! It really helped me

[@rcp3758]

why do i keep getting Error in plot.window(...) : need finite 'xlim' values while plotting?

[@cooljayhawk8774]

yup

[@sisaymarkoskahsay899]

Very interesting concept but how we can give lettering for the groups when there is significant variation among them

[@peterkalinde]

How can I change the p-value in the summary

[@jadmelhem3827]

YOU ARE A GOD

[@MrEkballo]

@thatRnerd, great video! I know this is about football teams, but if I was comparing how some substance affects the growth of bacteria and I would be given concentration just in words "Low, Medium and High" , I would have to substitute these words with 1 (for low), 2 (medium), 3 (high) , right? Or maybe 3,6,9 or 0,5,10? How to choose a right number to substitute words?

[@thatrnerd4265]

Because it is a class grouping, it is not going to effect the anova output for whichever values you pick. I would probably just go 1, 2, 3 or just use the words low, medium and high.

[@vaidehichennubhotla7887]

great video! Thanks ! Please tell me if this is all there is in a DOE.

[@thatrnerd4265]

For DOE you are often looking at multiple variables, and combinations between variables. In our example we are only considering weight, only one variable. So, there are often extra things you must do for more complex DOE.

[@princesscfn]

I keep getting this error accuracy.aov summary(accuracy.aov) Error in summary(accuracy.aov) : object 'accuracy.aov' not found i have followed all steps of the tutorial..

[@thatrnerd4265]

You need the 'tilde' button instead of the dash. It is the button above tab while holding shift. ~ It's kind of hard to see that in the video :)

[@princesscfn]

@@thatrnerd4265 LIFE SAVER !!! Thank you sooo much it's finally worked

[@abdouazzi405]

thank u! i wanna just more explanation about the table of aev

[@renanlopes2785]

I got an error when trying to plot this data : need finite 'xlim' values

[@cooljayhawk8774]

same

[@darthszarych5588]

Thanks this was helpful

[@chiragramoliya2580]

hello man, i need your help, for how to use diallelanalysisR package in R studio

[@akankshanikam2416]

Thank you for this. its really usefull!

[@milan9191]

When I try to plot and do anova I get an error saying "non-numergic argument to a binary operator" any ideas?

[@Sasha-bg2eg]

Just wanted to clarify if this is one way or two way ANOVA? Thankyou

[@thatrnerd4265]

Yes, this is a one-way ANOVA.

[@bryanlester1001]

When I use $, for example > weight

[@rockydelrio]

It should be a data frame for you to use $ sign

[@arielserravalle7909]

This error is caused by creating the linear model with VECTORS or DATA e.g. Vectors lin

[@sa-wi3lo]

How about preparing the data before starting the tutorial

[@TheBarOst]

would be easier: boxplot(data$weight ~ data$team) where data is our dataset

[@henrykelderman]

indeed as plot does not accept character data from team for me (version 4.0.4)

[@Nady2022]

How do i input Dacan Multiple Range Test (DMRT) command for means separation using R studio

[@mahmoudelnakeeb6108]

Plot coding didn't work with me, could you help me plz.

[@surbhiagrawal3951]

can we use melt(football) for combing the data?

[@lindasturm-flores6752]

Very helpful. Thank you!

[@calmyourself8380]

Thank you, Really appreciate it

[@bridgettsmith7206]

Thank you for the video

[@RabinRaut-h5v]

plot(weight~team, data=df) I am getting error in this command.

[@rabinraut1408]

I am liking you very video. Please help me sort out this issue. Thanks plot(weight ~ team, data = df) Error in plot.window(...) : need finite 'xlim' values In addition: Warning messages: 1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion 2: In min(x) : no non-missing arguments to min; returning Inf 3: In max(x) : no non-missing arguments to max; returning -Inf

[@chaimabenhmida3791]

This example is one way or two way?

[@slowburn3r949]

Does this method yield a Type III sum of squares?

[@thabisomasisi6426]

i don't see a link for the data set

[@alibaluee581]

Do you have any video for unbalanced design?

[@AnandKumar-dn2yd]

Hi I need ANOVA table for augumented RBD, can anyone help me in providing the code

[@yaweli2968]

Amazing video, thanks is the least I can say.

[@thatrnerd4265]

You're welcome! Thank you for the kind comment!

[@Daz2281]

Thank you!

[@omphileashley3142]

Appreciated.

[@stevenwilson5556]

R studio is awesome, but for the love of God, use an editor theme that's not pure eye bleeding painful white.

[@nathannickelson4353]

Thank god ... jeez this was easy peezy.

[@thatrnerd4265]

Glad this helped make it easier! R is amazing!

[@elpasemah]

thank you

[@sarithasanthashekharan7136]

Error found as could not find function team

[@acharyaaviseka85]

HI Thanks for the video and its a very useful but I have my data which has samples are in replicates and I have 50k row . how should I proceed to analyze and plot. kindly help. Thanks in Advance :) Avi

[@thatrnerd4265]

If the plot you are analyzing is a boxplot, then it is analyzed the same way. There will likely be more outliers with that many observations, but it still means the same thing. Middle line is median top and bottom of box is 25th and 75th percentiles.

[@rushikeshbulbule8120]

thanks

[@victoriabydone9840]

understood well

[@kidrenray8657]

Now all I need to do is understand everything else I’m missing, haha

[@paryazareie4935]

anova

[@EBSims]

you are the only person who deserves to go to heaven

[@hmmmmm8716]

I love you 💋

[@arsalansyed4709]

thank you!