Set Working Directory and Load in Data
setwd("H:/CHAYA")
mot <- read.csv("Longitudinal dataset_W1 to W6 (n=366).csv", fileEncoding="UTF-8-BOM")
mjmot <- mot[, grep(pattern = "MJmotive|^unique_id", colnames(mot))]
Load Libraries for Analysis
library(psych)
library(ggplot2)
library(corrplot) #plotting correlation matrices
library(GPArotation) #methods for factor rotation
library(nFactors) #methods for determining the number of factors
Look at Descriptive Statistics for Each Wave
#Wave 1
describe(mjmot[,c(87:103)])
## vars n mean sd median trimmed mad min max
## w1_MJmotive_enjoyment 1 364 3.71 1.07 4.00 3.81 1.24 1 5
## w1_MJmotive_conformity 2 364 1.33 0.64 1.00 1.17 0.00 1 5
## w1_MJmotive_coping 3 365 2.26 1.14 2.00 2.13 0.99 1 5
## w1_MJmotive_experiment 4 363 2.06 1.10 1.67 1.90 0.99 1 5
## w1_MJmotive_boredom 5 363 2.42 1.20 2.33 2.31 1.48 1 5
## w1_MJmotive_alcohol 6 365 1.71 0.87 1.33 1.57 0.49 1 5
## w1_MJmotive_celebration 7 363 3.06 1.17 3.00 3.04 1.48 1 5
## w1_MJmotive_alteredpercept 8 364 3.05 1.25 3.00 3.06 1.48 1 5
## w1_MJmotive_socialanxiety 9 364 2.74 1.28 2.67 2.68 1.48 1 5
## w1_MJmotive_relativelowrisk 10 362 3.14 1.26 3.00 3.17 1.48 1 5
## w1_MJmotive_sleep 11 362 3.38 1.30 3.50 3.46 1.73 1 5
## w1_MJmotive_availability 12 363 2.79 1.14 2.67 2.74 1.48 1 5
## w1_MJmotive_pain 13 363 3.03 1.34 3.00 3.03 1.48 1 5
## w1_MJmotive_nausea 14 365 2.18 1.29 1.67 1.99 0.99 1 5
## w1_MJmotive_attention 15 365 2.52 1.30 2.33 2.41 1.48 1 5
## w1_MJmotive_substitution 16 364 2.25 1.10 2.00 2.11 0.99 1 5
## w1_MJmotive_naturalmeds 17 363 2.74 1.36 2.67 2.67 1.48 1 5
## range skew kurtosis se
## w1_MJmotive_enjoyment 4 -0.60 -0.40 0.06
## w1_MJmotive_conformity 4 2.69 8.24 0.03
## w1_MJmotive_coping 4 0.84 -0.28 0.06
## w1_MJmotive_experiment 4 1.01 0.20 0.06
## w1_MJmotive_boredom 4 0.57 -0.75 0.06
## w1_MJmotive_alcohol 4 1.40 1.70 0.05
## w1_MJmotive_celebration 4 0.18 -1.00 0.06
## w1_MJmotive_alteredpercept 4 -0.06 -1.10 0.07
## w1_MJmotive_socialanxiety 4 0.32 -1.09 0.07
## w1_MJmotive_relativelowrisk 4 -0.10 -1.14 0.07
## w1_MJmotive_sleep 4 -0.30 -1.15 0.07
## w1_MJmotive_availability 4 0.28 -0.86 0.06
## w1_MJmotive_pain 4 0.02 -1.27 0.07
## w1_MJmotive_nausea 4 0.91 -0.45 0.07
## w1_MJmotive_attention 4 0.44 -1.05 0.07
## w1_MJmotive_substitution 4 0.91 0.23 0.06
## w1_MJmotive_naturalmeds 4 0.26 -1.22 0.07
#Wave 2
describe(mjmot[,c(2:18)])
## vars n mean sd median trimmed mad min max
## W2_MJmotive_enjoyment 1 339 3.33 1.16 3.33 3.39 1.48 1 5
## W2_MJmotive_conformity 2 339 1.23 0.55 1.00 1.09 0.00 1 4
## W2_MJmotive_coping 3 339 1.98 0.92 1.67 1.85 0.99 1 5
## W2_MJmotive_experiment 4 337 1.66 0.83 1.33 1.51 0.49 1 5
## W2_MJmotive_boredom 5 339 2.15 1.05 2.00 2.02 0.99 1 5
## W2_MJmotive_alcohol 6 339 1.56 0.79 1.33 1.39 0.49 1 5
## W2_MJmotive_celebration 7 339 2.72 1.15 2.33 2.66 0.99 1 5
## W2_MJmotive_alteredpercept 8 338 2.64 1.19 2.42 2.57 1.36 1 5
## W2_MJmotive_socialanxiety 9 338 2.29 1.10 2.00 2.18 1.48 1 5
## W2_MJmotive_relativelowrisk 10 338 2.72 1.31 2.67 2.66 1.48 1 5
## W2_MJmotive_sleep 11 336 3.25 1.22 3.33 3.28 1.48 1 5
## W2_MJmotive_availability 12 339 2.51 1.08 2.33 2.45 0.99 1 5
## W2_MJmotive_pain 13 339 2.75 1.25 2.67 2.70 1.48 1 5
## W2_MJmotive_nausea 14 339 1.97 1.14 1.67 1.79 0.99 1 5
## W2_MJmotive_attention 15 339 2.12 1.08 2.00 1.99 1.48 1 5
## W2_MJmotive_substitution 16 339 1.92 0.95 1.67 1.78 0.99 1 5
## W2_MJmotive_naturalmeds 17 338 2.50 1.30 2.33 2.39 1.48 1 5
## range skew kurtosis se
## W2_MJmotive_enjoyment 4 -0.36 -0.82 0.06
## W2_MJmotive_conformity 3 3.10 9.97 0.03
## W2_MJmotive_coping 4 1.02 0.44 0.05
## W2_MJmotive_experiment 4 1.54 2.30 0.05
## W2_MJmotive_boredom 4 0.83 -0.18 0.06
## W2_MJmotive_alcohol 4 1.78 3.10 0.04
## W2_MJmotive_celebration 4 0.48 -0.77 0.06
## W2_MJmotive_alteredpercept 4 0.39 -0.86 0.06
## W2_MJmotive_socialanxiety 4 0.65 -0.41 0.06
## W2_MJmotive_relativelowrisk 4 0.26 -1.17 0.07
## W2_MJmotive_sleep 4 -0.09 -1.19 0.07
## W2_MJmotive_availability 4 0.48 -0.71 0.06
## W2_MJmotive_pain 4 0.26 -1.12 0.07
## W2_MJmotive_nausea 4 1.13 0.21 0.06
## W2_MJmotive_attention 4 0.76 -0.32 0.06
## W2_MJmotive_substitution 4 1.09 0.81 0.05
## W2_MJmotive_naturalmeds 4 0.52 -0.97 0.07
#Wave 3
describe(mjmot[,c(19:35)])
## vars n mean sd median trimmed mad min max
## w3_MJmotive_enjoyment 1 322 3.29 1.18 3.33 3.34 1.48 1 5
## w3_MJmotive_conformity 2 322 1.21 0.51 1.00 1.09 0.00 1 5
## w3_MJmotive_coping 3 322 1.96 0.91 2.00 1.85 0.99 1 5
## w3_MJmotive_experiment 4 322 1.67 0.91 1.33 1.49 0.49 1 5
## w3_MJmotive_boredom 5 322 2.09 1.03 2.00 1.95 0.99 1 5
## w3_MJmotive_alcohol 6 322 1.50 0.74 1.33 1.35 0.49 1 5
## w3_MJmotive_celebration 7 322 2.61 1.16 2.33 2.53 0.99 1 5
## w3_MJmotive_alteredpercept 8 322 2.56 1.16 2.33 2.48 1.48 1 5
## w3_MJmotive_socialanxiety 9 321 2.26 1.13 2.00 2.13 1.48 1 5
## w3_MJmotive_relativelowrisk 10 322 2.73 1.33 2.67 2.67 1.98 1 5
## w3_MJmotive_sleep 11 322 3.11 1.25 3.00 3.13 1.48 1 5
## w3_MJmotive_availability 12 322 2.51 1.13 2.33 2.44 1.48 1 5
## w3_MJmotive_pain 13 322 2.67 1.26 2.33 2.60 1.48 1 5
## w3_MJmotive_nausea 14 322 1.89 1.09 1.50 1.69 0.74 1 5
## w3_MJmotive_attention 15 322 2.05 1.07 1.67 1.90 0.99 1 5
## w3_MJmotive_substitution 16 322 1.92 0.99 1.67 1.75 0.99 1 5
## w3_MJmotive_naturalmeds 17 322 2.48 1.26 2.33 2.36 1.48 1 5
## range skew kurtosis se
## w3_MJmotive_enjoyment 4 -0.28 -0.95 0.07
## w3_MJmotive_conformity 4 3.74 17.58 0.03
## w3_MJmotive_coping 4 0.94 0.33 0.05
## w3_MJmotive_experiment 4 1.69 2.66 0.05
## w3_MJmotive_boredom 4 1.02 0.44 0.06
## w3_MJmotive_alcohol 4 2.21 5.97 0.04
## w3_MJmotive_celebration 4 0.55 -0.63 0.06
## w3_MJmotive_alteredpercept 4 0.37 -0.86 0.06
## w3_MJmotive_socialanxiety 4 0.73 -0.38 0.06
## w3_MJmotive_relativelowrisk 4 0.19 -1.26 0.07
## w3_MJmotive_sleep 4 -0.06 -1.17 0.07
## w3_MJmotive_availability 4 0.42 -0.80 0.06
## w3_MJmotive_pain 4 0.36 -1.09 0.07
## w3_MJmotive_nausea 4 1.38 1.11 0.06
## w3_MJmotive_attention 4 0.96 0.10 0.06
## w3_MJmotive_substitution 4 1.24 1.20 0.06
## w3_MJmotive_naturalmeds 4 0.53 -0.85 0.07
#Wave 4
describe(mjmot[,c(36:52)])
## vars n mean sd median trimmed mad min max
## w4_MJmotive_enjoyment 1 302 3.18 1.09 3.33 3.21 0.99 1 5.00
## w4_MJmotive_conformity 2 302 1.32 0.64 1.00 1.16 0.00 1 4.67
## w4_MJmotive_coping 3 302 2.05 0.92 2.00 1.94 0.99 1 5.00
## w4_MJmotive_experiment 4 302 1.72 0.88 1.33 1.56 0.49 1 4.67
## w4_MJmotive_boredom 5 302 2.19 1.01 2.00 2.07 0.99 1 5.00
## w4_MJmotive_alcohol 6 302 1.58 0.78 1.33 1.43 0.49 1 4.00
## w4_MJmotive_celebration 7 302 2.57 1.11 2.33 2.50 0.99 1 5.00
## w4_MJmotive_alteredpercept 8 302 2.55 1.11 2.33 2.48 0.99 1 5.00
## w4_MJmotive_socialanxiety 9 302 2.23 1.12 2.00 2.09 1.48 1 5.00
## w4_MJmotive_relativelowrisk 10 302 2.56 1.24 2.33 2.47 1.48 1 5.00
## w4_MJmotive_sleep 11 302 2.98 1.20 3.00 2.97 1.48 1 5.00
## w4_MJmotive_availability 12 302 2.36 1.05 2.33 2.27 0.99 1 5.00
## w4_MJmotive_pain 13 302 2.48 1.13 2.33 2.40 0.99 1 5.00
## w4_MJmotive_nausea 14 302 1.94 1.07 1.67 1.75 0.99 1 5.00
## w4_MJmotive_attention 15 302 1.98 1.04 1.67 1.83 0.99 1 5.00
## w4_MJmotive_substitution 16 302 1.88 0.92 1.67 1.75 0.99 1 5.00
## w4_MJmotive_naturalmeds 17 302 2.33 1.22 2.00 2.19 1.48 1 5.00
## range skew kurtosis se
## w4_MJmotive_enjoyment 4.00 -0.21 -0.71 0.06
## w4_MJmotive_conformity 3.67 2.47 6.18 0.04
## w4_MJmotive_coping 4.00 0.88 0.34 0.05
## w4_MJmotive_experiment 3.67 1.36 1.33 0.05
## w4_MJmotive_boredom 4.00 0.86 0.15 0.06
## w4_MJmotive_alcohol 3.00 1.41 1.11 0.04
## w4_MJmotive_celebration 4.00 0.51 -0.65 0.06
## w4_MJmotive_alteredpercept 4.00 0.46 -0.62 0.06
## w4_MJmotive_socialanxiety 4.00 0.81 -0.19 0.06
## w4_MJmotive_relativelowrisk 4.00 0.42 -0.91 0.07
## w4_MJmotive_sleep 4.00 0.18 -1.05 0.07
## w4_MJmotive_availability 4.00 0.63 -0.24 0.06
## w4_MJmotive_pain 4.00 0.48 -0.74 0.06
## w4_MJmotive_nausea 4.00 1.28 0.90 0.06
## w4_MJmotive_attention 4.00 1.05 0.32 0.06
## w4_MJmotive_substitution 4.00 0.98 0.40 0.05
## w4_MJmotive_naturalmeds 4.00 0.71 -0.56 0.07
#Wave 5
describe(mjmot[,c(53:69)])
## vars n mean sd median trimmed mad min max
## w5_MJmotive_enjoyment 1 260 3.25 1.12 3.33 3.28 1.48 1 5
## w5_MJmotive_conformity 2 260 1.38 0.67 1.00 1.23 0.00 1 4
## w5_MJmotive_coping 3 260 2.21 0.99 2.00 2.10 0.99 1 5
## w5_MJmotive_experiment 4 260 1.86 1.04 1.33 1.67 0.49 1 5
## w5_MJmotive_boredom 5 260 2.27 1.08 2.00 2.16 0.99 1 5
## w5_MJmotive_alcohol 6 260 1.71 0.84 1.33 1.57 0.49 1 5
## w5_MJmotive_celebration 7 260 2.57 1.13 2.33 2.50 0.99 1 5
## w5_MJmotive_alteredpercept 8 260 2.67 1.22 2.33 2.60 1.48 1 5
## w5_MJmotive_socialanxiety 9 260 2.33 1.14 2.00 2.21 0.99 1 5
## w5_MJmotive_relativelowrisk 10 260 2.64 1.21 2.33 2.58 1.48 1 5
## w5_MJmotive_sleep 11 260 2.98 1.21 3.00 2.96 1.48 1 5
## w5_MJmotive_availability 12 260 2.40 1.07 2.33 2.33 0.99 1 5
## w5_MJmotive_pain 13 260 2.67 1.19 2.33 2.60 1.48 1 5
## w5_MJmotive_nausea 14 260 2.07 1.10 2.00 1.90 1.48 1 5
## w5_MJmotive_attention 15 260 2.14 1.16 2.00 1.98 1.48 1 5
## w5_MJmotive_substitution 16 260 1.96 0.98 1.67 1.81 0.99 1 5
## w5_MJmotive_naturalmeds 17 260 2.43 1.28 2.17 2.29 1.24 1 5
## range skew kurtosis se
## w5_MJmotive_enjoyment 4 -0.21 -0.94 0.07
## w5_MJmotive_conformity 3 1.89 2.91 0.04
## w5_MJmotive_coping 4 0.78 0.04 0.06
## w5_MJmotive_experiment 4 1.32 1.05 0.06
## w5_MJmotive_boredom 4 0.71 -0.31 0.07
## w5_MJmotive_alcohol 4 1.14 0.45 0.05
## w5_MJmotive_celebration 4 0.50 -0.68 0.07
## w5_MJmotive_alteredpercept 4 0.37 -0.96 0.08
## w5_MJmotive_socialanxiety 4 0.71 -0.46 0.07
## w5_MJmotive_relativelowrisk 4 0.31 -1.03 0.08
## w5_MJmotive_sleep 4 0.16 -1.11 0.08
## w5_MJmotive_availability 4 0.41 -0.84 0.07
## w5_MJmotive_pain 4 0.38 -0.91 0.07
## w5_MJmotive_nausea 4 1.00 0.25 0.07
## w5_MJmotive_attention 4 0.89 -0.20 0.07
## w5_MJmotive_substitution 4 1.02 0.56 0.06
## w5_MJmotive_naturalmeds 4 0.63 -0.77 0.08
#Wave 6
describe(mjmot[,c(70:86)])
## vars n mean sd median trimmed mad min max
## w6_MJmotive_enjoyment 1 251 3.25 1.17 3.33 3.30 0.99 1 5
## w6_MJmotive_conformity 2 251 1.38 0.76 1.00 1.18 0.00 1 5
## w6_MJmotive_coping 3 251 2.13 1.06 2.00 2.00 0.99 1 5
## w6_MJmotive_experiment 4 251 1.76 1.01 1.33 1.57 0.49 1 5
## w6_MJmotive_boredom 5 251 2.31 1.14 2.00 2.19 0.99 1 5
## w6_MJmotive_alcohol 6 251 1.62 0.86 1.33 1.46 0.49 1 5
## w6_MJmotive_celebration 7 251 2.54 1.15 2.33 2.47 0.99 1 5
## w6_MJmotive_alteredpercept 8 251 2.53 1.21 2.33 2.45 1.48 1 5
## w6_MJmotive_socialanxiety 9 251 2.23 1.16 2.00 2.09 1.48 1 5
## w6_MJmotive_relativelowrisk 10 251 2.49 1.28 2.33 2.39 1.48 1 5
## w6_MJmotive_sleep 11 251 2.83 1.22 2.67 2.79 1.48 1 5
## w6_MJmotive_availability 12 251 2.29 1.07 2.00 2.19 0.99 1 5
## w6_MJmotive_pain 13 251 2.43 1.23 2.00 2.32 1.48 1 5
## w6_MJmotive_nausea 14 251 1.91 1.12 1.33 1.71 0.49 1 5
## w6_MJmotive_attention 15 251 1.97 1.09 1.67 1.81 0.99 1 5
## w6_MJmotive_substitution 16 251 1.88 1.03 1.67 1.72 0.99 1 5
## w6_MJmotive_naturalmeds 17 251 2.25 1.24 2.00 2.09 1.48 1 5
## range skew kurtosis se
## w6_MJmotive_enjoyment 4 -0.37 -0.78 0.07
## w6_MJmotive_conformity 4 2.43 6.14 0.05
## w6_MJmotive_coping 4 0.92 0.14 0.07
## w6_MJmotive_experiment 4 1.32 0.84 0.06
## w6_MJmotive_boredom 4 0.71 -0.44 0.07
## w6_MJmotive_alcohol 4 1.61 2.52 0.05
## w6_MJmotive_celebration 4 0.49 -0.68 0.07
## w6_MJmotive_alteredpercept 4 0.35 -1.00 0.08
## w6_MJmotive_socialanxiety 4 0.72 -0.56 0.07
## w6_MJmotive_relativelowrisk 4 0.39 -1.12 0.08
## w6_MJmotive_sleep 4 0.21 -0.98 0.08
## w6_MJmotive_availability 4 0.68 -0.35 0.07
## w6_MJmotive_pain 4 0.56 -0.81 0.08
## w6_MJmotive_nausea 4 1.20 0.50 0.07
## w6_MJmotive_attention 4 0.99 0.03 0.07
## w6_MJmotive_substitution 4 1.13 0.71 0.06
## w6_MJmotive_naturalmeds 4 0.78 -0.53 0.08
Visual Correlation Matrix
corrplot(cor(mjmot[,c(87:103)], use="complete.obs"), order = "original", tl.col='black', tl.cex=.75)
corrplot(cor(mjmot[,c(2:18)], use="complete.obs"), order = "original", tl.col='black', tl.cex=.75)
corrplot(cor(mjmot[,c(19:35)], use="complete.obs"), order = "original", tl.col='black', tl.cex=.75)
corrplot(cor(mjmot[,c(36:52)], use="complete.obs"), order = "original", tl.col='black', tl.cex=.75)
corrplot(cor(mjmot[,c(53:69)], use="complete.obs"), order = "original", tl.col='black', tl.cex=.75)
corrplot(cor(mjmot[,c(70:86)], use="complete.obs"), order = "original", tl.col='black', tl.cex=.75)
Run ORTHOGONAL EFAs
EFAresult1 = factanal(~ ., data=mjmot[,c(87:103)], factors = 3, rotation = "varimax",
na.action = na.exclude, cut = .3) #note the formula specification allows NA
EFAresult2 = factanal(~ ., data=mjmot[,c(2:18)], factors = 3, rotation = "varimax",
na.action = na.exclude) #note the formula specification allows NA
EFAresult3 = factanal(~ ., data=mjmot[,c(19:35)], factors = 3, rotation = "varimax",
na.action = na.exclude) #note the formula specification allows NA
EFAresult4 = factanal(~ ., data=mjmot[,c(36:52)], factors = 3, rotation = "varimax",
na.action = na.exclude) #note the formula specification allows NA
EFAresult5 = factanal(~ ., data=mjmot[,c(53:69)], factors = 3, rotation = "varimax",
na.action = na.exclude) #note the formula specification allows NA
EFAresult6 = factanal(~ ., data=mjmot[,c(70:86)], factors = 3, rotation = "varimax",
na.action = na.exclude) #note the formula specification allows NA
View EFA Results
EFAresult1
##
## Call:
## factanal(x = ~., factors = 3, data = mjmot[, c(87:103)], na.action = na.exclude, rotation = "varimax", cut = 0.3)
##
## Uniquenesses:
## w1_MJmotive_enjoyment w1_MJmotive_conformity
## 0.401 0.646
## w1_MJmotive_coping w1_MJmotive_experiment
## 0.618 0.624
## w1_MJmotive_boredom w1_MJmotive_alcohol
## 0.308 0.728
## w1_MJmotive_celebration w1_MJmotive_alteredpercept
## 0.577 0.465
## w1_MJmotive_socialanxiety w1_MJmotive_relativelowrisk
## 0.363 0.525
## w1_MJmotive_sleep w1_MJmotive_availability
## 0.530 0.435
## w1_MJmotive_pain w1_MJmotive_nausea
## 0.350 0.403
## w1_MJmotive_attention w1_MJmotive_substitution
## 0.463 0.422
## w1_MJmotive_naturalmeds
## 0.414
##
## Loadings:
## Factor1 Factor2 Factor3
## w1_MJmotive_enjoyment 0.127 0.354 0.676
## w1_MJmotive_conformity 0.589
## w1_MJmotive_coping 0.364 0.473 0.159
## w1_MJmotive_experiment 0.234 0.439 0.358
## w1_MJmotive_boredom 0.725 0.401
## w1_MJmotive_alcohol 0.163 0.488
## w1_MJmotive_celebration 0.288 0.414 0.410
## w1_MJmotive_alteredpercept 0.322 0.197 0.627
## w1_MJmotive_socialanxiety 0.661 0.257 0.366
## w1_MJmotive_relativelowrisk 0.377 0.184 0.547
## w1_MJmotive_sleep 0.669 0.106 0.102
## w1_MJmotive_availability 0.579 0.479
## w1_MJmotive_pain 0.782 0.158 0.115
## w1_MJmotive_nausea 0.754 0.163
## w1_MJmotive_attention 0.687 0.243
## w1_MJmotive_substitution 0.609 0.380 0.250
## w1_MJmotive_naturalmeds 0.729 0.227
##
## Factor1 Factor2 Factor3
## SS loadings 4.003 2.515 2.210
## Proportion Var 0.235 0.148 0.130
## Cumulative Var 0.235 0.383 0.513
##
## Test of the hypothesis that 3 factors are sufficient.
## The chi square statistic is 252.85 on 88 degrees of freedom.
## The p-value is 7.41e-18
EFAresult2
##
## Call:
## factanal(x = ~., factors = 3, data = mjmot[, c(2:18)], na.action = na.exclude, rotation = "varimax")
##
## Uniquenesses:
## W2_MJmotive_enjoyment W2_MJmotive_conformity
## 0.272 0.534
## W2_MJmotive_coping W2_MJmotive_experiment
## 0.703 0.561
## W2_MJmotive_boredom W2_MJmotive_alcohol
## 0.443 0.557
## W2_MJmotive_celebration W2_MJmotive_alteredpercept
## 0.593 0.569
## W2_MJmotive_socialanxiety W2_MJmotive_relativelowrisk
## 0.461 0.511
## W2_MJmotive_sleep W2_MJmotive_availability
## 0.618 0.444
## W2_MJmotive_pain W2_MJmotive_nausea
## 0.314 0.394
## W2_MJmotive_attention W2_MJmotive_substitution
## 0.467 0.538
## W2_MJmotive_naturalmeds
## 0.406
##
## Loadings:
## Factor1 Factor2 Factor3
## W2_MJmotive_enjoyment 0.131 0.824 0.179
## W2_MJmotive_conformity 0.680
## W2_MJmotive_coping 0.365 0.255 0.314
## W2_MJmotive_experiment 0.286 0.308 0.512
## W2_MJmotive_boredom 0.495 0.556
## W2_MJmotive_alcohol 0.128 0.109 0.644
## W2_MJmotive_celebration 0.252 0.461 0.362
## W2_MJmotive_alteredpercept 0.367 0.508 0.193
## W2_MJmotive_socialanxiety 0.575 0.373 0.263
## W2_MJmotive_relativelowrisk 0.394 0.567 0.112
## W2_MJmotive_sleep 0.597 0.151
## W2_MJmotive_availability 0.475 0.574
## W2_MJmotive_pain 0.826
## W2_MJmotive_nausea 0.758 0.174
## W2_MJmotive_attention 0.705 0.145 0.123
## W2_MJmotive_substitution 0.545 0.353 0.200
## W2_MJmotive_naturalmeds 0.741 0.185 -0.103
##
## Factor1 Factor2 Factor3
## SS loadings 3.891 2.462 2.260
## Proportion Var 0.229 0.145 0.133
## Cumulative Var 0.229 0.374 0.507
##
## Test of the hypothesis that 3 factors are sufficient.
## The chi square statistic is 248.81 on 88 degrees of freedom.
## The p-value is 2.82e-17
EFAresult3
##
## Call:
## factanal(x = ~., factors = 3, data = mjmot[, c(19:35)], na.action = na.exclude, rotation = "varimax")
##
## Uniquenesses:
## w3_MJmotive_enjoyment w3_MJmotive_conformity
## 0.335 0.403
## w3_MJmotive_coping w3_MJmotive_experiment
## 0.611 0.497
## w3_MJmotive_boredom w3_MJmotive_alcohol
## 0.445 0.468
## w3_MJmotive_celebration w3_MJmotive_alteredpercept
## 0.468 0.446
## w3_MJmotive_socialanxiety w3_MJmotive_relativelowrisk
## 0.423 0.496
## w3_MJmotive_sleep w3_MJmotive_availability
## 0.530 0.408
## w3_MJmotive_pain w3_MJmotive_nausea
## 0.350 0.454
## w3_MJmotive_attention w3_MJmotive_substitution
## 0.483 0.485
## w3_MJmotive_naturalmeds
## 0.380
##
## Loadings:
## Factor1 Factor2 Factor3
## w3_MJmotive_enjoyment 0.136 0.800
## w3_MJmotive_conformity 0.767
## w3_MJmotive_coping 0.313 0.388 0.375
## w3_MJmotive_experiment 0.237 0.421 0.519
## w3_MJmotive_boredom 0.586 0.460
## w3_MJmotive_alcohol 0.199 0.696
## w3_MJmotive_celebration 0.308 0.613 0.248
## w3_MJmotive_alteredpercept 0.360 0.608 0.233
## w3_MJmotive_socialanxiety 0.567 0.456 0.219
## w3_MJmotive_relativelowrisk 0.373 0.603
## w3_MJmotive_sleep 0.678
## w3_MJmotive_availability 0.665 0.386
## w3_MJmotive_pain 0.794 0.123
## w3_MJmotive_nausea 0.721 0.146
## w3_MJmotive_attention 0.626 0.317 0.156
## w3_MJmotive_substitution 0.485 0.387 0.361
## w3_MJmotive_naturalmeds 0.763 0.192
##
## Factor1 Factor2 Factor3
## SS loadings 3.686 3.437 2.196
## Proportion Var 0.217 0.202 0.129
## Cumulative Var 0.217 0.419 0.548
##
## Test of the hypothesis that 3 factors are sufficient.
## The chi square statistic is 245.55 on 88 degrees of freedom.
## The p-value is 8.21e-17
EFAresult4
##
## Call:
## factanal(x = ~., factors = 3, data = mjmot[, c(36:52)], na.action = na.exclude, rotation = "varimax")
##
## Uniquenesses:
## w4_MJmotive_enjoyment w4_MJmotive_conformity
## 0.352 0.336
## w4_MJmotive_coping w4_MJmotive_experiment
## 0.561 0.428
## w4_MJmotive_boredom w4_MJmotive_alcohol
## 0.346 0.561
## w4_MJmotive_celebration w4_MJmotive_alteredpercept
## 0.481 0.430
## w4_MJmotive_socialanxiety w4_MJmotive_relativelowrisk
## 0.327 0.498
## w4_MJmotive_sleep w4_MJmotive_availability
## 0.526 0.405
## w4_MJmotive_pain w4_MJmotive_nausea
## 0.316 0.461
## w4_MJmotive_attention w4_MJmotive_substitution
## 0.415 0.505
## w4_MJmotive_naturalmeds
## 0.358
##
## Loadings:
## Factor1 Factor2 Factor3
## w4_MJmotive_enjoyment 0.150 0.790
## w4_MJmotive_conformity 0.809
## w4_MJmotive_coping 0.371 0.349 0.424
## w4_MJmotive_experiment 0.262 0.324 0.631
## w4_MJmotive_boredom 0.676 0.438
## w4_MJmotive_alcohol 0.119 0.237 0.607
## w4_MJmotive_celebration 0.310 0.603 0.241
## w4_MJmotive_alteredpercept 0.351 0.634 0.211
## w4_MJmotive_socialanxiety 0.646 0.470 0.186
## w4_MJmotive_relativelowrisk 0.418 0.537 0.199
## w4_MJmotive_sleep 0.656 0.209
## w4_MJmotive_availability 0.663 0.391
## w4_MJmotive_pain 0.820
## w4_MJmotive_nausea 0.710 0.158
## w4_MJmotive_attention 0.657 0.362 0.150
## w4_MJmotive_substitution 0.540 0.293 0.343
## w4_MJmotive_naturalmeds 0.791 0.104
##
## Factor1 Factor2 Factor3
## SS loadings 4.018 3.370 2.305
## Proportion Var 0.236 0.198 0.136
## Cumulative Var 0.236 0.435 0.570
##
## Test of the hypothesis that 3 factors are sufficient.
## The chi square statistic is 375.37 on 88 degrees of freedom.
## The p-value is 3.79e-37
EFAresult5
##
## Call:
## factanal(x = ~., factors = 3, data = mjmot[, c(53:69)], na.action = na.exclude, rotation = "varimax")
##
## Uniquenesses:
## w5_MJmotive_enjoyment w5_MJmotive_conformity
## 0.304 0.342
## w5_MJmotive_coping w5_MJmotive_experiment
## 0.604 0.477
## w5_MJmotive_boredom w5_MJmotive_alcohol
## 0.354 0.547
## w5_MJmotive_celebration w5_MJmotive_alteredpercept
## 0.377 0.339
## w5_MJmotive_socialanxiety w5_MJmotive_relativelowrisk
## 0.350 0.498
## w5_MJmotive_sleep w5_MJmotive_availability
## 0.478 0.390
## w5_MJmotive_pain w5_MJmotive_nausea
## 0.251 0.379
## w5_MJmotive_attention w5_MJmotive_substitution
## 0.323 0.517
## w5_MJmotive_naturalmeds
## 0.352
##
## Loadings:
## Factor1 Factor2 Factor3
## w5_MJmotive_enjoyment 0.182 0.814
## w5_MJmotive_conformity 0.110 0.102 0.797
## w5_MJmotive_coping 0.344 0.407 0.335
## w5_MJmotive_experiment 0.302 0.477 0.452
## w5_MJmotive_boredom 0.647 0.477
## w5_MJmotive_alcohol 0.137 0.186 0.632
## w5_MJmotive_celebration 0.336 0.685 0.203
## w5_MJmotive_alteredpercept 0.374 0.705 0.156
## w5_MJmotive_socialanxiety 0.624 0.469 0.201
## w5_MJmotive_relativelowrisk 0.442 0.540 0.122
## w5_MJmotive_sleep 0.673 0.261
## w5_MJmotive_availability 0.663 0.407
## w5_MJmotive_pain 0.856 0.117
## w5_MJmotive_nausea 0.735 0.148 0.244
## w5_MJmotive_attention 0.729 0.311 0.220
## w5_MJmotive_substitution 0.503 0.307 0.368
## w5_MJmotive_naturalmeds 0.794 0.108
##
## Factor1 Factor2 Factor3
## SS loadings 4.254 3.739 2.126
## Proportion Var 0.250 0.220 0.125
## Cumulative Var 0.250 0.470 0.595
##
## Test of the hypothesis that 3 factors are sufficient.
## The chi square statistic is 258.45 on 88 degrees of freedom.
## The p-value is 1.14e-18
EFAresult6
##
## Call:
## factanal(x = ~., factors = 3, data = mjmot[, c(70:86)], na.action = na.exclude, rotation = "varimax")
##
## Uniquenesses:
## w6_MJmotive_enjoyment w6_MJmotive_conformity
## 0.290 0.278
## w6_MJmotive_coping w6_MJmotive_experiment
## 0.510 0.385
## w6_MJmotive_boredom w6_MJmotive_alcohol
## 0.326 0.397
## w6_MJmotive_celebration w6_MJmotive_alteredpercept
## 0.491 0.329
## w6_MJmotive_socialanxiety w6_MJmotive_relativelowrisk
## 0.309 0.494
## w6_MJmotive_sleep w6_MJmotive_availability
## 0.501 0.368
## w6_MJmotive_pain w6_MJmotive_nausea
## 0.210 0.301
## w6_MJmotive_attention w6_MJmotive_substitution
## 0.313 0.434
## w6_MJmotive_naturalmeds
## 0.251
##
## Loadings:
## Factor1 Factor2 Factor3
## w6_MJmotive_enjoyment 0.137 0.830
## w6_MJmotive_conformity 0.191 0.823
## w6_MJmotive_coping 0.321 0.511 0.355
## w6_MJmotive_experiment 0.373 0.383 0.573
## w6_MJmotive_boredom 0.149 0.654 0.473
## w6_MJmotive_alcohol 0.249 0.729
## w6_MJmotive_celebration 0.266 0.616 0.243
## w6_MJmotive_alteredpercept 0.380 0.710 0.149
## w6_MJmotive_socialanxiety 0.637 0.485 0.223
## w6_MJmotive_relativelowrisk 0.434 0.531 0.189
## w6_MJmotive_sleep 0.642 0.287
## w6_MJmotive_availability 0.605 0.507
## w6_MJmotive_pain 0.874 0.147
## w6_MJmotive_nausea 0.786 0.114 0.263
## w6_MJmotive_attention 0.711 0.344 0.252
## w6_MJmotive_substitution 0.607 0.269 0.354
## w6_MJmotive_naturalmeds 0.852 0.138
##
## Factor1 Factor2 Factor3
## SS loadings 4.539 3.690 2.586
## Proportion Var 0.267 0.217 0.152
## Cumulative Var 0.267 0.484 0.636
##
## Test of the hypothesis that 3 factors are sufficient.
## The chi square statistic is 266.51 on 88 degrees of freedom.
## The p-value is 7.51e-20