groupmeans

Load the fisheriris data set.

load fisheriris

The column vector species contains iris flowers of three different species: setosa, versicolor, and virginica. The matrix meas contains four types of measurements for the flower: the length and width of sepals and petals in centimeters.

Perform a one-way MANOVA with species as the factor and the measurements in meas as the response variables.

maov = manova(species,meas)

maov = 
1-way manova

Y1,Y2,Y3,Y4 ~ 1 + Factor1

    Source     DF     TestStatistic    Value       F       DFNumerator    DFDenominator      pValue  
    _______    ___    _____________    ______    ______    ___________    _____________    __________

    Factor1      2       pillai        1.1919    53.466          8             290         9.7422e-53
    Error      147                                                                                   
    Total      149                                                                                   


  Properties, Methods

maov is a manova object that contains the results of the one-way MANOVA. The small p-value for species indicates that the flower species has a statistically significant effect on at least one of the flower measurements.

Calculate the marginal means for maov.

means = groupmeans(maov)

means=3×5 table
       Factor1         Mean        SE       Lower     Upper 
    ______________    ______    ________    ______    ______

    {'setosa'    }    2.5355    0.042807    2.4509    2.6201
    {'versicolor'}     3.573    0.042807    3.4884    3.6576
    {'virginica' }     4.285    0.042807    4.2004    4.3696

The 95% confidence intervals in the table output do not overlap, which is consistent with the small p-value for species.

Load the patients data set.

load patients

The variables Systolic and Diastolic contain data for patient systolic and diastolic blood pressure. The variables Smoker and SelfAssessedHealthStatus contain data for patient smoking status and self-assessed heath status.

Use the table function to create a table of factor values from the data in Systolic, Diastolic, Smoker, and SelfAssessedHealthStatus.

tbl = table(Systolic,Diastolic,Smoker,SelfAssessedHealthStatus,VariableNames=["Systolic" "Diastolic" "Smoker" "SelfAssessed"]);

Perform a two-way MANOVA to test the null hypothesis that smoking status does not have a statistically significant effect on systolic and diastolic blood pressure, and the null hypothesis that self-assessed health status does not have an effect on systolic and diastolic blood pressure.

maov = manova(tbl,["Systolic" "Diastolic"])

maov = 
2-way manova

Systolic,Diastolic ~ 1 + Smoker + SelfAssessed

       Source       DF    TestStatistic     Value         F       DFNumerator    DFDenominator      pValue  
    ____________    __    _____________    ________    _______    ___________    _____________    __________

    Smoker           1       pillai         0.67917     99.494          2              94         6.2384e-24
    SelfAssessed     3       pillai        0.053808    0.87552          6             190            0.51392
    Error           95                                                                                      
    Total           99                                                                                      


  Properties, Methods

maov is a manova object that contains the results of the two-way MANOVA. The small p-value for the Smoker term in the MANOVA model indicates that enough evidence exists to conclude that mean response vectors are statistically different across the factor values of Smoker. However, the large p-value for the SelfAssessed term indicates that not enough evidence exists to reject the null hypothesis that the mean response vectors are statistically the same across the values for SelfAssessed.

Calculate the marginal means for the values of the factor Smoker.

groupmeans(maov,"Smoker")

ans=2×5 table
    Smoker     Mean       SE       Lower     Upper 
    ______    ______    _______    ______    ______

    false     99.203    0.45685    98.296    100.11
    true      109.45    0.62574    108.21     110.7

The output shows that the marginal mean for non-smokers is lower than the marginal mean for smokers.

Load the carsmall data set.

load carsmall

The variable Model_Year contains data for the year a car was manufactured, and the variable Cylinders contains data for the number of engine cylinders in the car. The Acceleration and Displacement variables contain data for car acceleration and displacement.

Use the table function to create a table of factor values from the data in Model_Year, Cylinders, Acceleration, and Displacement.

tbl = table(Model_Year,Cylinders,Acceleration,Displacement,VariableNames=["Year" "Cylinders" "Acceleration" "Displacement"]);

Perform a two-way MANOVA using the table variables Year and Cylinders as factors, and the Acceleration and Displacement variables as response variables.

maov = manova(tbl,"Acceleration,Displacement ~ Cylinders + Year")

maov = 
2-way manova

Acceleration,Displacement ~ 1 + Year + Cylinders

     Source      DF    TestStatistic     Value        F       DFNumerator    DFDenominator      pValue  
    _________    __    _____________    ________    ______    ___________    _____________    __________

    Year          2       pillai        0.084893    2.1056          4             190           0.081708
    Cylinders     2       pillai         0.94174     42.27          4             190         2.5049e-25
    Error        95                                                                                     
    Total        99                                                                                     


  Properties, Methods

maov is a manova object that contains the results of the two-way MANOVA. The table output shows that the p-value for the MANOVA model term Year is too large to conclude that Year has a statistically significant effect on Acceleration and Displacement.

Calculate the means, standard errors, and 95% confidence intervals for each response variable, grouped by the factor values in Year.

means = groupmeans(maov,"Year",MarginalMean=false)

means=6×6 table
    Year     ResponseName      Mean       SE       Lower     Upper 
    ____    ______________    ______    _______    ______    ______

     70     "Acceleration"    14.102    0.45499    13.199    15.005
     70     "Displacement"    233.02     6.4983    220.12    245.92
     76     "Acceleration"    15.623    0.41203    14.805    16.441
     76     "Displacement"    218.43     5.8848    206.75    230.12
     82     "Acceleration"    15.083    0.53717    14.016    16.149
     82     "Displacement"    234.09     7.6721    218.86    249.32

The table output shows that the 95% confidence intervals have significant overlap for each response variable, supporting the conclusion that Year does not have a statistically significant effect on Acceleration and Displacement.