# Set the working directory to the location on your computer to which you downloaded the files from OSF. You will also need to download the datasets to this exact same location - Do NOT rename the files!!!
x <- getwd()
setwd(x)
library(knitr)
library(reshape2)
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(lattice)
library(magrittr)
library(lme4)
## Warning: package 'lme4' was built under R version 3.6.2
## Loading required package: Matrix
library(lmerTest)
## Warning: package 'lmerTest' was built under R version 3.6.2
## 
## Attaching package: 'lmerTest'
## The following object is masked from 'package:lme4':
## 
##     lmer
## The following object is masked from 'package:stats':
## 
##     step
library(knitr)
library (optimx)
## Warning: package 'optimx' was built under R version 3.6.3
library(rmarkdown)
opts_chunk$set(fig.path='figure/',fig.ext=NULL)
data <- read.table("bfijudgement.w.snsfreq.csv" ,header = T, sep = ",")
colnames(data) [1]<- "JID"

IG.Cond <- read.table("IGCon.csv", header = T, sep = ",")
TW.Cond <- read.table("TWCon.csv", header = T, sep = ",")

#All of this is just for the analysis that looks at the number of original Twitter posts
TW.2 <- read.table ("TW2O.csv", header = T, sep = ",")
TW.4 <- read.table ("TW4O.csv", header = T, sep = ",")
TW.6 <- read.table ("TW6O.csv", header = T, sep = ",")
TW.7 <- read.table ("TW7O.csv", header = T, sep = ",")
TW.9 <- read.table ("TW9O.csv", header = T, sep = ",")
TW.11 <- read.table ("TW11O.csv", header = T, sep = ",")
TW.14 <- read.table ("TW14O.csv", header = T, sep = ",")
TW.16 <- read.table ("TW16O.csv", header = T, sep = ",")
TW.18 <- read.table ("TW18O.csv", header = T, sep = ",")
TW.20 <- read.table ("TW20O.csv", header = T, sep = ",")
TW.22 <- read.table ("TW22O.csv", header = T, sep = ",")
TW.24 <- read.table ("TW24O.csv", header = T, sep = ",")
TW.25 <- read.table ("TW25O.csv", header = T, sep = ",")
TW.27 <- read.table ("TW27O.csv", header = T, sep = ",")
TW.29 <- read.table ("TW29O.csv", header = T, sep = ",")
TW.32 <- read.table ("TW32O.csv", header = T, sep = ",")
TW.34 <- read.table ("TW34O.csv", header = T, sep = ",")
TW.36 <- read.table ("TW36O.csv", header = T, sep = ",")
TW.38 <- read.table ("TW38O.csv", header = T, sep = ",")
TW.40 <- read.table ("TW40O.csv", header = T, sep = ",")
TW.42 <- read.table ("TW42O.csv", header = T, sep = ",")
TW.43 <- read.table ("TW43O.csv", header = T, sep = ",")
TW.45 <- read.table ("TW45O.csv", header = T, sep = ",")
TW.47 <- read.table ("TW47O.csv", header = T, sep = ",")
TW.50 <- read.table ("TW50O.csv", header = T, sep = ",")
TW.52 <- read.table ("TW52O.csv", header = T, sep = ",")
TW.54 <- read.table ("TW54O.csv", header = T, sep = ",")
TW.56 <- read.table ("TW56O.csv", header = T, sep = ",")
TW.58 <- read.table ("TW58O.csv", header = T, sep = ",")
TW.60 <- read.table ("TW60O.csv", header = T, sep = ",")
TW.61 <- read.table ("TW61O.csv", header = T, sep = ",")
TW.63 <- read.table ("TW63O.csv", header = T, sep = ",")
TW.65 <- read.table ("TW65O.csv", header = T, sep = ",")
TW.68 <- read.table ("TW68O.csv", header = T, sep = ",")
TW.70 <- read.table ("TW70O.csv", header = T, sep = ",")
TW.72 <- read.table ("TW72O.csv", header = T, sep = ",")
TW.74 <- read.table ("TW74O.csv", header = T, sep = ",")
TW.76 <- read.table ("TW76O.csv", header = T, sep = ",")
TW.78 <- read.table ("TW78O.csv", header = T, sep = ",")
TW.79 <- read.table ("TW79O.csv", header = T, sep = ",")
TW.81 <- read.table ("TW81O.csv", header = T, sep = ",")
TW.83 <- read.table ("TW83O.csv", header = T, sep = ",")
TW.86 <- read.table ("TW86O.csv", header = T, sep = ",")
TW.88 <- read.table ("TW88O.csv", header = T, sep = ",")
TW.90 <- read.table ("TW90O.csv", header = T, sep = ",")
TW.92 <- read.table ("TW92O.csv", header = T, sep = ",")
TW.94 <- read.table ("TW94O.csv", header = T, sep = ",")
TW.96 <- read.table ("TW96O.csv", header = T, sep = ",")
TW.97 <- read.table ("TW97O.csv", header = T, sep = ",")
TW.99 <- read.table ("TW99O.csv", header = T, sep = ",")
TW.101 <- read.table ("TW101O.csv", header = T, sep = ",")
TW.104 <- read.table ("TW104O.csv", header = T, sep = ",")
TW.106 <- read.table ("TW106O.csv", header = T, sep = ",")
TW.108 <- read.table ("TW108O.csv", header = T, sep = ",")
TW.110 <- read.table ("TW110O.csv", header = T, sep = ",")
TW.112 <- read.table ("TW112O.csv", header = T, sep = ",")
TW.114 <- read.table ("TW114O.csv", header = T, sep = ",")
TW.115 <- read.table ("TW115O.csv", header = T, sep = ",")
TW.117 <- read.table ("TW117O.csv", header = T, sep = ",")
TW.119 <- read.table ("TW119O.csv", header = T, sep = ",")
TW.122 <- read.table ("TW122O.csv", header = T, sep = ",")
TW.124 <- read.table ("TW124O.csv", header = T, sep = ",")
TW.126 <- read.table ("TW126O.csv", header = T, sep = ",")
TW.128 <- read.table ("TW128O.csv", header = T, sep = ",")
TW.130 <- read.table ("TW130O.csv", header = T, sep = ",")
TW.132 <- read.table ("TW132O.csv", header = T, sep = ",")
TW.133 <- read.table ("TW133O.csv", header = T, sep = ",")
TW.135 <- read.table ("TW135O.csv", header = T, sep = ",")
TW.137 <- read.table ("TW137O.csv", header = T, sep = ",")
TW.140 <- read.table ("TW140O.csv", header = T, sep = ",")
TW.142 <- read.table ("TW122O.csv", header = T, sep = ",")
TW.144 <- read.table ("TW124O.csv", header = T, sep = ",")
TW.146 <- read.table ("TW126O.csv", header = T, sep = ",")
TW.148 <- read.table ("TW128O.csv", header = T, sep = ",")
TW.150 <- read.table ("TW130O.csv", header = T, sep = ",")

## Each target has a realistic accuracy criterion composed of a self-rating
sr.targets = read.table("targets1.csv", fill= TRUE, header = T, sep = ",")

colnames(sr.targets) [1]<- 'TID'

## In order to properly organize the data, there must be a seperate dataframe for each targets' validity ratings.
sr.1   <- sr.targets %>% filter(TID == "1")

sr.2   <- sr.targets %>% filter(TID == "2")

sr.3   <- sr.targets %>% filter(TID == "3")

sr.4   <- sr.targets %>% filter(TID == "4")

sr.5   <- sr.targets %>% filter(TID == "5")
sr.6   <- sr.targets %>% filter(TID == "6")

sr.7   <- sr.targets %>% filter(TID == "7")

sr.8   <- sr.targets %>% filter(TID == "8")

sr.9   <- sr.targets %>% filter(TID == "9")

sr.10   <- sr.targets %>% filter(TID == "10")
sr.11  <- sr.targets %>% filter(TID == "11")

sr.12   <- sr.targets %>% filter(TID == "12")

sr.13   <- sr.targets %>% filter(TID == "13")

sr.14   <- sr.targets %>% filter(TID == "14")

sr.15   <- sr.targets %>% filter(TID == "15")
sr.16   <- sr.targets %>% filter(TID == "16")

sr.17   <- sr.targets %>% filter(TID == "17")

sr.18   <- sr.targets %>% filter(TID == "18")

sr.19   <- sr.targets %>% filter(TID == "19")

sr.20   <- sr.targets %>% filter(TID == "20")
sr.21   <- sr.targets %>% filter(TID == "21")

sr.22   <- sr.targets %>% filter(TID == "22")

sr.23   <- sr.targets %>% filter(TID == "23")

sr.24   <- sr.targets %>% filter(TID == "24")

sr.25   <- sr.targets %>% filter(TID == "25")
sr.26   <- sr.targets %>% filter(TID == "26")

sr.27   <- sr.targets %>% filter(TID == "27")

sr.28   <- sr.targets %>% filter(TID == "28")

sr.29   <- sr.targets %>% filter(TID == "29")

sr.30   <- sr.targets %>% filter(TID == "30")
sr.31   <- sr.targets %>% filter(TID == "31")

sr.32   <- sr.targets %>% filter(TID == "32")

sr.33   <- sr.targets %>% filter(TID == "33")

sr.34   <- sr.targets %>% filter(TID == "34")

sr.35   <- sr.targets %>% filter(TID == "35")
sr.36   <- sr.targets %>% filter(TID == "36")

sr.37   <- sr.targets %>% filter(TID == "37")

sr.38   <- sr.targets %>% filter(TID == "38")

sr.39   <- sr.targets %>% filter(TID == "39")

sr.40   <- sr.targets %>% filter(TID == "40")
sr.41   <- sr.targets %>% filter(TID == "41")

sr.42   <- sr.targets %>% filter(TID == "42")

sr.43   <- sr.targets %>% filter(TID == "43")

