User Tools

Site Tools


c:ms:2017:schedule:week15

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
c:ms:2017:schedule:week15 [2017/06/14 09:09] hkimscilc:ms:2017:schedule:week15 [2017/06/15 09:02] (current) hkimscil
Line 1: Line 1:
-====== t-test ====== +https://goo.gl/forms/x7NUmOxm2DxKmpd62 
-<code><- rnorm(10, mean=100, sd=15) + 
-x +====== two sample t-test ====== 
-t.test(x, mu=95)+<code>x1 <- rnorm(10, mean=100, sd=15) 
 +x1 
 +t.test(x1, mu=95)
 </code> </code>
 +  - What is the std error value?
 +  - What if you have a sample with a large size?
 +    * How do you get your se?
 +
 +<code>ster <- sd(x1)/sqrt(length(x1))
 +difx <- mean(x1)-95
 +tvalue <- difx/ster
 +
 +difx
 +ster
 +tvalue
 +
 +</code>
 +
 +<code>x2 <- rnorm(100, mean=100, sd=15)
 +head(x2)
 +tail(x2)
 +t.test(x2, mu=95)
 +
 +ster <- sd(x2)/sqrt(length(x2))
 +difx <- mean(x2)-95
 +tvalue <- difx/ster
 +
 +difx
 +ster
 +tvalue
 +
 +</code>
 +
 +<code>mtcars
 +mtcars$mpg
 +mtcars$am
 +
 +t.test(mtcars$mpg~mtcars$am)
 +</code>
 +
 +====== one sample t-test 2 ======
 +
 +<code>a = c(65, 78, 88, 55, 48, 95, 66, 57, 79, 81)
 +mean(a)
 +sd(a)
 +
 +t.test(a, mu=60)
 +</code>
 +====== Faired sample t-test ======
 +
 +<code>sleep
 +
 +sleep.wide <- data.frame(
 +    ID=1:10,
 +    group1=sleep$extra[1:10],
 +    group2=sleep$extra[11:20]
 +)
 +sleep.wide
 +
 +t.test(sleep.wide$group1 - sleep.wide$group2, mu=0, var.equal=TRUE)
 +t.test(sleep.wide$group1, sleep.wide$group2, paired=TRUE)
 +</code>
 +
 +{{:survey_ms_2017.sav}}
 +----
 +gender: male vs. female
 +
 +<code>T-TEST GROUPS=tend(1 2)
 +  /MISSING=ANALYSIS
 +  /VARIABLES=bigco smallco likeprog
 +  /CRITERIA=CI(.95).
 +</code>
 +
 +<code>T-TEST GROUPS=tend(1 2)
 +  /MISSING=ANALYSIS
 +  /VARIABLES=selfprom ethalright
 +  /CRITERIA=CI(.95).
 +</code>
 +
 +----
 +tendency - hum/soc vs. natsc/eng
 +
 +<code>T-TEST GROUPS=gender(1 2)
 +  /MISSING=ANALYSIS
 +  /VARIABLES=bigco smallco likeprog
 +  /CRITERIA=CI(.95).
 +</code>
 +
 +<code>T-TEST GROUPS=gender(1 2)
 +  /MISSING=ANALYSIS
 +  /VARIABLES=selfprom ethalright
 +  /CRITERIA=CI(.95).
 +</code>
 +====== ANOVA ======
 +[[:r/oneway_anova|Oneway ANOVA]]
 +[[:r/twoway_anova|Twoway ANOVA]]
 +====== Regression ======
 +[[:/r/linear_regression|Linear Regression]]
 +[[:r/multiple_regression|Multiple Regression]]
 +
 +====== class survey ======
 +
 +<code>
 +UNIANOVA bigco BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +
 +UNIANOVA smallco BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +
 +UNIANOVA allowance BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +
 +UNIANOVA workhr BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +
 +UNIANOVA selfprom BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +
 +UNIANOVA ethalright BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +
 +UNIANOVA phmin BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +
 +UNIANOVA phpeo BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +
 +UNIANOVA textmin BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +
 +UNIANOVA textpeo BY gender year2
 +  /METHOD=SSTYPE(3)
 +  /INTERCEPT=INCLUDE
 +  /PLOT=PROFILE(gender*year2)
 +  /EMMEANS=TABLES(gender) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(year2) COMPARE ADJ(LSD)
 +  /EMMEANS=TABLES(gender*year2) 
 +  /PRINT=OPOWER ETASQ DESCRIPTIVE
 +  /CRITERIA=ALPHA(.05)
 +  /DESIGN=gender year2 gender*year2.
 +</code>
 +
  
c/ms/2017/schedule/week15.1497400785.txt.gz · Last modified: 2017/06/14 09:09 by hkimscil

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki