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