sr.44   <- sr.targets %>% filter(TID == "44")

sr.45   <- sr.targets %>% filter(TID == "45")
sr.46   <- sr.targets %>% filter(TID == "46")

sr.47   <- sr.targets %>% filter(TID == "47")

sr.48   <- sr.targets %>% filter(TID == "48")

sr.49   <- sr.targets %>% filter(TID == "49")

sr.50   <- sr.targets %>% filter(TID == "50")
sr.51   <- sr.targets %>% filter(TID == "51")

sr.52   <- sr.targets %>% filter(TID == "52")

sr.53   <- sr.targets %>% filter(TID == "53")

sr.54   <- sr.targets %>% filter(TID == "54")

sr.55   <- sr.targets %>% filter(TID == "55")
sr.56   <- sr.targets %>% filter(TID == "56")

sr.57   <- sr.targets %>% filter(TID == "57")

sr.58   <- sr.targets %>% filter(TID == "58")

sr.59   <- sr.targets %>% filter(TID == "59")

sr.60   <- sr.targets %>% filter(TID == "60")
sr.61   <- sr.targets %>% filter(TID == "61")

sr.62   <- sr.targets %>% filter(TID == "62")

sr.63   <- sr.targets %>% filter(TID == "63")

sr.64   <- sr.targets %>% filter(TID == "64")

sr.65   <- sr.targets %>% filter(TID == "65")
sr.66   <- sr.targets %>% filter(TID == "66")

sr.67   <- sr.targets %>% filter(TID == "67")

sr.68   <- sr.targets %>% filter(TID == "68")

sr.69   <- sr.targets %>% filter(TID == "69")

sr.70   <- sr.targets %>% filter(TID == "70")
sr.71   <- sr.targets %>% filter(TID == "71")

sr.72   <- sr.targets %>% filter(TID == "72")

sr.73   <- sr.targets %>% filter(TID == "73")

sr.74   <- sr.targets %>% filter(TID == "74")

sr.75   <- sr.targets %>% filter(TID == "75")
sr.76   <- sr.targets %>% filter(TID == "76")

sr.77   <- sr.targets %>% filter(TID == "77")

sr.78   <- sr.targets %>% filter(TID == "78")

sr.79   <- sr.targets %>% filter(TID == "79")

sr.80   <- sr.targets %>% filter(TID == "80")
sr.81   <- sr.targets %>% filter(TID == "81")

sr.82   <- sr.targets %>% filter(TID == "82")

sr.83   <- sr.targets %>% filter(TID == "83")

sr.84   <- sr.targets %>% filter(TID == "84")

sr.85   <- sr.targets %>% filter(TID == "85")
sr.86   <- sr.targets %>% filter(TID == "86")

sr.87   <- sr.targets %>% filter(TID == "87")

sr.88   <- sr.targets %>% filter(TID == "88")

sr.89   <- sr.targets %>% filter(TID == "89")

sr.90   <- sr.targets %>% filter(TID == "90")
sr.91   <- sr.targets %>% filter(TID == "91")

sr.92   <- sr.targets %>% filter(TID == "92")

sr.93   <- sr.targets %>% filter(TID == "93")

sr.94   <- sr.targets %>% filter(TID == "94")

sr.95   <- sr.targets %>% filter(TID == "95")
sr.96   <- sr.targets %>% filter(TID == "96")

sr.97   <- sr.targets %>% filter(TID == "97")

sr.98   <- sr.targets %>% filter(TID == "98")

sr.99   <- sr.targets %>% filter(TID == "99")

sr.100   <- sr.targets %>% filter(TID == "100")
sr.101   <- sr.targets %>% filter(TID == "101")

sr.102   <- sr.targets %>% filter(TID == "102")

sr.103   <- sr.targets %>% filter(TID == "103")

sr.104   <- sr.targets %>% filter(TID == "104")

sr.105   <- sr.targets %>% filter(TID == "105")
sr.106   <- sr.targets %>% filter(TID == "106")

sr.107   <- sr.targets %>% filter(TID == "107")

sr.108   <- sr.targets %>% filter(TID == "108")

sr.109   <- sr.targets %>% filter(TID == "109")

sr.110   <- sr.targets %>% filter(TID == "110")
sr.111   <- sr.targets %>% filter(TID == "111")
sr.112   <- sr.targets %>% filter(TID == "112")

sr.113   <- sr.targets %>% filter(TID == "113")

sr.114   <- sr.targets %>% filter(TID == "114")

sr.115   <- sr.targets %>% filter(TID == "115")

sr.116   <- sr.targets %>% filter(TID == "116")
sr.117   <- sr.targets %>% filter(TID == "117")

sr.118   <- sr.targets %>% filter(TID == "118")

sr.119   <- sr.targets %>% filter(TID == "119")

sr.120   <- sr.targets %>% filter(TID == "120")

sr.121   <- sr.targets %>% filter(TID == "121")
sr.122   <- sr.targets %>% filter(TID == "122")

sr.123   <- sr.targets %>% filter(TID == "123")

sr.124   <- sr.targets %>% filter(TID == "124")

sr.125   <- sr.targets %>% filter(TID == "125")

sr.126   <- sr.targets %>% filter(TID == "126")
sr.127   <- sr.targets %>% filter(TID == "127")

sr.128   <- sr.targets %>% filter(TID == "128")

sr.129   <- sr.targets %>% filter(TID == "129")

sr.130   <- sr.targets %>% filter(TID == "130")
sr.131   <- sr.targets %>% filter(TID == "131")
sr.132   <- sr.targets %>% filter(TID == "132")

sr.133   <- sr.targets %>% filter(TID == "133")

sr.134   <- sr.targets %>% filter(TID == "134")

sr.135   <- sr.targets %>% filter(TID == "135")

sr.136   <- sr.targets %>% filter(TID == "136")
sr.137   <- sr.targets %>% filter(TID == "137")

sr.138   <- sr.targets %>% filter(TID == "138")

sr.139   <- sr.targets %>% filter(TID == "139")

sr.140   <- sr.targets %>% filter(TID == "140")
sr.141   <- sr.targets %>% filter(TID == "141")
sr.142   <- sr.targets %>% filter(TID == "142")

sr.143   <- sr.targets %>% filter(TID == "143")

sr.144   <- sr.targets %>% filter(TID == "144")

sr.145   <- sr.targets %>% filter(TID == "145")

sr.146   <- sr.targets %>% filter(TID == "146")
sr.147   <- sr.targets %>% filter(TID == "147")

sr.148   <- sr.targets %>% filter(TID == "148")

sr.149   <- sr.targets %>% filter(TID == "149")

sr.150   <- sr.targets %>% filter(TID == "150")

