First of all let me began by thanking you for giving such a wonderful presentation on SPD in r and sir I need some clarity that while comparing interaction treatments means in SPD, usually we will look in two directions know sir, i.e., (1) two sub-plot treatment means at a given main-plot treatment and (2) two main plot treatment means at same or different sub-plot treatment, do we get in that way of interaction comparisons in this package?. Please do respond sir.
Thank you for your compliment. The above two points which you have mentioned specify the main effects for both factors whatever I perceived from your statement. Interaction effect in split plot specify mean values for both main plot and sub plot factor however the difference from factorial experiment is that the factor with less importance is kept in main plot with containing larger experimental error.
@@AGRONInfoTech I agree with Mr. devadas. there are two different LSD value to calculate while we are comparing interaction treatments means in SPD, (1) comparision. To compare two subplot treatment means for the same whole plot treatment. (e.g. super:4w v.s. super:3w). (2) comparision. To compare two whole plot means at the same or different sub plot treatments (e.g. shaheen:4w v.s. super:2w) or (e.g. shaheen:2w v.s. super:2w). please check the webiste for more detail LSD calculation of SPD: www.ndsu.edu/faculty/horsley/splplot.pdf As your R code in the video, it seems that the result of "out3" conformed to the (1) comparision. that means, the value of LSD (out3)=4.132411 Only can used in comparing subfactor"age" for the same mainfactor "var". It's Not suitable for all means of comparision. We need to calculate another LSD value for the (2) comparision. My question is, how to calculate the LSD value for the (2) situation??
@@AGRONInfoTech thank you so much. You're help is invaluable. I'm facing two other problems with my data. First I take mesuarements in 3 different times, so I have the following factors: factor light (4), factor variety (4), factor NaCl treatment (2), replications (4), and with time I have another factor, right?. So I have 4x4x2x4x3 which lead me with 384 observations 😮💨. I was thinking in analyzing each time separately to see what is the factor that influences the most, then eliminate that factor and evaluate within my 3 different times. What do you think of this approach? Also, I'm worried because I might and unbalanced design since some plants in my experiment die earlier. How does that influences in the analysis? (Sorry to bother you with all theses questions)
Thank you very much and it was great Experience, and I would like to know which CV value to use in split plot, the main plot factor, the sub plot factor or the interaction?, because it lists like CVa,and CVb ,I would like to thank you and wish to follow your other posts.
There are actually two CVs, one for the sub plot CV (b) and one for whole plots CV (a). Presentation of the former is usually sufficient. However you can also use the whole plot CV. CV for whole plots is usually lower than CV for sub plot because larger plot tend to have lower CV. Sub plot CV also reflects variation for interaction term.
OMG! thank you so much! its works perfectly and I understand clearly on this method! i have a question, what is the code for add the average number of each bar? I want to put each X bar of average on top of each error line.. thank you
Thanks. You can use text(data, model$average, cex=1) to add average values above the bars. Here model may be the output of aov() function and after dollar sign type the component representing average values.
Thanks for your helpful tutorial. Could you tell how to run LSD test for multi-variables data set? If I run LSD test one by one, tit will takes long time due to there are only four groups but 15 variables.
You can check the normality of the data before proceeding for this analysis. However for homogeneity of variances it depends what type of data you are working on. For example if you have a dataset that was repeated over years, location, sites or harvests then you need to first look at the homogeneity of variances across the repeating factor. I hope this make sense.
Very good demonstration. For the interaction plot, how do i sort them? For example, I want them in the order of Shaheen:2W, Super:2W, Shaheen:3W, Super:3W, Shaheen:4W, Super:4W.
Thank you. I shall create video on it later on however, if you need it urgently then you can send me data and the Rscript you have at the below email address. I shall amend it and will share it with you. agron.infotech@gmail.com
I have to work on augmented factorial and split anlysis. I shall keep augmented split plot analysis in loop. Shall upload video on this aspect later on. Thanks for the suggestion.
hello. I am trying, but I m still having the ANOVA table with "'Nan'" for block, Ea, and Eb in F value and Pr columns. Then, when I continued with the LSD test, I m still getting this message '"Error in qt(1 - alpha/2, DFerror) :Non-numeric argument to mathematical function". How can I solve the issue, please.
Thank you for nice presentation to learn the SPD analysis. Can you please share the R Codes for this analysis. Please let me know the R and R Studio difference. Can we done the Split Plot Design Analysis in both of this. Thanks.
Dear Dr. Baviskar the link for R script or commands to perform split plot analysis is added in description of this video. Please download it. However for the difference in R and R studio, I shall simply say as the difference between window 98 and XP. You need to install R and then R studio.
This is the model I would like to develop in R. I am trying to use the Location:Treatment interaction df and MSE in the LSD separation test but don't know how to do. Your R code showed the error message from my previous questions.Combined RCBD (treatments fixed, locations random). SOV Df Expected Mean Square F-test Location l-1 Non-valid Rep(Loc) l(r-1) Non-valid Treatment t-1 1222−++ΣtrlrLTτσσ T MS/Loc*Trt MS Location*Treatment (l-1)(t-1) 22LTrσσ+ Loc*Trt MS/Error MS Error l(r-1)(t-1) 2σ Total lrt-1
Can you send me data containing at least one variable with treatment, location and replications levels at agron.infotech@gmail.com I shall run the same model and see what actual problem you are facing and how to fix it.
When main or interaction effect is non significant then there is no need to go for mean comparison test however I have just performed it to show the way in case if you have significant interaction. Ea and Eb are the error terms they don't have probability values. Only rep, factor A, factor B and interaction contain p value.
Thanks for your feed back. I shall also include 3 and 4 factor factorial analysis soon. However you can go through 2 factor factorial analysis under rcbd and crd in videos. ru-vid.com/video/%D0%B2%D0%B8%D0%B4%D0%B5%D0%BE-ATAwYOM4liM.html
Use below command for split split plot analysis ssp.plot(block, pplot, splot, ssplot, Y) Where block = replictions pplot = Main plot factor splot = Sub plot factor ssplot = sub-subplot factor Y = response variable (e.g yield) You need to load library agricolae before using this command require(agricolae) Also use attach function before using ssp.plot function attach(datafilename)
In video tutorial follow instructions at 3:22 glb and Eb commands are used to represent Error df and Error MS for the interaction term. And then use out1 command to get LSD test. If you still did not get point let me know.
thank you for your video, but i still confused, some script is error symbol in Rstudio, start from "install.packages("agricolae" dependecies = TRUE)" and "sp.plot", please help if you have new script, thank you
Can you copy and send me complete commands that you used during analysis in R at agron.infotech@gmail.com I shall probably lookover what actual issue you are facing and also the data you have used for at least one response variable or yield.
May be you have not attached the components of the data file. Using below command may fix the it. attach(data) If still have some issues please let me know
@@AGRONInfoTech Thanks for your response. I have created a dummy csv file (drive.google.com/file/d/1aQuOODDLtkUI5K97f-HvXutxfX79v38Z/view?usp=sharing). Let's say the sequence I need is 20:a, 20:b:40:a,60:a,80:a,100:a,40:b,60:b,80:b,100:b.
@@AGRONInfoTech RStudio Desktop 2022.02.0+443. Prairie Trillium Release. Here is the error code - > library(agricolae) Error: package or namespace load failed for ‘agricolae’ in loadNamespace(i, c(lib.loc, .libPaths()), versionCheck = vI[[i]]): there is no package called ‘klaR’ In addition: Warning messages: 1: package ‘lmerTest’ was built under R version 3.6.2 2: package ‘lme4’ was built under R version 3.6.2 3: package ‘agricolae’ was built under R version 3.6.2
Have you successfully installed agricolae package keeping dependencies as TRUE? Try using below command in console to reinstall this package. install.packages("agricolae", dependencies = T)
@@AGRONInfoTech Thank you very much for your suggestion. I tried the command but still not working. Error: package or namespace load failed for ‘agricolae’ in loadNamespace(i, c(lib.loc, .libPaths()), versionCheck = vI[[i]]): there is no package called ‘klaR’ In addition: Warning message: package ‘agricolae’ was built under R version 3.6.2
Hi, I have two years of collected data from split-plot experiments. I am a new user of R, So, I am looking for help from anyone to combine the analysis of my experiment with full mainplot error and subplot error. Thank you!
Thank you for the Wonderful video!! I have some question, accordind to the website: www.ndsu.edu/faculty/horsley/splplot.pdf there are two different LSD value to calculate while we are comparing interaction treatments means in SPD, (1) comparision: To compare two subplot treatment means for the same whole plot treatment. (e.g. super:4w v.s. super:3w). (2) comparision: To compare two whole plot means at the same or different sub plot treatments (e.g. shaheen:4w v.s. super:2w) or (e.g. shaheen:2w v.s. super:2w). In this video, the result of "out3" conformed to the (1) comparision, that means, the value of LSD (out3)=4.132411 Only can used in comparing sub-factor"age" for the same main-factor "var". IT's Not Suitable for all interaction treatments means of comparision. We need to calculate another LSD value for the (2) comparision. Please let me know if i am wrong, and how to calculate the LSD value for the (2) situation?? Looking for your reply, thank you very much!!
Hello. I have tried testing with duncan.test for the interaction analysis but it didn't work. For my data, I typed as: duncan.test (data, response, main-factor:sub-factor, Edf_b, EMS_b, alpha=0.05, group=TRUE, console=TRUE), it returned an error as: Error in names(DUNCAN)
Before analysis of variance you need to convert factor variables from character to factor by using below command: data$main-factor = as.factor(data$main-factor) data$sub-factor = as.factor(data$sub-factor) Then run ANOVA followed by duncan.test Hope this will fix the issue. If still have an error share me your data at agron.infotech@gmail.com
Thank you for your reply. I followed what you has suggested but it showed the same error. It only work with LSD.test. I will send the data to your email. Thank you.