-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi. I have been using SAS for some time, and now I have discovered R. I am very happy with it, but I have not found out how to perform some of the multiple comparisons I was used to do in SAS. With the SAS/STAT, I generally used the MEANS (for comparison of arithmetic means) and the LSMEANS (for adjusted means) statements of the GLM procedure (I think it is equivalent to lm in R). They provided a lot of tests: LSD, Duncan, Tukey-Kramer, Bonferroni, Scheff??, SNK, etc. However, in R I have only discovered Tukey-HSD. I have searched for information about this, but I was not successful. I wonder if anybody knows where I could learn about this. I would like to use these tests in R, and also obtain the adjusted means like these produced by the LSMEANS statement. Thank you. Felipe - ----------------------------------oOo---------------------------------- Felipe Mart??nez-Pastor (BSc, PhD) Animal Reproduction and Obstetrics Veterinary Clinic Hospital 24071-Le??n (Spain) Phone: 987 291 430 / 987 291 000 + 5203 Fax: 987 295 203 e-mail: dbcfmp at unileon.es - --- This message has been electronically signed using GPG, a free system to sign and encrypt documents. If you want to learn more, visit http://www.gnupg.org -----BEGIN PGP SIGNATURE----- iEYEARECAAYFAkMpUVYACgkQWtdQtNzjBl4rywCePiSxJw8/N6HzdZ7C+YMHf2K6 YSYAniw5GSo0ihrt4+OabHJ4c2PKNHkp =aF38 -----END PGP SIGNATURE-----
Felipe <felipe <at> unileon.es> writes:> With the SAS/STAT, I generally used the MEANS (for comparison of > arithmetic means) and the LSMEANS (for adjusted means) statements of the > GLM procedure (I think it is equivalent to lm in R). They provided a lot > of tests: LSD, Duncan, Tukey-Kramer, Bonferroni, Scheff伱仼, SNK, etc. > However, in R I have only discovered Tukey-HSD.Package multcomp with the workhorse-function simint comes close to what you want. Dieter
-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Thank you, I think multcomp is very near to what I was looking for. However, I am still looking for a mean to obtain least-squares (adjusted) means and std. errors of these means, and performing comparisons among these means, as the LSMEANS do in SAS. I have read other messages and have looked at car, effects and Design manuals, but I am not sure if this is what I am looking for. Any clue? Felipe Dieter Menne wrote: | Felipe <felipe <at> unileon.es> writes: | | |>With the SAS/STAT, I generally used the MEANS (for comparison of |>arithmetic means) and the LSMEANS (for adjusted means) statements of the |>GLM procedure (I think it is equivalent to lm in R). They provided a lot |>of tests: LSD, Duncan, Tukey-Kramer, Bonferroni, Scheff伱仼, SNK, etc. |>However, in R I have only discovered Tukey-HSD. | | | Package multcomp with the workhorse-function simint comes close to what you | want. | | | Dieter | | __ -----BEGIN PGP SIGNATURE----- iEYEARECAAYFAkMpo6UACgkQWtdQtNzjBl72rwCcCuw5qxD1BWsensDI71RzhNgL MUcAnA0Iq4tfoKSr/ymIV1nEZHZijvLW =pavA -----END PGP SIGNATURE-----
Do RSiteSearch("lsmeans") and go from there.
Andy
> From: Felipe
>
> -----BEGIN PGP SIGNED MESSAGE-----
> Hash: SHA1
>
> Thank you, I think multcomp is very near to what I was looking for.
> However, I am still looking for a mean to obtain least-squares
> (adjusted) means and std. errors of these means, and performing
> comparisons among these means, as the LSMEANS do in SAS. I have read
> other messages and have looked at car, effects and Design
> manuals, but I
> am not sure if this is what I am looking for. Any clue?
>
> Felipe
>
> Dieter Menne wrote:
> | Felipe <felipe <at> unileon.es> writes:
> |
> |
> |>With the SAS/STAT, I generally used the MEANS (for comparison of
> |>arithmetic means) and the LSMEANS (for adjusted means)
> statements of the
> |>GLM procedure (I think it is equivalent to lm in R). They
> provided a lot
> |>of tests: LSD, Duncan, Tukey-Kramer, Bonferroni, Scheff??, SNK, etc.
> |>However, in R I have only discovered Tukey-HSD.
> |
> |
> | Package multcomp with the workhorse-function simint comes close to
> what you
> | want.
> |
> |
> | Dieter
> |
> | __
> -----BEGIN PGP SIGNATURE-----
>
> iEYEARECAAYFAkMpo6UACgkQWtdQtNzjBl72rwCcCuw5qxD1BWsensDI71RzhNgL
> MUcAnA0Iq4tfoKSr/ymIV1nEZHZijvLW
> =pavA
> -----END PGP SIGNATURE-----
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>