# This is the BFI normative profile, which is the averaged realistic accuracy criteria from the study
normbfi = read.table("NormativeBFI2.csv", fill=TRUE, header = T, sep = ",")
longcrit.1.bfi2 <- melt(sr.1, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.2.bfi2 <- melt(sr.2, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.3.bfi2 <- melt(sr.3, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.4.bfi2 <- melt(sr.4, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.5.bfi2 <- melt(sr.5, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.6.bfi2 <- melt(sr.6, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.7.bfi2 <- melt(sr.7, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.8.bfi2 <- melt(sr.8, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.9.bfi2 <- melt(sr.9, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.10.bfi2 <- melt(sr.10, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.11.bfi2 <- melt(sr.11, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.12.bfi2 <- melt(sr.12, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.13.bfi2 <- melt(sr.13, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.14.bfi2 <- melt(sr.14, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.15.bfi2 <- melt(sr.15, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.16.bfi2 <- melt(sr.16, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.17.bfi2 <- melt(sr.17, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.18.bfi2 <- melt(sr.18, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.19.bfi2 <- melt(sr.19, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.20.bfi2 <- melt(sr.20, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.21.bfi2 <- melt(sr.21, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.22.bfi2 <- melt(sr.22, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.23.bfi2 <- melt(sr.23, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.24.bfi2 <- melt(sr.24, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.25.bfi2 <- melt(sr.25, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.26.bfi2 <- melt(sr.26, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.27.bfi2 <- melt(sr.27, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.28.bfi2 <- melt(sr.28, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.29.bfi2 <- melt(sr.29, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.30.bfi2 <- melt(sr.30, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.31.bfi2 <- melt(sr.31, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.32.bfi2 <- melt(sr.32, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.33.bfi2 <- melt(sr.33, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.34.bfi2 <- melt(sr.34, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.35.bfi2 <- melt(sr.35, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.36.bfi2 <- melt(sr.36, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.37.bfi2 <- melt(sr.37, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.38.bfi2 <- melt(sr.38, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.39.bfi2 <- melt(sr.39, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.40.bfi2 <- melt(sr.40, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.41.bfi2 <- melt(sr.41, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.42.bfi2 <- melt(sr.42, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.43.bfi2 <- melt(sr.43, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.44.bfi2 <- melt(sr.44, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.45.bfi2 <- melt(sr.45, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.46.bfi2 <- melt(sr.46, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.47.bfi2 <- melt(sr.47, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.48.bfi2 <- melt(sr.48, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.49.bfi2 <- melt(sr.49, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.50.bfi2 <- melt(sr.50, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.51.bfi2 <- melt(sr.51, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.52.bfi2 <- melt(sr.52, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.53.bfi2 <- melt(sr.53, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.54.bfi2 <- melt(sr.54, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.55.bfi2 <- melt(sr.55, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.56.bfi2 <- melt(sr.56, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.57.bfi2 <- melt(sr.57, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.58.bfi2 <- melt(sr.58, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.59.bfi2 <- melt(sr.59, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.60.bfi2 <- melt(sr.60, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.61.bfi2 <- melt(sr.61, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.62.bfi2 <- melt(sr.62, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.63.bfi2 <- melt(sr.63, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.64.bfi2 <- melt(sr.64, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.65.bfi2 <- melt(sr.65, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.66.bfi2 <- melt(sr.66, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.67.bfi2 <- melt(sr.67, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.68.bfi2 <- melt(sr.68, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.69.bfi2 <- melt(sr.69, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.70.bfi2 <- melt(sr.70, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.71.bfi2 <- melt(sr.71, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.72.bfi2 <- melt(sr.72, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.73.bfi2 <- melt(sr.73, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.74.bfi2 <- melt(sr.74, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.75.bfi2 <- melt(sr.75, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.76.bfi2 <- melt(sr.76, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.77.bfi2 <- melt(sr.77, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.78.bfi2 <- melt(sr.78, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.79.bfi2 <- melt(sr.79, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.80.bfi2 <- melt(sr.80, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.81.bfi2 <- melt(sr.81, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.82.bfi2 <- melt(sr.82, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.83.bfi2 <- melt(sr.83, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.84.bfi2 <- melt(sr.84, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.85.bfi2 <- melt(sr.85, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.86.bfi2 <- melt(sr.86, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.87.bfi2 <- melt(sr.87, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.88.bfi2 <- melt(sr.88, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.89.bfi2 <- melt(sr.89, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.90.bfi2 <- melt(sr.90, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.91.bfi2 <- melt(sr.91, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.92.bfi2 <- melt(sr.92, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.93.bfi2 <- melt(sr.93, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.94.bfi2 <- melt(sr.94, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.95.bfi2 <- melt(sr.95, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.96.bfi2 <- melt(sr.96, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.97.bfi2 <- melt(sr.97, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.98.bfi2 <- melt(sr.98, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.99.bfi2 <- melt(sr.99, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.100.bfi2 <- melt(sr.100, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.101.bfi2 <- melt(sr.101, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.102.bfi2 <- melt(sr.102, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.103.bfi2 <- melt(sr.103, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.104.bfi2 <- melt(sr.104, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.105.bfi2 <- melt(sr.105, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.106.bfi2 <- melt(sr.106, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.107.bfi2 <- melt(sr.107, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.108.bfi2 <- melt(sr.108, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.109.bfi2 <- melt(sr.109, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.110.bfi2 <- melt(sr.110, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.111.bfi2 <- melt(sr.111, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.112.bfi2 <- melt(sr.112, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.113.bfi2 <- melt(sr.113, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.114.bfi2 <- melt(sr.114, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.115.bfi2 <- melt(sr.115, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.116.bfi2 <- melt(sr.116, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.117.bfi2 <- melt(sr.117, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.118.bfi2 <- melt(sr.118, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.119.bfi2 <- melt(sr.119, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.120.bfi2 <- melt(sr.120, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.121.bfi2 <- melt(sr.121, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.122.bfi2 <- melt(sr.122, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.123.bfi2 <- melt(sr.123, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.124.bfi2 <- melt(sr.124, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.125.bfi2 <- melt(sr.125, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.126.bfi2 <- melt(sr.126, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.127.bfi2 <- melt(sr.127, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.128.bfi2 <- melt(sr.128, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.129.bfi2 <- melt(sr.129, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.130.bfi2 <- melt(sr.130, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.131.bfi2 <- melt(sr.131, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.132.bfi2 <- melt(sr.132, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.133.bfi2 <- melt(sr.133, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.134.bfi2 <- melt(sr.134, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.135.bfi2 <- melt(sr.135, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.136.bfi2 <- melt(sr.136, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.137.bfi2 <- melt(sr.137, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.138.bfi2 <- melt(sr.138, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.139.bfi2 <- melt(sr.139, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.140.bfi2 <- melt(sr.140, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.141.bfi2 <- melt(sr.141, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.142.bfi2 <- melt(sr.142, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.143.bfi2 <- melt(sr.143, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.144.bfi2 <- melt(sr.144, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.145.bfi2 <- melt(sr.145, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.146.bfi2 <- melt(sr.146, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.147.bfi2 <- melt(sr.147, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.148.bfi2 <- melt(sr.148, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.149.bfi2 <- melt(sr.149, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longcrit.150.bfi2 <- melt(sr.150, id.vars = c(1), measure.vars = c(2:31), variable.name = "Sr.item", value.name = "Sr")
longnorm.bfi2 <- melt(normbfi, measure.vars = c(1:30), variable.name = "Norm.item", value.name = "Norm")
longratings.1.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2:31), variable.name = "Rating.item", value.name = "Rating")
longratings.1.bfi2 <- longratings.1.bfi2[order(longratings.1.bfi2$JID, longratings.1.bfi2$Rating.item),]

longratings.2.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(32:61), variable.name = "Rating.item", value.name = "Rating")
longratings.2.bfi2 <- longratings.2.bfi2[order(longratings.2.bfi2$JID, longratings.2.bfi2$Rating.item),]

longratings.3.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(62:91), variable.name = "Rating.item", value.name = "Rating")
longratings.3.bfi2 <- longratings.3.bfi2[order(longratings.3.bfi2$JID, longratings.3.bfi2$Rating.item),]

longratings.4.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(92:121), variable.name = "Rating.item", value.name = "Rating")
longratings.4.bfi2 <- longratings.4.bfi2[order(longratings.4.bfi2$JID, longratings.4.bfi2$Rating.item),]

longratings.5.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(122:151), variable.name = "Rating.item", value.name = "Rating")
longratings.5.bfi2 <- longratings.5.bfi2[order(longratings.5.bfi2$JID, longratings.5.bfi2$Rating.item),]

longratings.6.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(152:181), variable.name = "Rating.item", value.name = "Rating")
longratings.6.bfi2 <- longratings.6.bfi2[order(longratings.6.bfi2$JID, longratings.6.bfi2$Rating.item),]

longratings.7.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(182:211), variable.name = "Rating.item", value.name = "Rating")
longratings.7.bfi2 <- longratings.7.bfi2[order(longratings.7.bfi2$JID, longratings.7.bfi2$Rating.item),]

longratings.8.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(212:241), variable.name = "Rating.item", value.name = "Rating")
longratings.8.bfi2 <- longratings.8.bfi2[order(longratings.8.bfi2$JID, longratings.8.bfi2$Rating.item),]

longratings.9.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(242:271), variable.name = "Rating.item", value.name = "Rating")
longratings.9.bfi2 <- longratings.9.bfi2[order(longratings.9.bfi2$JID, longratings.9.bfi2$Rating.item),]

longratings.10.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(272:301), variable.name = "Rating.item", value.name = "Rating")
longratings.10.bfi2 <- longratings.10.bfi2[order(longratings.10.bfi2$JID, longratings.10.bfi2$Rating.item),]

longratings.11.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(302:331), variable.name = "Rating.item", value.name = "Rating")
longratings.11.bfi2 <- longratings.11.bfi2[order(longratings.11.bfi2$JID, longratings.11.bfi2$Rating.item),]

longratings.12.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(332:361), variable.name = "Rating.item", value.name = "Rating")
longratings.12.bfi2 <- longratings.12.bfi2[order(longratings.12.bfi2$JID, longratings.12.bfi2$Rating.item),]

longratings.13.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(362:391), variable.name = "Rating.item", value.name = "Rating")
longratings.13.bfi2 <- longratings.13.bfi2[order(longratings.13.bfi2$JID, longratings.13.bfi2$Rating.item),]

longratings.14.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(392:421), variable.name = "Rating.item", value.name = "Rating")
longratings.14.bfi2 <- longratings.14.bfi2[order(longratings.14.bfi2$JID, longratings.14.bfi2$Rating.item),]

longratings.15.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(422:451), variable.name = "Rating.item", value.name = "Rating")
longratings.15.bfi2 <- longratings.15.bfi2[order(longratings.15.bfi2$JID, longratings.15.bfi2$Rating.item),]

longratings.16.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(452:481), variable.name = "Rating.item", value.name = "Rating")
longratings.16.bfi2 <- longratings.16.bfi2[order(longratings.16.bfi2$JID, longratings.16.bfi2$Rating.item),]

longratings.17.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(482:511), variable.name = "Rating.item", value.name = "Rating")
longratings.17.bfi2 <- longratings.17.bfi2[order(longratings.17.bfi2$JID, longratings.17.bfi2$Rating.item),]

longratings.18.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(512:541), variable.name = "Rating.item", value.name = "Rating")
longratings.18.bfi2 <- longratings.18.bfi2[order(longratings.18.bfi2$JID, longratings.18.bfi2$Rating.item),]

longratings.19.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(542:571), variable.name = "Rating.item", value.name = "Rating")
longratings.19.bfi2 <- longratings.19.bfi2[order(longratings.19.bfi2$JID, longratings.19.bfi2$Rating.item),]

longratings.20.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(572:601), variable.name = "Rating.item", value.name = "Rating")
longratings.20.bfi2 <- longratings.20.bfi2[order(longratings.20.bfi2$JID, longratings.20.bfi2$Rating.item),]

longratings.21.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(602:631), variable.name = "Rating.item", value.name = "Rating")
longratings.21.bfi2 <- longratings.21.bfi2[order(longratings.21.bfi2$JID, longratings.21.bfi2$Rating.item),]

longratings.22.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(632:661), variable.name = "Rating.item", value.name = "Rating")
longratings.22.bfi2 <- longratings.22.bfi2[order(longratings.22.bfi2$JID, longratings.22.bfi2$Rating.item),]

longratings.23.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(662:691), variable.name = "Rating.item", value.name = "Rating")
longratings.23.bfi2 <- longratings.23.bfi2[order(longratings.23.bfi2$JID, longratings.23.bfi2$Rating.item),]

longratings.24.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(692:721), variable.name = "Rating.item", value.name = "Rating")
longratings.24.bfi2 <- longratings.24.bfi2[order(longratings.24.bfi2$JID, longratings.24.bfi2$Rating.item),]

longratings.25.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(722:751), variable.name = "Rating.item", value.name = "Rating")
longratings.25.bfi2 <- longratings.25.bfi2[order(longratings.25.bfi2$JID, longratings.25.bfi2$Rating.item),]

longratings.26.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(752:781), variable.name = "Rating.item", value.name = "Rating")
longratings.26.bfi2 <- longratings.26.bfi2[order(longratings.26.bfi2$JID, longratings.26.bfi2$Rating.item),]

longratings.27.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(782:811), variable.name = "Rating.item", value.name = "Rating")
longratings.27.bfi2 <- longratings.27.bfi2[order(longratings.27.bfi2$JID, longratings.27.bfi2$Rating.item),]

longratings.28.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(812:841), variable.name = "Rating.item", value.name = "Rating")
longratings.28.bfi2 <- longratings.28.bfi2[order(longratings.28.bfi2$JID, longratings.28.bfi2$Rating.item),]

longratings.29.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(842:871), variable.name = "Rating.item", value.name = "Rating")
longratings.29.bfi2 <- longratings.29.bfi2[order(longratings.29.bfi2$JID, longratings.29.bfi2$Rating.item),]

longratings.30.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(872:901), variable.name = "Rating.item", value.name = "Rating")
longratings.30.bfi2 <- longratings.30.bfi2[order(longratings.30.bfi2$JID, longratings.30.bfi2$Rating.item),]

longratings.31.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(902:931), variable.name = "Rating.item", value.name = "Rating")
longratings.31.bfi2 <- longratings.31.bfi2[order(longratings.31.bfi2$JID, longratings.31.bfi2$Rating.item),]

longratings.32.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(932:961), variable.name = "Rating.item", value.name = "Rating")
longratings.32.bfi2 <- longratings.32.bfi2[order(longratings.32.bfi2$JID, longratings.32.bfi2$Rating.item),]

longratings.33.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(962:991), variable.name = "Rating.item", value.name = "Rating")
longratings.33.bfi2 <- longratings.33.bfi2[order(longratings.33.bfi2$JID, longratings.33.bfi2$Rating.item),]

longratings.34.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(992:1021), variable.name = "Rating.item", value.name = "Rating")
longratings.34.bfi2 <- longratings.34.bfi2[order(longratings.34.bfi2$JID, longratings.34.bfi2$Rating.item),]

longratings.35.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1022:1051), variable.name = "Rating.item", value.name = "Rating")
longratings.35.bfi2 <- longratings.35.bfi2[order(longratings.35.bfi2$JID, longratings.35.bfi2$Rating.item),]

longratings.36.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1052:1081), variable.name = "Rating.item", value.name = "Rating")
longratings.36.bfi2 <- longratings.36.bfi2[order(longratings.36.bfi2$JID, longratings.36.bfi2$Rating.item),]

longratings.37.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1082:1111), variable.name = "Rating.item", value.name = "Rating")
longratings.37.bfi2 <- longratings.37.bfi2[order(longratings.37.bfi2$JID, longratings.37.bfi2$Rating.item),]

longratings.38.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1112:1141), variable.name = "Rating.item", value.name = "Rating")
longratings.38.bfi2 <- longratings.38.bfi2[order(longratings.38.bfi2$JID, longratings.38.bfi2$Rating.item),]

longratings.39.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1142:1171), variable.name = "Rating.item", value.name = "Rating")
longratings.39.bfi2 <- longratings.39.bfi2[order(longratings.39.bfi2$JID, longratings.39.bfi2$Rating.item),]

longratings.40.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1172:1201), variable.name = "Rating.item", value.name = "Rating")
longratings.40.bfi2 <- longratings.40.bfi2[order(longratings.40.bfi2$JID, longratings.40.bfi2$Rating.item),]

longratings.41.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1202:1231), variable.name = "Rating.item", value.name = "Rating")
longratings.41.bfi2 <- longratings.1.bfi2[order(longratings.41.bfi2$JID, longratings.41.bfi2$Rating.item),]

longratings.42.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1232:1261), variable.name = "Rating.item", value.name = "Rating")
longratings.42.bfi2 <- longratings.42.bfi2[order(longratings.42.bfi2$JID, longratings.42.bfi2$Rating.item),]

longratings.43.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1262:1291), variable.name = "Rating.item", value.name = "Rating")
longratings.43.bfi2 <- longratings.43.bfi2[order(longratings.43.bfi2$JID, longratings.43.bfi2$Rating.item),]

longratings.44.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1292:1321), variable.name = "Rating.item", value.name = "Rating")
longratings.44.bfi2 <- longratings.44.bfi2[order(longratings.44.bfi2$JID, longratings.44.bfi2$Rating.item),]

longratings.45.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1322:1351), variable.name = "Rating.item", value.name = "Rating")
longratings.45.bfi2 <- longratings.45.bfi2[order(longratings.45.bfi2$JID, longratings.45.bfi2$Rating.item),]

longratings.46.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1352:1381), variable.name = "Rating.item", value.name = "Rating")
longratings.46.bfi2 <- longratings.46.bfi2[order(longratings.46.bfi2$JID, longratings.46.bfi2$Rating.item),]

longratings.47.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1382:1411), variable.name = "Rating.item", value.name = "Rating")
longratings.47.bfi2 <- longratings.47.bfi2[order(longratings.47.bfi2$JID, longratings.47.bfi2$Rating.item),]

longratings.48.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1412:1441), variable.name = "Rating.item", value.name = "Rating")
longratings.48.bfi2 <- longratings.48.bfi2[order(longratings.48.bfi2$JID, longratings.48.bfi2$Rating.item),]

longratings.49.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1442:1471), variable.name = "Rating.item", value.name = "Rating")
longratings.49.bfi2 <- longratings.49.bfi2[order(longratings.49.bfi2$JID, longratings.49.bfi2$Rating.item),]

longratings.50.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1472:1501), variable.name = "Rating.item", value.name = "Rating")
longratings.50.bfi2 <- longratings.50.bfi2[order(longratings.50.bfi2$JID, longratings.50.bfi2$Rating.item),]

longratings.51.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1502:1531), variable.name = "Rating.item", value.name = "Rating")
longratings.51.bfi2 <- longratings.51.bfi2[order(longratings.51.bfi2$JID, longratings.51.bfi2$Rating.item),]

longratings.52.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1532:1561), variable.name = "Rating.item", value.name = "Rating")
longratings.52.bfi2 <- longratings.52.bfi2[order(longratings.52.bfi2$JID, longratings.52.bfi2$Rating.item),]

longratings.53.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1562:1591), variable.name = "Rating.item", value.name = "Rating")
longratings.53.bfi2 <- longratings.53.bfi2[order(longratings.53.bfi2$JID, longratings.53.bfi2$Rating.item),]

longratings.54.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1592:1621), variable.name = "Rating.item", value.name = "Rating")
longratings.54.bfi2 <- longratings.54.bfi2[order(longratings.54.bfi2$JID, longratings.54.bfi2$Rating.item),]

longratings.55.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1622:1651), variable.name = "Rating.item", value.name = "Rating")
longratings.55.bfi2 <- longratings.55.bfi2[order(longratings.55.bfi2$JID, longratings.55.bfi2$Rating.item),]

longratings.56.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1652:1681), variable.name = "Rating.item", value.name = "Rating")
longratings.56.bfi2 <- longratings.56.bfi2[order(longratings.56.bfi2$JID, longratings.56.bfi2$Rating.item),]

longratings.57.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1682:1711), variable.name = "Rating.item", value.name = "Rating")
longratings.57.bfi2 <- longratings.57.bfi2[order(longratings.57.bfi2$JID, longratings.57.bfi2$Rating.item),]

longratings.58.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1712:1741), variable.name = "Rating.item", value.name = "Rating")
longratings.58.bfi2 <- longratings.58.bfi2[order(longratings.58.bfi2$JID, longratings.58.bfi2$Rating.item),]

longratings.59.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1742:1771), variable.name = "Rating.item", value.name = "Rating")
longratings.59.bfi2 <- longratings.59.bfi2[order(longratings.59.bfi2$JID, longratings.59.bfi2$Rating.item),]

longratings.60.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1772:1801), variable.name = "Rating.item", value.name = "Rating")
longratings.60.bfi2 <- longratings.60.bfi2[order(longratings.60.bfi2$JID, longratings.60.bfi2$Rating.item),]

longratings.61.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1802:1831), variable.name = "Rating.item", value.name = "Rating")
longratings.61.bfi2 <- longratings.61.bfi2[order(longratings.61.bfi2$JID, longratings.61.bfi2$Rating.item),]

longratings.62.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1832:1861), variable.name = "Rating.item", value.name = "Rating")
longratings.62.bfi2 <- longratings.62.bfi2[order(longratings.62.bfi2$JID, longratings.62.bfi2$Rating.item),]

longratings.63.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1862:1891), variable.name = "Rating.item", value.name = "Rating")
longratings.63.bfi2 <- longratings.63.bfi2[order(longratings.63.bfi2$JID, longratings.63.bfi2$Rating.item),]

longratings.64.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1892:1921), variable.name = "Rating.item", value.name = "Rating")
longratings.64.bfi2 <- longratings.64.bfi2[order(longratings.64.bfi2$JID, longratings.64.bfi2$Rating.item),]

longratings.65.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1922:1951), variable.name = "Rating.item", value.name = "Rating")
longratings.65.bfi2 <- longratings.65.bfi2[order(longratings.65.bfi2$JID, longratings.65.bfi2$Rating.item),]

longratings.66.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1952:1981), variable.name = "Rating.item", value.name = "Rating")
longratings.66.bfi2 <- longratings.66.bfi2[order(longratings.66.bfi2$JID, longratings.66.bfi2$Rating.item),]

longratings.67.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(1982:2011), variable.name = "Rating.item", value.name = "Rating")
longratings.67.bfi2 <- longratings.67.bfi2[order(longratings.67.bfi2$JID, longratings.67.bfi2$Rating.item),]

longratings.68.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2012:2041), variable.name = "Rating.item", value.name = "Rating")
longratings.68.bfi2 <- longratings.68.bfi2[order(longratings.68.bfi2$JID, longratings.68.bfi2$Rating.item),]

longratings.69.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2042:2071), variable.name = "Rating.item", value.name = "Rating")
longratings.69.bfi2 <- longratings.69.bfi2[order(longratings.69.bfi2$JID, longratings.69.bfi2$Rating.item),]

longratings.70.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2072:2101), variable.name = "Rating.item", value.name = "Rating")
longratings.70.bfi2 <- longratings.70.bfi2[order(longratings.70.bfi2$JID, longratings.70.bfi2$Rating.item),]

longratings.71.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2102:2131), variable.name = "Rating.item", value.name = "Rating")
longratings.71.bfi2 <- longratings.71.bfi2[order(longratings.71.bfi2$JID, longratings.71.bfi2$Rating.item),]

longratings.72.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2132:2161), variable.name = "Rating.item", value.name = "Rating")
longratings.72.bfi2 <- longratings.72.bfi2[order(longratings.72.bfi2$JID, longratings.72.bfi2$Rating.item),]

longratings.73.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2162:2191), variable.name = "Rating.item", value.name = "Rating")
longratings.73.bfi2 <- longratings.73.bfi2[order(longratings.73.bfi2$JID, longratings.73.bfi2$Rating.item),]

longratings.74.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2192:2221), variable.name = "Rating.item", value.name = "Rating")
longratings.74.bfi2 <- longratings.74.bfi2[order(longratings.74.bfi2$JID, longratings.74.bfi2$Rating.item),]

longratings.75.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2222:2251), variable.name = "Rating.item", value.name = "Rating")
longratings.75.bfi2 <- longratings.75.bfi2[order(longratings.75.bfi2$JID, longratings.75.bfi2$Rating.item),]

longratings.76.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2252:2281), variable.name = "Rating.item", value.name = "Rating")
longratings.76.bfi2 <- longratings.76.bfi2[order(longratings.76.bfi2$JID, longratings.76.bfi2$Rating.item),]

longratings.77.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2282:2311), variable.name = "Rating.item", value.name = "Rating")
longratings.77.bfi2 <- longratings.77.bfi2[order(longratings.77.bfi2$JID, longratings.77.bfi2$Rating.item),]

longratings.78.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2312:2341), variable.name = "Rating.item", value.name = "Rating")
longratings.78.bfi2 <- longratings.78.bfi2[order(longratings.78.bfi2$JID, longratings.78.bfi2$Rating.item),]

longratings.79.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2342:2371), variable.name = "Rating.item", value.name = "Rating")
longratings.79.bfi2 <- longratings.79.bfi2[order(longratings.79.bfi2$JID, longratings.79.bfi2$Rating.item),]

longratings.80.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2372:2401), variable.name = "Rating.item", value.name = "Rating")
longratings.80.bfi2 <- longratings.80.bfi2[order(longratings.80.bfi2$JID, longratings.80.bfi2$Rating.item),]

longratings.81.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2402:2431), variable.name = "Rating.item", value.name = "Rating")
longratings.81.bfi2 <- longratings.81.bfi2[order(longratings.81.bfi2$JID, longratings.81.bfi2$Rating.item),]

longratings.82.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2432:2461), variable.name = "Rating.item", value.name = "Rating")
longratings.82.bfi2 <- longratings.82.bfi2[order(longratings.82.bfi2$JID, longratings.82.bfi2$Rating.item),]

longratings.83.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2462:2491), variable.name = "Rating.item", value.name = "Rating")
longratings.83.bfi2 <- longratings.83.bfi2[order(longratings.83.bfi2$JID, longratings.83.bfi2$Rating.item),]

longratings.84.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2492:2521), variable.name = "Rating.item", value.name = "Rating")
longratings.84.bfi2 <- longratings.84.bfi2[order(longratings.84.bfi2$JID, longratings.84.bfi2$Rating.item),]

longratings.85.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2522:2551), variable.name = "Rating.item", value.name = "Rating")
longratings.85.bfi2 <- longratings.85.bfi2[order(longratings.85.bfi2$JID, longratings.85.bfi2$Rating.item),]

longratings.86.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2552:2581), variable.name = "Rating.item", value.name = "Rating")
longratings.86.bfi2 <- longratings.86.bfi2[order(longratings.86.bfi2$JID, longratings.86.bfi2$Rating.item),]

longratings.87.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2582:2611), variable.name = "Rating.item", value.name = "Rating")
longratings.87.bfi2 <- longratings.87.bfi2[order(longratings.87.bfi2$JID, longratings.87.bfi2$Rating.item),]

longratings.88.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2612:2641), variable.name = "Rating.item", value.name = "Rating")
longratings.88.bfi2 <- longratings.88.bfi2[order(longratings.88.bfi2$JID, longratings.88.bfi2$Rating.item),]

longratings.89.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2462:2671), variable.name = "Rating.item", value.name = "Rating")
longratings.89.bfi2 <- longratings.89.bfi2[order(longratings.89.bfi2$JID, longratings.89.bfi2$Rating.item),]

longratings.90.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2672:2701), variable.name = "Rating.item", value.name = "Rating")
longratings.90.bfi2 <- longratings.90.bfi2[order(longratings.90.bfi2$JID, longratings.90.bfi2$Rating.item),]

longratings.91.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2702:2731), variable.name = "Rating.item", value.name = "Rating")
longratings.91.bfi2 <- longratings.91.bfi2[order(longratings.91.bfi2$JID, longratings.91.bfi2$Rating.item),]

longratings.92.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2732:2761), variable.name = "Rating.item", value.name = "Rating")
longratings.92.bfi2 <- longratings.92.bfi2[order(longratings.92.bfi2$JID, longratings.92.bfi2$Rating.item),]

