docs/WithinSubjectsVignette.md
In cwendorf/easi-dev: Estimation Approach to Statistical Inference
title: "Estimation Approach to Statistical Inference"
author: "Craig A. Wendorf"
date: "2020-10-18"
output:
rmarkdown::html_vignette:
keep_md: TRUE
vignette: >
%\VignetteIndexEntry{Within Subjects Vignette}
%\VignetteEngine{knitr::rmarkdown}
%\VignetteEncoding{UTF-8}
Estimation Approach to Statistical Inference
Within Subjects Vignette
Three Time Period Example Data
Time1 <- c(5,6,6,7,8)
Time2 <- c(7,7,8,8,9)
Time3 <- c(8,8,9,9,9)
WithinData <- data.frame(Time1,Time2,Time3)
WithinData
## Time1 Time2 Time3
## 1 5 7 8
## 2 6 7 8
## 3 6 8 9
## 4 7 8 9
## 5 8 9 9
Analyses of Multiple Variables
Confidence Intervals for the Means
ciMeans(Time1,Time2,Time3)
## N M SD SE LL UL
## Time1 5 6.4 1.140 0.510 4.984 7.816
## Time2 5 7.8 0.837 0.374 6.761 8.839
## Time3 5 8.6 0.548 0.245 7.920 9.280
ciMeans(Time1,Time2,Time3,conf.level=.99)
## N M SD SE LL UL
## Time1 5 6.4 1.140 0.510 4.052 8.748
## Time2 5 7.8 0.837 0.374 6.077 9.523
## Time3 5 8.6 0.548 0.245 7.472 9.728
Plots of Confidence Intervals for the Means
cipMeans(Time1,Time2,Time3)
cipMeans(Time1,Time2,Time3,conf.level=.99,mu=6)
Significance Tests for the Means
nhstMeans(Time1,Time2,Time3)
## Diff SE t df p
## Time1 6.4 0.510 12.551 4 0
## Time2 7.8 0.374 20.846 4 0
## Time3 8.6 0.245 35.109 4 0
nhstMeans(Time1,Time2,Time3,mu=6)
## Diff SE t df p
## Time1 0.4 0.510 0.784 4 0.477
## Time2 1.8 0.374 4.811 4 0.009
## Time3 2.6 0.245 10.614 4 0.000
Analyses of a Variable Comparison
Confidence Interval for the Mean Difference
ciDifference(Time1,Time2)
## Diff SE df LL UL
## -1.400 0.245 4.000 -2.080 -0.720
ciDifference(Time1,Time2,conf.level=.99)
## Diff SE df LL UL
## -1.400 0.245 4.000 -2.080 -0.720
ciDifference(Time3,Time1)
## Diff SE df LL UL
## 2.200 0.374 4.000 1.161 3.239
Plots of Confidence Intervals for the Mean Difference
cipDifference(Time1,Time2)
cipDifference(Time1,Time2,conf.level=.99)
Significance Test for the Mean Difference
nhstDifference(Time1,Time2)
## Diff SE t df p
## -1.400 0.245 -5.715 4.000 0.005
nhstDifference(Time1,Time2,mu=-2)
## Diff SE t df p
## -1.400 0.245 2.449 4.000 0.070
Analyses of Pairwise Comparisons
Confidence Intervals for the Pairwise Comparisons
ciPairwise(Time1,Time2,Time3)
## Diff SE df LL UL
## Time1 v Time2 -1.4 0.245 4 -2.080 -0.720
## Time1 v Time3 -2.2 0.374 4 -3.239 -1.161
## Time2 v Time3 -0.8 0.200 4 -1.355 -0.245
ciPairwise(Time1,Time2,Time3,conf.level=.99)
## Diff SE df LL UL
## Time1 v Time2 -1.4 0.245 4 -2.528 -0.272
## Time1 v Time3 -2.2 0.374 4 -3.923 -0.477
## Time2 v Time3 -0.8 0.200 4 -1.721 0.121
Plot of the Confidence Intervals for the Pairwise Comparisons
cipPairwise(Time1,Time2,Time3)
cipPairwise(Time1,Time2,Time3,mu=-2,conf.level=.99)
Significance Tests of the Pairwise Comparisons
nhstPairwise(Time1,Time2,Time3)
## Diff SE t df p
## Time1 v Time2 -1.4 0.245 -5.715 4 0.005
## Time1 v Time3 -2.2 0.374 -5.880 4 0.004
## Time2 v Time3 -0.8 0.200 -4.000 4 0.016
nhstPairwise(Time1,Time2,Time3,mu=-2)
## Diff SE t df p
## Time1 v Time2 -1.4 0.245 2.449 4 0.070
## Time1 v Time3 -2.2 0.374 -0.535 4 0.621
## Time2 v Time3 -0.8 0.200 6.000 4 0.004
Analyses of a Set of Contrasts
Confidence Intervals for the Set of Contrasts
ciContrasts(Time1,Time2,Time3,contrasts=contr.sum)
## Est SE LL UL
## (Intercept) 7.6 0.115 7.334 7.866
## Variable1 -1.2 0.163 -1.577 -0.823
## Variable2 0.2 0.163 -0.177 0.577
ciContrasts(Time1,Time2,Time3,contrasts=contr.sum,conf.level=.99)
## Est SE LL UL
## (Intercept) 7.6 0.115 7.213 7.987
## Variable1 -1.2 0.163 -1.748 -0.652
## Variable2 0.2 0.163 -0.348 0.748
ciContrasts(Time1,Time2,Time3,contrasts=contr.treatment)
## Est SE LL UL
## (Intercept) 6.4 0.200 5.939 6.861
## Variable2 1.4 0.283 0.748 2.052
## Variable3 2.2 0.283 1.548 2.852
ciContrasts(Time1,Time2,Time3,contrasts=contr.poly)
## Est SE LL UL
## (Intercept) 7.600 0.115 7.334 7.866
## Variable.L 1.556 0.200 1.094 2.017
## Variable.Q -0.245 0.200 -0.706 0.216
ciContrasts(Time1,Time2,Time3,contrasts=contr.helmert)
## Est SE LL UL
## (Intercept) 7.6 0.115 7.334 7.866
## Variable1 0.7 0.141 0.374 1.026
## Variable2 0.5 0.082 0.312 0.688
ciContrasts(Time1,Time2,Time3,contrasts=contr.SAS)
## Est SE LL UL
## (Intercept) 8.6 0.200 8.139 9.061
## Variable1 -2.2 0.283 -2.852 -1.548
## Variable2 -0.8 0.283 -1.452 -0.148
Plot of the Confidence Intervals for the Set of Contrasts
cipContrasts(Time1,Time2,Time3,contrasts=contr.sum)
cipContrasts(Time1,Time2,Time3,contrasts=contr.sum,conf.level=.99)
cipContrasts(Time1,Time2,Time3,contrasts=contr.treatment)
cipContrasts(Time1,Time2,Time3,contrasts=contr.poly)
cipContrasts(Time1,Time2,Time3,contrasts=contr.helmert)
cipContrasts(Time1,Time2,Time3,contrasts=contr.SAS)
Significance Tests of the Set of Contrasts
nhstContrasts(Time1,Time2,Time3,contrasts=contr.sum)
## Diff SE t p
## (Intercept) 7.6 0.115 65.818 0.000
## Variable1 -1.2 0.163 -7.348 0.000
## Variable2 0.2 0.163 1.225 0.256
nhstContrasts(Time1,Time2,Time3,contrasts=contr.treatment)
## Diff SE t p
## (Intercept) 6.4 0.200 32.000 0.000
## Variable2 1.4 0.283 4.950 0.001
## Variable3 2.2 0.283 7.778 0.000
nhstContrasts(Time1,Time2,Time3,contrasts=contr.poly)
## Diff SE t p
## (Intercept) 7.600 0.115 65.818 0.000
## Variable.L 1.556 0.200 7.778 0.000
## Variable.Q -0.245 0.200 -1.225 0.256
nhstContrasts(Time1,Time2,Time3,contrasts=contr.helmert)
## Diff SE t p
## (Intercept) 7.6 0.115 65.818 0.000
## Variable1 0.7 0.141 4.950 0.001
## Variable2 0.5 0.082 6.124 0.000
nhstContrasts(Time1,Time2,Time3,contrasts=contr.SAS)
## Diff SE t p
## (Intercept) 8.6 0.200 43.000 0.000
## Variable1 -2.2 0.283 -7.778 0.000
## Variable2 -0.8 0.283 -2.828 0.022
cwendorf/easi-dev documentation built on Oct. 23, 2020, 4:28 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
title: "Estimation Approach to Statistical Inference" author: "Craig A. Wendorf" date: "2020-10-18" output: rmarkdown::html_vignette: keep_md: TRUE vignette: > %\VignetteIndexEntry{Within Subjects Vignette} %\VignetteEngine{knitr::rmarkdown} %\VignetteEncoding{UTF-8}
Estimation Approach to Statistical Inference
Within Subjects Vignette
Three Time Period Example Data
Time1 <- c(5,6,6,7,8)
Time2 <- c(7,7,8,8,9)
Time3 <- c(8,8,9,9,9)
WithinData <- data.frame(Time1,Time2,Time3)
WithinData
## Time1 Time2 Time3
## 1 5 7 8
## 2 6 7 8
## 3 6 8 9
## 4 7 8 9
## 5 8 9 9
Analyses of Multiple Variables
Confidence Intervals for the Means
ciMeans(Time1,Time2,Time3)
## N M SD SE LL UL
## Time1 5 6.4 1.140 0.510 4.984 7.816
## Time2 5 7.8 0.837 0.374 6.761 8.839
## Time3 5 8.6 0.548 0.245 7.920 9.280
ciMeans(Time1,Time2,Time3,conf.level=.99)
## N M SD SE LL UL
## Time1 5 6.4 1.140 0.510 4.052 8.748
## Time2 5 7.8 0.837 0.374 6.077 9.523
## Time3 5 8.6 0.548 0.245 7.472 9.728
Plots of Confidence Intervals for the Means
cipMeans(Time1,Time2,Time3)
cipMeans(Time1,Time2,Time3,conf.level=.99,mu=6)
Significance Tests for the Means
nhstMeans(Time1,Time2,Time3)
## Diff SE t df p
## Time1 6.4 0.510 12.551 4 0
## Time2 7.8 0.374 20.846 4 0
## Time3 8.6 0.245 35.109 4 0
nhstMeans(Time1,Time2,Time3,mu=6)
## Diff SE t df p
## Time1 0.4 0.510 0.784 4 0.477
## Time2 1.8 0.374 4.811 4 0.009
## Time3 2.6 0.245 10.614 4 0.000
Analyses of a Variable Comparison
Confidence Interval for the Mean Difference
ciDifference(Time1,Time2)
## Diff SE df LL UL
## -1.400 0.245 4.000 -2.080 -0.720
ciDifference(Time1,Time2,conf.level=.99)
## Diff SE df LL UL
## -1.400 0.245 4.000 -2.080 -0.720
ciDifference(Time3,Time1)
## Diff SE df LL UL
## 2.200 0.374 4.000 1.161 3.239
Plots of Confidence Intervals for the Mean Difference
cipDifference(Time1,Time2)
cipDifference(Time1,Time2,conf.level=.99)
Significance Test for the Mean Difference
nhstDifference(Time1,Time2)
## Diff SE t df p
## -1.400 0.245 -5.715 4.000 0.005
nhstDifference(Time1,Time2,mu=-2)
## Diff SE t df p
## -1.400 0.245 2.449 4.000 0.070
Analyses of Pairwise Comparisons
Confidence Intervals for the Pairwise Comparisons
ciPairwise(Time1,Time2,Time3)
## Diff SE df LL UL
## Time1 v Time2 -1.4 0.245 4 -2.080 -0.720
## Time1 v Time3 -2.2 0.374 4 -3.239 -1.161
## Time2 v Time3 -0.8 0.200 4 -1.355 -0.245
ciPairwise(Time1,Time2,Time3,conf.level=.99)
## Diff SE df LL UL
## Time1 v Time2 -1.4 0.245 4 -2.528 -0.272
## Time1 v Time3 -2.2 0.374 4 -3.923 -0.477
## Time2 v Time3 -0.8 0.200 4 -1.721 0.121
Plot of the Confidence Intervals for the Pairwise Comparisons
cipPairwise(Time1,Time2,Time3)
cipPairwise(Time1,Time2,Time3,mu=-2,conf.level=.99)
Significance Tests of the Pairwise Comparisons
nhstPairwise(Time1,Time2,Time3)
## Diff SE t df p
## Time1 v Time2 -1.4 0.245 -5.715 4 0.005
## Time1 v Time3 -2.2 0.374 -5.880 4 0.004
## Time2 v Time3 -0.8 0.200 -4.000 4 0.016
nhstPairwise(Time1,Time2,Time3,mu=-2)
## Diff SE t df p
## Time1 v Time2 -1.4 0.245 2.449 4 0.070
## Time1 v Time3 -2.2 0.374 -0.535 4 0.621
## Time2 v Time3 -0.8 0.200 6.000 4 0.004
Analyses of a Set of Contrasts
Confidence Intervals for the Set of Contrasts
ciContrasts(Time1,Time2,Time3,contrasts=contr.sum)
## Est SE LL UL
## (Intercept) 7.6 0.115 7.334 7.866
## Variable1 -1.2 0.163 -1.577 -0.823
## Variable2 0.2 0.163 -0.177 0.577
ciContrasts(Time1,Time2,Time3,contrasts=contr.sum,conf.level=.99)
## Est SE LL UL
## (Intercept) 7.6 0.115 7.213 7.987
## Variable1 -1.2 0.163 -1.748 -0.652
## Variable2 0.2 0.163 -0.348 0.748
ciContrasts(Time1,Time2,Time3,contrasts=contr.treatment)
## Est SE LL UL
## (Intercept) 6.4 0.200 5.939 6.861
## Variable2 1.4 0.283 0.748 2.052
## Variable3 2.2 0.283 1.548 2.852
ciContrasts(Time1,Time2,Time3,contrasts=contr.poly)
## Est SE LL UL
## (Intercept) 7.600 0.115 7.334 7.866
## Variable.L 1.556 0.200 1.094 2.017
## Variable.Q -0.245 0.200 -0.706 0.216
ciContrasts(Time1,Time2,Time3,contrasts=contr.helmert)
## Est SE LL UL
## (Intercept) 7.6 0.115 7.334 7.866
## Variable1 0.7 0.141 0.374 1.026
## Variable2 0.5 0.082 0.312 0.688
ciContrasts(Time1,Time2,Time3,contrasts=contr.SAS)
## Est SE LL UL
## (Intercept) 8.6 0.200 8.139 9.061
## Variable1 -2.2 0.283 -2.852 -1.548
## Variable2 -0.8 0.283 -1.452 -0.148
Plot of the Confidence Intervals for the Set of Contrasts
cipContrasts(Time1,Time2,Time3,contrasts=contr.sum)
cipContrasts(Time1,Time2,Time3,contrasts=contr.sum,conf.level=.99)
cipContrasts(Time1,Time2,Time3,contrasts=contr.treatment)
cipContrasts(Time1,Time2,Time3,contrasts=contr.poly)
cipContrasts(Time1,Time2,Time3,contrasts=contr.helmert)
cipContrasts(Time1,Time2,Time3,contrasts=contr.SAS)
Significance Tests of the Set of Contrasts
nhstContrasts(Time1,Time2,Time3,contrasts=contr.sum)
## Diff SE t p
## (Intercept) 7.6 0.115 65.818 0.000
## Variable1 -1.2 0.163 -7.348 0.000
## Variable2 0.2 0.163 1.225 0.256
nhstContrasts(Time1,Time2,Time3,contrasts=contr.treatment)
## Diff SE t p
## (Intercept) 6.4 0.200 32.000 0.000
## Variable2 1.4 0.283 4.950 0.001
## Variable3 2.2 0.283 7.778 0.000
nhstContrasts(Time1,Time2,Time3,contrasts=contr.poly)
## Diff SE t p
## (Intercept) 7.600 0.115 65.818 0.000
## Variable.L 1.556 0.200 7.778 0.000
## Variable.Q -0.245 0.200 -1.225 0.256
nhstContrasts(Time1,Time2,Time3,contrasts=contr.helmert)
## Diff SE t p
## (Intercept) 7.6 0.115 65.818 0.000
## Variable1 0.7 0.141 4.950 0.001
## Variable2 0.5 0.082 6.124 0.000
nhstContrasts(Time1,Time2,Time3,contrasts=contr.SAS)
## Diff SE t p
## (Intercept) 8.6 0.200 43.000 0.000
## Variable1 -2.2 0.283 -7.778 0.000
## Variable2 -0.8 0.283 -2.828 0.022
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.