An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20111126/2a02c3fc/attachment.pl>
SPSS -> R
7 messages · Kristi Shoemaker, Frank E Harrell Jr, michael.weylandt at gmail.com (R. Michael Weylandt +2 more
If you know SPSS already why not learn R modeling syntax and do this yourself? If ALPHA(.05) implies that you are using stepwise variable selection note that this is an invalid statistical technique. Frank Kristi Shoemaker wrote
I'm an SPSS user trying to make the transition to R. Can someone help me translate the following SPSS code into R?: GLM Total_tp1 Total_tp2 WITH Age Sex ? /WSFACTOR=Time 2 Repeated ? /METHOD=SSTYPE(3) ? /CRITERIA=ALPHA(.05) ? /WSDESIGN= Time ? /DESIGN= Age Sex Age*Sex. Also. can anyone recommend any resources to help SPSS users learn to things in R? Thanks, -kristi [[alternative HTML version deleted]]
______________________________________________ R-help@ mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
----- Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/SPSS-R-tp4110995p4111006.html Sent from the R help mailing list archive at Nabble.com.
Perhaps this website and the associated book will be of help: http://r4stats.com/ Michael
On Nov 26, 2011, at 11:08 AM, Kristi Shoemaker <kristi.shoemaker at yahoo.com> wrote:
I'm an SPSS user trying to make the transition to R. Can someone help me translate the following SPSS code into R?: GLM Total_tp1 Total_tp2 WITH Age Sex ? /WSFACTOR=Time 2 Repeated ? /METHOD=SSTYPE(3) ? /CRITERIA=ALPHA(.05) ? /WSDESIGN= Time ? /DESIGN= Age Sex Age*Sex. Also. can anyone recommend any resources to help SPSS users learn to things in R? Thanks, -kristi [[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Dear Kristi, I assume that this is a repeated-measures ANOVA with one within-subjects factor (Time) and two between-subjects factors (Age and Sex, which are crossed). If Age is numeric, and not a factor, then the "type-III" tests that you requested don't test sensible hypotheses. In any event, if my guess is right about the design, then you can use the Anova() function in the car package for an equivalent analysis. See the repeated-measures example in ?Anova (for the O'Brien and Kaiser data). You've already had an answer to the more general question. I hope this helps, John -------------------------------- John Fox Senator William McMaster Professor of Social Statistics Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox
-----Original Message----- From: r-help-bounces at r-project.org [mailto:r-help-bounces at r- project.org] On Behalf Of Kristi Shoemaker Sent: November-26-11 11:08 AM To: r-help at r-project.org Subject: [R] SPSS -> R I'm an SPSS user trying to make the transition to R. Can someone help me translate the following SPSS code into R?: GLM Total_tp1 Total_tp2 WITH Age Sex /WSFACTOR=Time 2 Repeated /METHOD=SSTYPE(3) /CRITERIA=ALPHA(.05) /WSDESIGN= Time /DESIGN= Age Sex Age*Sex. Also. can anyone recommend any resources to help SPSS users learn to things in R? Thanks, -kristi [[alternative HTML version deleted]]
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20111127/482090b5/attachment.pl>
2 days later
At 16:41 27/11/2011, Kristi Shoemaker wrote:
Hi John, Your assumptions are correct and those examples were very helpful, thanks! I think I'm almost there, but I'm screwing up something with the within-subjects factor (the example has two, I only have one). See below. Why am I not seeing my within-subjects factors in the ANOVA report?
I am by no means an expert in this but my reading of the documentation suggests that the idata and idesign parameters are for multivariate models and that is not what you specified.
Dat2
Subj Age Sex Time HippV 1 s01 32.9 1 0w 6.50098 2 s01 32.9 1 6w 6.91793 3 s02 35.1 0 0w 7.32480 4 s02 35.1 0 6w 7.56012 5 s03 34.4 0 0w 6.51385 6 s03 34.4 0 6w 6.56875 9 s05 39.9 1 0w 6.92855 10 s05 39.9 1 6w 6.94926 11 s06 29.5 1 0w 6.99383 12 s06 29.5 1 6w 7.10568 13 s07 45.9 1 0w 6.94380 14 s07 45.9 1 6w 7.08190 15 s08 20.3 1 0w 7.76881 16 s08 20.3 1 6w 7.72725 17 s09 26.9 0 0w 5.37566 18 s09 26.9 0 6w 5.74887 21 s11 22.0 0 0w 7.12992 22 s11 22.0 0 6w 7.16237 23 s12 31.0 1 0w 6.70629 24 s12 31.0 1 6w 6.80872 25 s13 50.1 1 0w 7.22649 26 s13 50.1 1 6w 7.58900 27 s14 22.2 0 0w 5.97577 28 s14 22.2 0 6w 5.80801 29 s16 20.4 1 0w 7.99554 30 s16 20.4 1 6w 8.09260 31 s20 24.0 0 0w 7.01014 32 s20 24.0 0 6w 6.87821 33 s21 32.4 0 0w 5.90883 34 s21 32.4 0 6w 5.95392 37 s24 22.2 0 0w 6.15474 38 s24 22.2 0 6w 6.00906 41 s27 22.2 1 0w 7.88765 42 s27 22.2 1 6w 7.76038 49 s35 24.3 0 0w 6.05998 50 s35 24.3 0 6w 6.07399 51 s36 23.5 0 0w 7.83182 52 s36 23.5 0 6w 7.60268 53 s38 59.7 1 0w 7.39672 54 s38 59.7 1 6w 6.98291 55 s39 40.5 0 0w 7.31330 56 s39 40.5 0 6w 7.50559 57 s40 24.2 1 0w 8.54958 58 s40 24.2 1 6w 8.65016 59 s41 23.6 1 0w 7.76049 60 s41 23.6 1 6w 7.58946 61 s42 53.3 0 0w 7.03388 62 s42 53.3 0 6w 7.48384 63 s43 34.4 0 0w 6.86967 64 s43 34.4 0 6w 6.81076 65 s44 44.8 0 0w 7.33779 66 s44 44.8 0 6w 7.86175 67 s45 40.2 1 0w 6.55963 68 s45 40.2 1 6w 6.52577 71 s47 26.5 0 0w 7.09418 72 s47 26.5 0 6w 6.92850 73 s48 22.7 0 0w 6.77078 74 s48 22.7 0 6w 6.67289 75 s50 36.1 1 0w 7.47208 76 s50 36.1 1 6w 7.55876
library(car)
Loading required package: MASS Loading required package: nnet Loading required package: survival Loading required package: splines
timepoints <- factor(c("0w", "6w"),
+ levels=c("0w","6w"))
idata=data.frame(timepoints) mod.ok <- lm(HippV ~ Age*Sex,data=Dat2) (av.ok <- Anova(mod.ok, idata=idata, idesign=~timepoints))
Anova Table (Type II tests)
Response: HippV
Sum Sq Df F value Pr(>F)
Age 0.0234 1 0.071 0.7909202
Sex 5.3082 1 16.128 0.0001781 ***
Age:Sex 4.7772 1 14.514 0.0003478 ***
Residuals 18.4314 56
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
________________________________ From: John Fox <jfox at mcmaster.ca> Cc: r-help at r-project.org Sent: Saturday, November 26, 2011 7:19 PM Subject: RE: [R] SPSS -> R Dear Kristi, I assume that this is a repeated-measures ANOVA with one within-subjects factor (Time) and two between-subjects factors (Age and Sex, which are crossed). If Age is numeric, and not a factor, then the "type-III" tests that you requested don't test sensible hypotheses. In any event, if my guess is right about the design, then you can use the Anova() function in the car package for an equivalent analysis. See the repeated-measures example in ?Anova (for the O'Brien and Kaiser data). You've already had an answer to the more general question. I hope this helps, John -------------------------------- John Fox Senator William McMaster Professor of Social Statistics Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox -----Original Message----- From: r-help-bounces at r-project.org [mailto:r-help-bounces at r- project.org] On Behalf Of Kristi Shoemaker Sent: November-26-11 11:08 AM To: r-help at r-project.org Subject: [R] SPSS -> R I'm an SPSS user trying to make the transition to R. Can someone help me translate the following SPSS code into R?: GLM Total_tp1 Total_tp2 WITH Age Sex /WSFACTOR=Time 2 Repeated /METHOD=SSTYPE(3) /CRITERIA=ALPHA(.05) /WSDESIGN= Time /DESIGN= Age Sex Age*Sex. Also. can anyone recommend any resources to help SPSS users learn to things in R? Thanks, -kristi [[alternative HTML version deleted]] [[alternative HTML version deleted]]
Michael Dewey info at aghmed.fsnet.co.uk http://www.aghmed.fsnet.co.uk/home.html
Dear Michael, The original poster and I ended up pursuing this question off-list and, I think, resolving the issue, discovering in the process that SPSS apparently reports "type-II" tests that are really "type-III" tests. (Yes, one fits a multivariate linear model prior to calling Anova() with the idata and idesign arguments.) Best, John -------------------------------- John Fox Senator William McMaster Professor of Social Statistics Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox
-----Original Message----- From: r-help-bounces at r-project.org [mailto:r-help-bounces at r- project.org] On Behalf Of Michael Dewey Sent: November-30-11 6:49 AM To: Kristi Shoemaker; r-help at r-project.org Subject: Re: [R] SPSS -> R At 16:41 27/11/2011, Kristi Shoemaker wrote:
Hi John, Your assumptions are correct and those examples were very helpful,
thanks!
I think I'm almost there, but I'm screwing up something with the within-subjects factor (the example has two, I only have one). See below. Why am I not seeing my within-subjects factors in the ANOVA
report? I am by no means an expert in this but my reading of the documentation suggests that the idata and idesign parameters are for multivariate models and that is not what you specified.
Dat2
Subj Age Sex Time HippV 1 s01 32.9 1 0w 6.50098 2 s01 32.9 1 6w 6.91793 3 s02 35.1 0 0w 7.32480 4 s02 35.1 0 6w 7.56012 5 s03 34.4 0 0w 6.51385 6 s03 34.4 0 6w 6.56875 9 s05 39.9 1 0w 6.92855 10 s05 39.9 1 6w 6.94926 11 s06 29.5 1 0w 6.99383 12 s06 29.5 1 6w 7.10568 13 s07 45.9 1 0w 6.94380 14 s07 45.9 1 6w 7.08190 15 s08 20.3 1 0w 7.76881 16 s08 20.3 1 6w 7.72725 17 s09 26.9 0 0w 5.37566 18 s09 26.9 0 6w 5.74887 21 s11 22.0 0 0w 7.12992 22 s11 22.0 0 6w 7.16237 23 s12 31.0 1 0w 6.70629 24 s12 31.0 1 6w 6.80872 25 s13 50.1 1 0w 7.22649 26 s13 50.1 1 6w 7.58900 27 s14 22.2 0 0w 5.97577 28 s14 22.2 0 6w 5.80801 29 s16 20.4 1 0w 7.99554 30 s16 20.4 1 6w 8.09260 31 s20 24.0 0 0w 7.01014 32 s20 24.0 0 6w 6.87821 33 s21 32.4 0 0w 5.90883 34 s21 32.4 0 6w 5.95392 37 s24 22.2 0 0w 6.15474 38 s24 22.2 0 6w 6.00906 41 s27 22.2 1 0w 7.88765 42 s27 22.2 1 6w 7.76038 49 s35 24.3 0 0w 6.05998 50 s35 24.3 0 6w 6.07399 51 s36 23.5 0 0w 7.83182 52 s36 23.5 0 6w 7.60268 53 s38 59.7 1 0w 7.39672 54 s38 59.7 1 6w 6.98291 55 s39 40.5 0 0w 7.31330 56 s39 40.5 0 6w 7.50559 57 s40 24.2 1 0w 8.54958 58 s40 24.2 1 6w 8.65016 59 s41 23.6 1 0w 7.76049 60 s41 23.6 1 6w 7.58946 61 s42 53.3 0 0w 7.03388 62 s42 53.3 0 6w 7.48384 63 s43 34.4 0 0w 6.86967 64 s43 34.4 0 6w 6.81076 65 s44 44.8 0 0w 7.33779 66 s44 44.8 0 6w 7.86175 67 s45 40.2 1 0w 6.55963 68 s45 40.2 1 6w 6.52577 71 s47 26.5 0 0w 7.09418 72 s47 26.5 0 6w 6.92850 73 s48 22.7 0 0w 6.77078 74 s48 22.7 0 6w 6.67289 75 s50 36.1 1 0w 7.47208 76 s50 36.1 1 6w 7.55876
library(car)
Loading required package: MASS Loading required package: nnet Loading required package: survival Loading required package: splines
timepoints <- factor(c("0w", "6w"),
+ levels=c("0w","6w"))
idata=data.frame(timepoints) mod.ok <- lm(HippV ~ Age*Sex,data=Dat2) (av.ok <- Anova(mod.ok, idata=idata, idesign=~timepoints))
Anova Table (Type II tests)
Response: HippV
Sum Sq Df F value Pr(>F)
Age 0.0234 1 0.071 0.7909202
Sex 5.3082 1 16.128 0.0001781 ***
Age:Sex 4.7772 1 14.514 0.0003478 ***
Residuals 18.4314 56
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
________________________________ From: John Fox <jfox at mcmaster.ca> Cc: r-help at r-project.org Sent: Saturday, November 26, 2011 7:19 PM Subject: RE: [R] SPSS -> R Dear Kristi, I assume that this is a repeated-measures ANOVA with one within-subjects factor (Time) and two between-subjects factors (Age
and
Sex, which are crossed). If Age is numeric, and not a factor, then the "type-III" tests that you requested don't test sensible hypotheses. In any event, if my guess is right about the design, then you can use the Anova() function in the car package for an equivalent analysis. See
the
repeated-measures example in ?Anova (for the O'Brien and Kaiser data). You've already had an answer to the more general question. I hope this helps, John -------------------------------- John Fox Senator William McMaster Professor of Social Statistics Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox
-----Original Message----- From: r-help-bounces at r-project.org [mailto:r-help-bounces at r- project.org] On Behalf Of Kristi Shoemaker Sent: November-26-11 11:08 AM To: r-help at r-project.org Subject: [R] SPSS -> R I'm an SPSS user trying to make the transition to R. Can someone help me translate the following SPSS code into R?: GLM Total_tp1 Total_tp2 WITH Age Sex /WSFACTOR=Time 2 Repeated /METHOD=SSTYPE(3) /CRITERIA=ALPHA(.05) /WSDESIGN= Time /DESIGN= Age Sex Age*Sex. Also. can anyone recommend any resources to help SPSS users learn
to
things in R?
Thanks,
-kristi
[[alternative HTML version deleted]]
[[alternative HTML version deleted]]
Michael Dewey info at aghmed.fsnet.co.uk http://www.aghmed.fsnet.co.uk/home.html
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting- guide.html and provide commented, minimal, self-contained, reproducible code.