longratings.93.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2762:2791), variable.name = "Rating.item", value.name = "Rating")
longratings.93.bfi2 <- longratings.93.bfi2[order(longratings.93.bfi2$JID, longratings.93.bfi2$Rating.item),]

longratings.94.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2792:2821), variable.name = "Rating.item", value.name = "Rating")
longratings.94.bfi2 <- longratings.94.bfi2[order(longratings.94.bfi2$JID, longratings.94.bfi2$Rating.item),]

longratings.95.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2822:2851), variable.name = "Rating.item", value.name = "Rating")
longratings.95.bfi2 <- longratings.95.bfi2[order(longratings.95.bfi2$JID, longratings.95.bfi2$Rating.item),]

longratings.96.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2852:2881), variable.name = "Rating.item", value.name = "Rating")
longratings.96.bfi2 <- longratings.96.bfi2[order(longratings.96.bfi2$JID, longratings.96.bfi2$Rating.item),]

longratings.97.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2882:2911), variable.name = "Rating.item", value.name = "Rating")
longratings.97.bfi2 <- longratings.97.bfi2[order(longratings.97.bfi2$JID, longratings.97.bfi2$Rating.item),]

longratings.98.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2912:2941), variable.name = "Rating.item", value.name = "Rating")
longratings.98.bfi2 <- longratings.98.bfi2[order(longratings.98.bfi2$JID, longratings.98.bfi2$Rating.item),]

