Saves results from an unbalanced analysis of variance, by `AUNBALANCED`

, in R data frames (R.W. Payne & Z. Zhang).

### Options

`TERM` = formula |
Treatment term whose means, effects etc. are to be saved; must be set if any of these are to be saved, unless there is only one treatment term |
---|---|

`COMBINATIONS` = string token |
Factor combinations for which to form predicted means (`present` , `estimable` ); default `esti` |

`ADJUSTMENT` = string token |
Type of adjustment to be made when predicting means (`marginal` , `equal` , `observed` ); default `marg` |

`LSDLEVEL` = scalar |
Significance level (%) for least significant differences and multiple comparisons; default 5 |

`RMETHOD` = string token |
Type of residuals to form (`simple` , `standardized` ); default `simp` |

`MCOMPARISON` = string token |
Method to use to make multiple comparisons between the means (`flsd` , `fstudentizedlsd` , `bonferroni` , `sidak` ); default `*` i.e. none |

`SAVE` = identifier |
Save structure (from `AUNBALANCED` ) containing details of the analysis for which further output is required; if omitted, output is from the most recent use of `AUNBALANCED` |

### Parameters

`INFORMATION` = string tokens |
What to save (`aovtable` , `effects` , `means` , `semeans` , `esemeans` , `sedmeans` , `lsdmeans` , `replications` , `fittedvalues` ); must be set |
---|---|

`OUTFILENAME` = texts |
Name of the R (`.rda` ) file to create for each set of information; must be set |

`COLUMNNAMES` = texts |
Specifies names for the columns in the file; if this is not set, suitable names are chosen automatically |

`EXIT` = scalars |
Records the exit status, 0 if the information was saved successfully, 1 otherwise |

### Description

`AU2RDA`

puts results from an analysis by `AUNBALANCED`

into R data frames. By default the results are from the most recent `AUNBALANCED`

analysis, but you can use the `SAVE`

option to specify the save structure from some other analysis.

The `INFORMATION`

parameter specifies the information to save, with settings:

`aovtable` |
analysis of variance table; |
---|---|

`effects` |
table of treatment effects; |

`means` |
table of predicted means; |

`semeans` |
table of standard errors of means; |

`esemeans` |
table of effective standard errors of means; |

`sedmeans` |
symmetric matrix of standard errors of differences of means; |

`lsdmeans` |
symmetric matrix of least significant differences of means; |

`replications` |
replication table; |

`fittedvalues` |
y-variate, fitted values and residuals. |

The `OUTFILENAME`

parameter specifies the name of the file to save the information. If this does not have a `.rda`

extension, one will be appended to the name automatically.

The `COLUMNNAMES`

parameter allows you to specify a text containing names to use for the columns in the file. Thus may be useful if you want to avoid name clashes when you are saving several sets of output.

The `EXIT`

parameter can save a scalar to record the exit status for each set of information. This contains zero if the information was saved successfully, and one otherwise. (For example, it will not be possible to save missing-value estimates if no responses were missing.)

The `TERM`

option specifies the treatment term whose means, effects etc. are to be saved. This must be set if any of these are to be saved, unless there is only one treatment term.

The `LSDLEVEL`

option specifies the significance level (%) to use in the calculation of least significant differences, or the experiment-wise significance level for multiple comparisons (default 5%).

The `MCOMPARISON`

option specifies the method to use for multiple comparisons. The settings are `FLSD`

(Fisher’s Least Significant Difference), `FSTUDENTIZEDLSD`

(Fisher’s Protected Least Significant Difference, using the Studentized Range statistic rather than Student’s t), `Bonferroni`

and `Sidak`

. By default, no multiple comparisons are done. The results are saved in a textual column, with the means. This contains identical letters (a, b etc.) alongside the sets of means that do not differ significantly. For more details see procedure `AMCOMPARISON`

, which is used to do the calculations.

The `RMETHOD`

option controls whether the residuals are simple residuals (like those printed by `ANOVA`

– the default) or whether they are standardized according to their variances.

Options: `TERM`

, `COMBINATIONS`

, `ADJUSTMENT`

, `LSDLEVEL`

, `RMETHOD`

, `MCOMPARISON`

, `SAVE`

.

Parameters: `INFORMATION`

, `OUTFILENAME`

, `COLUMNNAMES`

, `EXIT`

.

### Action with `RESTRICT`

If the `Y`

variate in the analysis is restricted, that restriction will carry over into the fitted-values spreadsheet.

### See also

Commands for: Analysis of variance and Input and output.

### Example

CAPTION 'AU2RDA example',\ 'Data from Genstat 5 Release 1 Reference Manual, page 340.';\ STYLE=meta,plain FACTOR [NVALUES=36; LEVELS=3; VALUES=12(1...3)] Block FACTOR [NVALUES=36; LABELS=!t(baresoil,emerald,emergo)] Leachate & [LABELS=!t('1','1/4','1/16','1/64')] Dilution VARIATE [NVALUES=36] Nhatch,Nnohatch READ Leachate,Dilution,Nhatch,Nnohatch 1 2 109 318 3 4 54 350 3 1 * 415 2 2 783 212 3 3 652 1375 2 4 490 816 1 3 95 1219 2 1 1012 66 1 4 166 943 3 2 1059 313 1 1 257 1006 2 3 1058 234 2 4 507 1119 1 2 194 840 1 3 175 1707 1 1 326 609 3 4 142 980 2 3 286 230 3 2 546 313 2 2 * 301 2 1 2471 112 3 3 76 489 1 4 208 503 3 1 * 325 1 1 322 913 1 2 255 2246 3 2 1774 1446 2 2 999 193 2 4 388 1836 3 4 221 1800 1 3 220 1902 2 1 2821 187 3 1 1486 463 3 3 717 1473 1 4 143 941 2 3 968 550 : CALCULATE Logit%h = LOG(Nhatch/Nnohatch) BLOCKSTRUCTURE Block TREATMENTSTRUCTURE Leachate*Dilution AUNBALANCED [PSE=differences,alldifferences; FPROBABILITY=yes] Logit%h " save the information " AU2RDA [TERM=Leachate.Dilution] INFORMATION=\ aovtable,means,sedmeans,replication; OUTFILENAME=\ 'aovtable.rda','means.rda','sedmeans.rda','rep.rda' " show the contents of the file " FOR file='aovtable.rda','means.rda','sedmeans.rda','rep.rda' IMPORT file; ISAVE=cols PRINT cols[] DELETE [REDEFINE=yes] cols " tidy up by deleting file (delete the line below if you want to keep it) " FDELETE file ENDFOR