How to correctly specify WithinDesign for fitrm
16 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
I am trying to use matlab function fitrm. But the input 'WithinDesign' does not seem to make any sense to me. It seems as anything with the same length as the number of repeated measures is acceptable and the final output seem independent of the this WithinDesign. Take for instance:
load repeatedmeas
rm = fitrm(between,'y1-y8 ~ Group*Gender','WithinDesign',within)
ranova(rm)
and
load repeatedmeas
rm = fitrm(between,'y1-y8 ~ Group*Gender','WithinDesign',[1:8])
ranova(rm)
and
load repeatedmeas
rm = fitrm(between,'y1-y8 ~ Group*Gender','WithinDesign',zeros(8,1))
ranova(rm)
all produces the exact same output, as shown in the attached image.
Could someone please explain to me why they are all the same. Since its all the same what is the point of specifying it?. Also could some explain the design matrix 'within', What is it really telling the program to do, any reference to SPSS would be great.
Thanks.
0 commentaires
Réponse acceptée
Jeff Miller
le 29 Mar 2018
WithinDesign lets you store specific numeric values associated with the levels of the within-Ss factors in the RepeatedMeasuresModel object that fitlm produces as output. For example, these values might be the exact dosages of some drug if you were comparing 8 dosages. As you have noted, these don't affect the ANOVA output per se, but they can affect subsequent things you do with the output object.
For example, you get different plots with:
rm8 = fitrm(between,'y1-y8 ~ Group*Gender','WithinDesign',[1:8])
plot(rm8)
rm0 = fitrm(between,'y1-y8 ~ Group*Gender','WithinDesign',zeros(8,1))
plot(rm0)
The ANOVA model doesn't care what numerical values are associated with those 8 dosages, but it is still nice to show them in the plot. As another example, you would also need the dosage values if you wanted to compute a linear trend across the 8 dosages.
As far as I know, SPSS doesn't integrate the X values into its repeated measures ANOVAs in this way (but I don't know SPSS that well).
2 commentaires
Thomas61197
le 31 Mar 2020
if y1-y8 are variables containing the response values to different dosages then 'y1-y8 ~ Group*Gender' means the response to the different dosages depends on the independent variables Group and Gender. If you would like to know what the effect is of the different dosages on the response, the model specification would look like this: 'response ~ dosage' where response is a dependent variable containing numerical values and where dosage is a independent catagorical variable containing the 8 dosages, y1-y8. Here 'between' would be a table with the amount of rows corrosponding to the number of obervations and containing the two variables response and dosage.
My question is: what do i put here?
rm = fitrm(between,'response ~ dosage','WithinDesign', here)
Plus de réponses (1)
Ahmed Ramadan
le 26 Juil 2018
To correctly specify within design, you need to modify your ranova command as follows.
ranovatbl = ranova(rm,'WithinModel','w1+w2')
On that page, you'll see different p-values for the two cases:
- ranovatbl = ranova(rm)
- [ranovatbl,A,C,D] = ranova(rm,'WithinModel','w1+w2')
1 commentaire
Adam Danz
le 22 Mai 2019
That only specifies the within-subject variable names. The question was regarding the within-subject model input to fitrm(). An example of that is here:
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!