longratings.99.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2942:2971), variable.name = "Rating.item", value.name = "Rating")
longratings.99.bfi2 <- longratings.99.bfi2[order(longratings.99.bfi2$JID, longratings.99.bfi2$Rating.item),]

longratings.100.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(2972:3001), variable.name = "Rating.item", value.name = "Rating")
longratings.100.bfi2 <- longratings.100.bfi2[order(longratings.100.bfi2$JID, longratings.100.bfi2$Rating.item),]

longratings.101.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3002:3031), variable.name = "Rating.item", value.name = "Rating")
longratings.101.bfi2 <- longratings.101.bfi2[order(longratings.101.bfi2$JID, longratings.101.bfi2$Rating.item),]

longratings.102.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3032:3061), variable.name = "Rating.item", value.name = "Rating")
longratings.102.bfi2 <- longratings.102.bfi2[order(longratings.102.bfi2$JID, longratings.102.bfi2$Rating.item),]

longratings.103.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3062:3091), variable.name = "Rating.item", value.name = "Rating")
longratings.103.bfi2 <- longratings.103.bfi2[order(longratings.103.bfi2$JID, longratings.103.bfi2$Rating.item),]

longratings.104.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3092:3121), variable.name = "Rating.item", value.name = "Rating")
longratings.104.bfi2 <- longratings.104.bfi2[order(longratings.104.bfi2$JID, longratings.104.bfi2$Rating.item),]

longratings.105.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3122:3151), variable.name = "Rating.item", value.name = "Rating")
longratings.105.bfi2 <- longratings.105.bfi2[order(longratings.105.bfi2$JID, longratings.105.bfi2$Rating.item),]

longratings.106.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3152:3181), variable.name = "Rating.item", value.name = "Rating")
longratings.106.bfi2 <- longratings.106.bfi2[order(longratings.106.bfi2$JID, longratings.106.bfi2$Rating.item),]

longratings.107.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3182:3211), variable.name = "Rating.item", value.name = "Rating")
longratings.107.bfi2 <- longratings.107.bfi2[order(longratings.107.bfi2$JID, longratings.107.bfi2$Rating.item),]

longratings.108.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3212:3241), variable.name = "Rating.item", value.name = "Rating")
longratings.108.bfi2 <- longratings.108.bfi2[order(longratings.108.bfi2$JID, longratings.108.bfi2$Rating.item),]

longratings.109.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3242:3271), variable.name = "Rating.item", value.name = "Rating")
longratings.109.bfi2 <- longratings.109.bfi2[order(longratings.109.bfi2$JID, longratings.109.bfi2$Rating.item),]

longratings.110.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3272:3301), variable.name = "Rating.item", value.name = "Rating")
longratings.110.bfi2 <- longratings.110.bfi2[order(longratings.110.bfi2$JID, longratings.110.bfi2$Rating.item),]

longratings.111.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3302:3331), variable.name = "Rating.item", value.name = "Rating")
longratings.111.bfi2 <- longratings.111.bfi2[order(longratings.111.bfi2$JID, longratings.111.bfi2$Rating.item),]

longratings.112.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3332:3361), variable.name = "Rating.item", value.name = "Rating")
longratings.112.bfi2 <- longratings.112.bfi2[order(longratings.112.bfi2$JID, longratings.112.bfi2$Rating.item),]

longratings.113.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3362:3391), variable.name = "Rating.item", value.name = "Rating")
longratings.113.bfi2 <- longratings.113.bfi2[order(longratings.113.bfi2$JID, longratings.113.bfi2$Rating.item),]

longratings.114.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3392:3421), variable.name = "Rating.item", value.name = "Rating")
longratings.114.bfi2 <- longratings.114.bfi2[order(longratings.114.bfi2$JID, longratings.114.bfi2$Rating.item),]

longratings.115.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3422:3451), variable.name = "Rating.item", value.name = "Rating")
longratings.115.bfi2 <- longratings.115.bfi2[order(longratings.115.bfi2$JID, longratings.115.bfi2$Rating.item),]

longratings.116.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3452:3481), variable.name = "Rating.item", value.name = "Rating")
longratings.116.bfi2 <- longratings.116.bfi2[order(longratings.116.bfi2$JID, longratings.116.bfi2$Rating.item),]

longratings.117.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3482:3511), variable.name = "Rating.item", value.name = "Rating")
longratings.117.bfi2 <- longratings.117.bfi2[order(longratings.117.bfi2$JID, longratings.117.bfi2$Rating.item),]

longratings.118.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3512:3541), variable.name = "Rating.item", value.name = "Rating")
longratings.118.bfi2 <- longratings.118.bfi2[order(longratings.118.bfi2$JID, longratings.118.bfi2$Rating.item),]

longratings.119.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3542:3571), variable.name = "Rating.item", value.name = "Rating")
longratings.119.bfi2 <- longratings.119.bfi2[order(longratings.119.bfi2$JID, longratings.119.bfi2$Rating.item),]

longratings.120.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3572:3601), variable.name = "Rating.item", value.name = "Rating")
longratings.120.bfi2 <- longratings.120.bfi2[order(longratings.120.bfi2$JID, longratings.120.bfi2$Rating.item),]

longratings.121.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3602:3631), variable.name = "Rating.item", value.name = "Rating")
longratings.121.bfi2 <- longratings.121.bfi2[order(longratings.121.bfi2$JID, longratings.121.bfi2$Rating.item),]

longratings.122.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3632:3661), variable.name = "Rating.item", value.name = "Rating")
longratings.122.bfi2 <- longratings.112.bfi2[order(longratings.112.bfi2$JID, longratings.122.bfi2$Rating.item),]

longratings.123.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3662:3691), variable.name = "Rating.item", value.name = "Rating")
longratings.123.bfi2 <- longratings.123.bfi2[order(longratings.123.bfi2$JID, longratings.123.bfi2$Rating.item),]

longratings.124.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3692:3721), variable.name = "Rating.item", value.name = "Rating")
longratings.124.bfi2 <- longratings.124.bfi2[order(longratings.124.bfi2$JID, longratings.124.bfi2$Rating.item),]

longratings.125.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3722:3751), variable.name = "Rating.item", value.name = "Rating")
longratings.125.bfi2 <- longratings.125.bfi2[order(longratings.125.bfi2$JID, longratings.125.bfi2$Rating.item),]

longratings.126.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3752:3781), variable.name = "Rating.item", value.name = "Rating")
longratings.126.bfi2 <- longratings.126.bfi2[order(longratings.126.bfi2$JID, longratings.126.bfi2$Rating.item),]

longratings.127.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3782:3811), variable.name = "Rating.item", value.name = "Rating")
longratings.127.bfi2 <- longratings.127.bfi2[order(longratings.127.bfi2$JID, longratings.127.bfi2$Rating.item),]

longratings.128.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3812:3841), variable.name = "Rating.item", value.name = "Rating")
longratings.128.bfi2 <- longratings.128.bfi2[order(longratings.128.bfi2$JID, longratings.128.bfi2$Rating.item),]

longratings.129.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3842:3871), variable.name = "Rating.item", value.name = "Rating")
longratings.129.bfi2 <- longratings.129.bfi2[order(longratings.129.bfi2$JID, longratings.129.bfi2$Rating.item),]

longratings.130.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3872:3901), variable.name = "Rating.item", value.name = "Rating")
longratings.130.bfi2 <- longratings.130.bfi2[order(longratings.130.bfi2$JID, longratings.130.bfi2$Rating.item),]

longratings.131.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3902:3931), variable.name = "Rating.item", value.name = "Rating")
longratings.131.bfi2 <- longratings.131.bfi2[order(longratings.131.bfi2$JID, longratings.131.bfi2$Rating.item),]

longratings.132.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3932:3961), variable.name = "Rating.item", value.name = "Rating")
longratings.132.bfi2 <- longratings.132.bfi2[order(longratings.132.bfi2$JID, longratings.132.bfi2$Rating.item),]

longratings.133.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3962:3991), variable.name = "Rating.item", value.name = "Rating")
longratings.133.bfi2 <- longratings.133.bfi2[order(longratings.133.bfi2$JID, longratings.133.bfi2$Rating.item),]

longratings.134.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(3992:4021), variable.name = "Rating.item", value.name = "Rating")
longratings.134.bfi2 <- longratings.134.bfi2[order(longratings.134.bfi2$JID, longratings.134.bfi2$Rating.item),]

longratings.135.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4022:4051), variable.name = "Rating.item", value.name = "Rating")
longratings.135.bfi2 <- longratings.135.bfi2[order(longratings.135.bfi2$JID, longratings.135.bfi2$Rating.item),]

longratings.136.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4052:4081), variable.name = "Rating.item", value.name = "Rating")
longratings.136.bfi2 <- longratings.136.bfi2[order(longratings.136.bfi2$JID, longratings.136.bfi2$Rating.item),]

longratings.137.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4082:4111), variable.name = "Rating.item", value.name = "Rating")
longratings.137.bfi2 <- longratings.137.bfi2[order(longratings.137.bfi2$JID, longratings.137.bfi2$Rating.item),]

longratings.138.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4112:4141), variable.name = "Rating.item", value.name = "Rating")
longratings.138.bfi2 <- longratings.138.bfi2[order(longratings.138.bfi2$JID, longratings.138.bfi2$Rating.item),]

longratings.139.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4142:4171), variable.name = "Rating.item", value.name = "Rating")
longratings.139.bfi2 <- longratings.139.bfi2[order(longratings.139.bfi2$JID, longratings.139.bfi2$Rating.item),]

longratings.140.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4172:4201), variable.name = "Rating.item", value.name = "Rating")
longratings.140.bfi2 <- longratings.140.bfi2[order(longratings.140.bfi2$JID, longratings.140.bfi2$Rating.item),]

longratings.141.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4202:4231), variable.name = "Rating.item", value.name = "Rating")
longratings.141.bfi2 <- longratings.141.bfi2[order(longratings.141.bfi2$JID, longratings.141.bfi2$Rating.item),]

longratings.142.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4232:4261), variable.name = "Rating.item", value.name = "Rating")
longratings.142.bfi2 <- longratings.142.bfi2[order(longratings.142.bfi2$JID, longratings.142.bfi2$Rating.item),]

longratings.143.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4262:4291), variable.name = "Rating.item", value.name = "Rating")
longratings.143.bfi2 <- longratings.143.bfi2[order(longratings.143.bfi2$JID, longratings.143.bfi2$Rating.item),]

longratings.144.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4292:4321), variable.name = "Rating.item", value.name = "Rating")
longratings.144.bfi2 <- longratings.144.bfi2[order(longratings.144.bfi2$JID, longratings.144.bfi2$Rating.item),]

longratings.145.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4322:4351), variable.name = "Rating.item", value.name = "Rating")
longratings.145.bfi2 <- longratings.145.bfi2[order(longratings.145.bfi2$JID, longratings.145.bfi2$Rating.item),]

longratings.146.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4352:4381), variable.name = "Rating.item", value.name = "Rating")
longratings.146.bfi2 <- longratings.146.bfi2[order(longratings.146.bfi2$JID, longratings.146.bfi2$Rating.item),]

longratings.147.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4382:4411), variable.name = "Rating.item", value.name = "Rating")
longratings.147.bfi2 <- longratings.147.bfi2[order(longratings.147.bfi2$JID, longratings.147.bfi2$Rating.item),]

longratings.148.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4412:4441), variable.name = "Rating.item", value.name = "Rating")
longratings.148.bfi2 <- longratings.148.bfi2[order(longratings.148.bfi2$JID, longratings.148.bfi2$Rating.item),]

longratings.149.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4442:4471), variable.name = "Rating.item", value.name = "Rating")
longratings.149.bfi2 <- longratings.149.bfi2[order(longratings.149.bfi2$JID, longratings.149.bfi2$Rating.item),]

longratings.150.bfi2 <- melt(data, id.vars = c(1, 4502:4504), measure.vars = c(4472:4501), variable.name = "Rating.item", value.name = "Rating")
longratings.150.bfi2 <- longratings.150.bfi2[order(longratings.150.bfi2$JID, longratings.150.bfi2$Rating.item),]



SAM.2TW.bfi2 <- cbind(longratings.2.bfi2,longnorm.bfi2,longcrit.2.bfi2, TW.Cond, TW.2)

SAM.4TW.bfi2 <- cbind(longratings.4.bfi2,longnorm.bfi2,longcrit.4.bfi2, TW.Cond, TW.4)

SAM.6TW.bfi2 <- cbind(longratings.6.bfi2,longnorm.bfi2,longcrit.6.bfi2, TW.Cond, TW.6)
SAM.7TW.bfi2 <- cbind(longratings.7.bfi2,longnorm.bfi2,longcrit.7.bfi2, TW.Cond, TW.7)

SAM.9TW.bfi2 <- cbind(longratings.9.bfi2,longnorm.bfi2,longcrit.9.bfi2, TW.Cond, TW.9)

SAM.11TW.bfi2 <- cbind(longratings.11.bfi2,longnorm.bfi2,longcrit.11.bfi2, TW.Cond, TW.11)

SAM.14TW.bfi2 <- cbind(longratings.14.bfi2,longnorm.bfi2,longcrit.14.bfi2, TW.Cond, TW.14)

SAM.16TW.bfi2 <- cbind(longratings.16.bfi2,longnorm.bfi2,longcrit.16.bfi2, TW.Cond, TW.16)

SAM.18TW.bfi2 <- cbind(longratings.18.bfi2,longnorm.bfi2,longcrit.18.bfi2, TW.Cond, TW.18)

SAM.20TW.bfi2 <- cbind(longratings.20.bfi2,longnorm.bfi2,longcrit.20.bfi2, TW.Cond, TW.20)

SAM.22TW.bfi2 <- cbind(longratings.22.bfi2,longnorm.bfi2,longcrit.22.bfi2, TW.Cond, TW.22)

SAM.24TW.bfi2 <- cbind(longratings.24.bfi2,longnorm.bfi2,longcrit.24.bfi2, TW.Cond, TW.24)
SAM.25TW.bfi2 <- cbind(longratings.25.bfi2,longnorm.bfi2,longcrit.25.bfi2, TW.Cond, TW.25)

SAM.27TW.bfi2 <- cbind(longratings.27.bfi2,longnorm.bfi2,longcrit.27.bfi2, TW.Cond, TW.27)

SAM.29TW.bfi2 <- cbind(longratings.29.bfi2,longnorm.bfi2,longcrit.29.bfi2, TW.Cond, TW.29)

SAM.32TW.bfi2 <- cbind(longratings.32.bfi2,longnorm.bfi2,longcrit.32.bfi2, TW.Cond, TW.32)

SAM.34TW.bfi2 <- cbind(longratings.34.bfi2,longnorm.bfi2,longcrit.34.bfi2, TW.Cond, TW.34)

SAM.36TW.bfi2 <- cbind(longratings.36.bfi2,longnorm.bfi2,longcrit.36.bfi2, TW.Cond, TW.36)

SAM.38TW.bfi2 <- cbind(longratings.38.bfi2,longnorm.bfi2,longcrit.38.bfi2, TW.Cond, TW.38)

SAM.40TW.bfi2 <- cbind(longratings.40.bfi2,longnorm.bfi2,longcrit.40.bfi2, TW.Cond, TW.40)

SAM.42TW.bfi2 <- cbind(longratings.42.bfi2,longnorm.bfi2,longcrit.42.bfi2, TW.Cond, TW.42)
SAM.43TW.bfi2 <- cbind(longratings.43.bfi2,longnorm.bfi2,longcrit.43.bfi2, TW.Cond, TW.43)

SAM.45TW.bfi2 <- cbind(longratings.45.bfi2,longnorm.bfi2,longcrit.45.bfi2, TW.Cond, TW.45)

SAM.47TW.bfi2 <- cbind(longratings.47.bfi2,longnorm.bfi2,longcrit.47.bfi2, TW.Cond, TW.47)

SAM.50TW.bfi2 <- cbind(longratings.50.bfi2,longnorm.bfi2,longcrit.50.bfi2, TW.Cond, TW.50)

SAM.52TW.bfi2 <- cbind(longratings.52.bfi2,longnorm.bfi2,longcrit.52.bfi2, TW.Cond, TW.52)

SAM.54TW.bfi2 <- cbind(longratings.54.bfi2,longnorm.bfi2,longcrit.54.bfi2, TW.Cond, TW.54)

SAM.56TW.bfi2 <- cbind(longratings.56.bfi2,longnorm.bfi2,longcrit.56.bfi2, TW.Cond, TW.56)

SAM.58TW.bfi2 <- cbind(longratings.58.bfi2,longnorm.bfi2,longcrit.58.bfi2, TW.Cond, TW.58)

SAM.60TW.bfi2 <- cbind(longratings.60.bfi2,longnorm.bfi2,longcrit.60.bfi2, TW.Cond, TW.60)
SAM.61TW.bfi2 <- cbind(longratings.61.bfi2,longnorm.bfi2,longcrit.61.bfi2, TW.Cond, TW.61)

SAM.63TW.bfi2 <- cbind(longratings.63.bfi2,longnorm.bfi2,longcrit.63.bfi2, TW.Cond, TW.63)

SAM.65TW.bfi2 <- cbind(longratings.65.bfi2,longnorm.bfi2,longcrit.65.bfi2, TW.Cond, TW.65)

SAM.68TW.bfi2 <- cbind(longratings.68.bfi2,longnorm.bfi2,longcrit.68.bfi2, TW.Cond, TW.68)

SAM.70TW.bfi2 <- cbind(longratings.70.bfi2,longnorm.bfi2,longcrit.70.bfi2, TW.Cond, TW.70)

SAM.72TW.bfi2 <- cbind(longratings.72.bfi2,longnorm.bfi2,longcrit.72.bfi2, TW.Cond, TW.72)

SAM.74TW.bfi2 <- cbind(longratings.74.bfi2,longnorm.bfi2,longcrit.74.bfi2, TW.Cond, TW.74)

SAM.76TW.bfi2 <- cbind(longratings.76.bfi2,longnorm.bfi2,longcrit.76.bfi2, TW.Cond, TW.76)

SAM.78TW.bfi2 <- cbind(longratings.78.bfi2,longnorm.bfi2,longcrit.78.bfi2, TW.Cond, TW.78)
SAM.79TW.bfi2 <- cbind(longratings.79.bfi2,longnorm.bfi2,longcrit.79.bfi2, TW.Cond, TW.79)

SAM.81TW.bfi2 <- cbind(longratings.81.bfi2,longnorm.bfi2,longcrit.81.bfi2, TW.Cond, TW.81)

SAM.83TW.bfi2 <- cbind(longratings.83.bfi2,longnorm.bfi2,longcrit.83.bfi2, TW.Cond, TW.83)

SAM.86TW.bfi2 <- cbind(longratings.86.bfi2,longnorm.bfi2,longcrit.86.bfi2, TW.Cond, TW.86)

SAM.88TW.bfi2 <- cbind(longratings.88.bfi2,longnorm.bfi2,longcrit.88.bfi2, TW.Cond, TW.88)

SAM.90TW.bfi2 <- cbind(longratings.90.bfi2,longnorm.bfi2,longcrit.90.bfi2, TW.Cond, TW.90)

SAM.92TW.bfi2 <- cbind(longratings.92.bfi2,longnorm.bfi2,longcrit.92.bfi2, TW.Cond, TW.92)

SAM.94TW.bfi2 <- cbind(longratings.94.bfi2,longnorm.bfi2,longcrit.94.bfi2, TW.Cond, TW.94)

SAM.96TW.bfi2 <- cbind(longratings.96.bfi2,longnorm.bfi2,longcrit.96.bfi2, TW.Cond, TW.96)
SAM.97TW.bfi2 <- cbind(longratings.97.bfi2,longnorm.bfi2,longcrit.97.bfi2, TW.Cond, TW.97)

SAM.99TW.bfi2 <- cbind(longratings.99.bfi2,longnorm.bfi2,longcrit.99.bfi2, TW.Cond, TW.99)

SAM.101TW.bfi2 <- cbind(longratings.101.bfi2,longnorm.bfi2,longcrit.101.bfi2, TW.Cond, TW.101)

SAM.104TW.bfi2 <- cbind(longratings.104.bfi2,longnorm.bfi2,longcrit.104.bfi2, TW.Cond, TW.104)

SAM.106TW.bfi2 <- cbind(longratings.106.bfi2,longnorm.bfi2,longcrit.106.bfi2, TW.Cond, TW.106)

SAM.108TW.bfi2 <- cbind(longratings.108.bfi2,longnorm.bfi2,longcrit.108.bfi2, TW.Cond, TW.108)

SAM.110TW.bfi2 <- cbind(longratings.110.bfi2,longnorm.bfi2,longcrit.110.bfi2, TW.Cond, TW.110)

SAM.112TW.bfi2 <- cbind(longratings.112.bfi2,longnorm.bfi2,longcrit.112.bfi2, TW.Cond, TW.112)

SAM.114TW.bfi2 <- cbind(longratings.114.bfi2,longnorm.bfi2,longcrit.114.bfi2, TW.Cond, TW.114)
SAM.115TW.bfi2 <- cbind(longratings.115.bfi2,longnorm.bfi2,longcrit.115.bfi2, TW.Cond, TW.115)

SAM.117TW.bfi2 <- cbind(longratings.117.bfi2,longnorm.bfi2,longcrit.117.bfi2, TW.Cond, TW.117)

SAM.119TW.bfi2 <- cbind(longratings.119.bfi2,longnorm.bfi2,longcrit.119.bfi2, TW.Cond, TW.119)


SAM.122TW.bfi2 <- cbind(longratings.122.bfi2,longnorm.bfi2,longcrit.122.bfi2, TW.Cond, TW.122)

SAM.124TW.bfi2 <- cbind(longratings.124.bfi2,longnorm.bfi2,longcrit.124.bfi2, TW.Cond, TW.124)

SAM.126TW.bfi2 <- cbind(longratings.126.bfi2,longnorm.bfi2,longcrit.126.bfi2, TW.Cond, TW.126)

SAM.128TW.bfi2 <- cbind(longratings.128.bfi2,longnorm.bfi2,longcrit.128.bfi2, TW.Cond, TW.128)

SAM.130TW.bfi2 <- cbind(longratings.130.bfi2,longnorm.bfi2,longcrit.130.bfi2, TW.Cond, TW.130)

SAM.132TW.bfi2 <- cbind(longratings.132.bfi2,longnorm.bfi2,longcrit.132.bfi2, TW.Cond, TW.132)
SAM.133TW.bfi2 <- cbind(longratings.133.bfi2,longnorm.bfi2,longcrit.133.bfi2, TW.Cond, TW.133)

SAM.135TW.bfi2 <- cbind(longratings.135.bfi2,longnorm.bfi2,longcrit.135.bfi2, TW.Cond, TW.135)

SAM.137TW.bfi2 <- cbind(longratings.137.bfi2,longnorm.bfi2,longcrit.137.bfi2, TW.Cond, TW.137)

SAM.140TW.bfi2 <- cbind(longratings.140.bfi2,longnorm.bfi2,longcrit.140.bfi2, TW.Cond, TW.140)

SAM.142TW.bfi2 <- cbind(longratings.142.bfi2,longnorm.bfi2,longcrit.142.bfi2, TW.Cond, TW.142)

SAM.144TW.bfi2 <- cbind(longratings.144.bfi2,longnorm.bfi2,longcrit.144.bfi2, TW.Cond, TW.144)

SAM.146TW.bfi2 <- cbind(longratings.146.bfi2,longnorm.bfi2,longcrit.146.bfi2, TW.Cond, TW.146)

SAM.148TW.bfi2 <- cbind(longratings.148.bfi2,longnorm.bfi2,longcrit.148.bfi2, TW.Cond, TW.148)

SAM.150TW.bfi2 <- cbind(longratings.150.bfi2,longnorm.bfi2,longcrit.150.bfi2, TW.Cond, TW.150)


SAM.bfi.TW <- rbind(SAM.2TW.bfi2, SAM.4TW.bfi2,
                    SAM.6TW.bfi2,
                    SAM.7TW.bfi2,
                    SAM.9TW.bfi2,
                    SAM.11TW.bfi2,
                    SAM.14TW.bfi2,
                    SAM.16TW.bfi2,
                    SAM.18TW.bfi2,
                    SAM.20TW.bfi2,
                    SAM.22TW.bfi2,
                    SAM.24TW.bfi2,
                    SAM.25TW.bfi2,
                    SAM.27TW.bfi2,
                    SAM.29TW.bfi2,
                    SAM.32TW.bfi2,
                    SAM.34TW.bfi2,
                    SAM.36TW.bfi2,
                    SAM.38TW.bfi2,
                    SAM.40TW.bfi2,
                    SAM.42TW.bfi2,
                    SAM.43TW.bfi2,
                    SAM.45TW.bfi2,
                    SAM.47TW.bfi2,
                    SAM.50TW.bfi2,
                    SAM.52TW.bfi2,
                    SAM.54TW.bfi2,
                    SAM.56TW.bfi2,
                    SAM.58TW.bfi2,
                    SAM.60TW.bfi2,
                    SAM.61TW.bfi2,
                    SAM.63TW.bfi2,
                    SAM.65TW.bfi2,
                    SAM.68TW.bfi2,
                    SAM.70TW.bfi2,
                    SAM.72TW.bfi2,
                    SAM.74TW.bfi2,
                    SAM.76TW.bfi2,
                    SAM.78TW.bfi2,
                    SAM.79TW.bfi2,
                    SAM.81TW.bfi2,
                    SAM.83TW.bfi2,
                    SAM.86TW.bfi2,
                    SAM.88TW.bfi2,
                    SAM.90TW.bfi2,
                    SAM.92TW.bfi2,
                    SAM.94TW.bfi2,
                    SAM.96TW.bfi2,
                    SAM.97TW.bfi2,
                    SAM.99TW.bfi2,
                    SAM.101TW.bfi2,
                    SAM.104TW.bfi2,
                    SAM.106TW.bfi2,
                    SAM.108TW.bfi2,
                    SAM.110TW.bfi2,
                    SAM.112TW.bfi2,
                    SAM.114TW.bfi2,
                    SAM.115TW.bfi2,
                    SAM.117TW.bfi2,
                    SAM.119TW.bfi2,
                    SAM.122TW.bfi2,
                    SAM.124TW.bfi2,
                    SAM.126TW.bfi2,
                    SAM.128TW.bfi2,
                    SAM.130TW.bfi2,
                    SAM.132TW.bfi2,
                    SAM.133TW.bfi2,
                    SAM.135TW.bfi2,
                    SAM.137TW.bfi2,
                    SAM.140TW.bfi2,
                    SAM.142TW.bfi2,
                    SAM.144TW.bfi2,
                    SAM.146TW.bfi2, 
                    SAM.148TW.bfi2, 
                    SAM.150TW.bfi2)




# Subtract the Normative profile from the predictors and then mean center predictors for each personality measure

SAM.bfi.TW$Norm.MC <- SAM.bfi.TW$Norm - mean(SAM.bfi.TW$Norm, na.rm=TRUE)

SAM.bfi.TW$SRminusNorm <- SAM.bfi.TW$Sr - SAM.bfi.TW$Norm
SAM.bfi.TW$SR.MC <- SAM.bfi.TW$SRminusNorm - (mean(SAM.bfi.TW$SRminusNorm, na.rm=TRUE))
summary(bfi2.TW <- lmer(Rating ~ 1 + SR.MC + Norm.MC +(1+SR.MC +Norm.MC | JID) +(1 + SR.MC + Norm.MC | TID), data = SAM.bfi.TW, control = lmerControl(optimizer = 'bobyqa')))
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: Rating ~ 1 + SR.MC + Norm.MC + (1 + SR.MC + Norm.MC | JID) +  
##     (1 + SR.MC + Norm.MC | TID)
##    Data: SAM.bfi.TW
## Control: lmerControl(optimizer = "bobyqa")
## 
## REML criterion at convergence: 70355.7
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.1753 -0.7546  0.0183  0.7459  3.3081 
## 
## Random effects:
##  Groups   Name        Variance Std.Dev. Corr       
##  JID      (Intercept) 0.096303 0.31033             
##           SR.MC       0.012250 0.11068  -0.32      
##           Norm.MC     0.137508 0.37082  -0.19  0.71
##  TID      (Intercept) 0.001404 0.03747             
##           SR.MC       0.013593 0.11659  0.59       
##           Norm.MC     0.048579 0.22041  0.78  0.63 
##  Residual             1.098907 1.04829             
## Number of obs: 23508, groups:  JID, 271; TID, 75
## 
## Fixed effects:
##              Estimate Std. Error        df t value Pr(>|t|)    
## (Intercept)   3.21080    0.02060 263.61053 155.839  < 2e-16 ***
## SR.MC         0.08557    0.01623  95.18750   5.273 8.39e-07 ***
## Norm.MC       0.26845    0.03573 142.84646   7.513 5.77e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##         (Intr) SR.MC 
## SR.MC   -0.017       
## Norm.MC  0.006  0.579
names(SAM.bfi.TW)
##  [1] "JID"         "SNSFreq"     "TwitFreq"    "IGFreq"      "Rating.item"
##  [6] "Rating"      "Norm.item"   "Norm"        "TID"         "Sr.item"    
## [11] "Sr"          "ï..ConType"  "ï..Original" "Norm.MC"     "SRminusNorm"
## [16] "SR.MC"
names(SAM.bfi.TW)[names(SAM.bfi.TW) == "ï..Original"] <- "Original"

#original posts
summary(twitter.original <- lmer(Rating ~ 1 + SR.MC*Original + Norm.MC*Original  +(1 + SR.MC + Norm.MC | JID) + (1 + SR.MC + Norm.MC | TID), data = SAM.bfi.TW, control = lmerControl(optimizer ='optimx', optCtrl=list(method='nlminb'))))
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: 
## Rating ~ 1 + SR.MC * Original + Norm.MC * Original + (1 + SR.MC +  
##     Norm.MC | JID) + (1 + SR.MC + Norm.MC | TID)
##    Data: SAM.bfi.TW
## Control: 
## lmerControl(optimizer = "optimx", optCtrl = list(method = "nlminb"))
## 
## REML criterion at convergence: 70378
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.1714 -0.7557  0.0183  0.7465  3.3089 
## 
## Random effects:
##  Groups   Name        Variance Std.Dev. Corr       
##  JID      (Intercept) 0.095713 0.30938             
##           SR.MC       0.012232 0.11060  -0.32      
##           Norm.MC     0.137609 0.37096  -0.19  0.71
##  TID      (Intercept) 0.001219 0.03491             
##           SR.MC       0.013795 0.11745  0.64       
##           Norm.MC     0.049218 0.22185  0.79  0.63 
##  Residual             1.098927 1.04830             
## Number of obs: 23508, groups:  JID, 271; TID, 75
## 
## Fixed effects:
##                    Estimate Std. Error         df t value Pr(>|t|)    
## (Intercept)       3.232e+00  2.305e-02  2.366e+02 140.212  < 2e-16 ***
## SR.MC             8.342e-02  2.323e-02  8.128e+01   3.591 0.000562 ***
## Original         -5.943e-03  2.992e-03  5.100e+01  -1.987 0.052356 .  
## Norm.MC           2.849e-01  4.867e-02  1.025e+02   5.853 5.84e-08 ***
## SR.MC:Original    5.197e-04  4.715e-03  6.962e+01   0.110 0.912540    
## Original:Norm.MC -4.719e-03  9.401e-03  7.097e+01  -0.502 0.617251    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) SR.MC  Orignl Nrm.MC SR.MC:
## SR.MC        0.066                            
## Original    -0.458 -0.173                     
## Norm.MC      0.093  0.573 -0.205              
## SR.MC:Orgnl -0.109 -0.712  0.248 -0.384       
## Orgnl:Nr.MC -0.137 -0.403  0.304 -0.676  0.568