《共分散構造分析》
#パッケージインストール
install.packages(c("lavaan","semTools","semPlot"))

#探索的因子分析
X <- read.csv("sem.csv", header=T)
factanal(X,factors=13,rotation="promax")


#確認的因子分析
K1.model.cfa <- '
f1 =~ sannka1 + sannka2 + sannka3
f2 =~ meisei1 + meisei2 
f3 =~ jiko1 + jiko2
f4 =~ syuui1
f5 =~ shinnki1 + shinnki2 + shinnki3
f6 =~ keikenn1 + keikenn2 + keikenn3 + nouryoku1 + nouryoku3
f7 =~ kannshinn1 + kannshinn2
f8 =~ timeido1 + timeido2 + timeido3
f9 =~ pops1 + pops2 + pops3
f10 =~ anime2
f11 =~ ennka3
f12 =~ jazz1 + jazz2 + jazz3
f13 =~ ninnzuu3'

library(lavaan)
fit <-cfa(model=K1.model.cfa,data=X,estimator="ML")
summary(fit, fit.measures=TRUE,standardized=T,rsquare=T)


#共分散構造分析
K1.model.sem<-'
f1 =~ sannka1 + sannka2 + sannka3
f2 =~ meisei1 + meisei2 
f3 =~ jiko1 + jiko2
f4 =~ syuui1
f5 =~ shinnki1 + shinnki2 + shinnki3
f6 =~ keikenn1 + keikenn2 + keikenn3 + nouryoku1 + nouryoku3
f7 =~ kannshinn1 + kannshinn2
f8 =~ timeido1 + timeido2 + timeido3
f9 =~ pops1 + pops2 + pops3
f10 =~ anime2
f11 =~ ennka3
f12 =~ jazz1 + jazz2 + jazz3
f13 =~ ninnzuu3
f1~f2
f1~f3
f1~f4
f1~f5
f1~f6
f1~f7
f1~f8
f1~f9
f1~f10
f1~f11
f1~f12
f1~f13'

fit <-sem(model=K1.model.sem, data=X, estimator="ML")
summary(object=fit)
summary(fit, fit.meas)


《コンジョイント分析》
#パッケージインストール
install.packages("conjoint")


《参加意図①》
library(conjoint)

experiment<-expand.grid(
Base=c("日本","海外"), 
Cheer=c("応援している","応援していない"), 
Post=c("投稿している","投稿していない"))

design <-
caFactorialDesign(data=experiment,type="orthogonal")
design

X <- read.csv("sannka1.csv", header=T)
tprefm<- X
tprof<-caEncodedDesign(design)
tlevn<-as.matrix(c("日本","海外","応援している","応援していない","投稿している","投稿していない"))
Conjoint(tprefm,tprof,tlevn)


《参加意図②》
library(conjoint)

experiment<-expand.grid(
Base=c("明確である","明確でない"), 
Cheer=c("フィードバックあり","フィードバックなし"), 
Post=c("経済的報酬","体験的報酬","報酬なし"))

design <-
caFactorialDesign(data=experiment,type="orthogonal")
design

X <- read.csv("sannka2.csv", header=T)
tprefm<- X
tprof<-caEncodedDesign(design)
tlevn<-as.matrix(c("明確である","明確でない","フィードバックあり","フィードバックなし","経済的報酬","体験的報酬","報酬なし"))
Conjoint(tprefm,tprof,tlevn)


《参加意図③》
library(conjoint)

experiment<-expand.grid(
Base=c("定期開催あり","定期開催なし"), 
Cheer=c("モバイル端末から参加可能","モバイル端末から参加不可能"), 
Post=c("選択","ワード","文"))

design <-
caFactorialDesign(data=experiment,type="orthogonal")
design

X <- read.csv("sannka3.csv", header=T)
tprefm<- X
tprof<-caEncodedDesign(design)
tlevn<-as.matrix(c("定期開催あり","定期開催なし","モバイル端末から参加可能","モバイル端末から参加不可能","選択","ワード","文"))
Conjoint(tprefm,tprof,tlevn)


《視聴意図①》
library(conjoint)

experiment<-expand.grid(
Base=c("日本","海外"), 
Cheer=c("応援している","応援していない"), 
Post=c("投稿している","投稿していない"))

design <-
caFactorialDesign(data=experiment,type="orthogonal")
design

X <- read.csv("shityou1.csv", header=T)
tprefm<- X
tprof<-caEncodedDesign(design)
tlevn<-as.matrix(c("日本","海外","応援している","応援していない","投稿している","投稿していない"))
Conjoint(tprefm,tprof,tlevn)

《視聴意図②》
library(conjoint)

experiment<-expand.grid(
Base=c("明確である","明確でない"), 
Cheer=c("フィードバックあり","フィードバックなし"), 
Post=c("経済的報酬","体験的報酬","報酬なし"))

design <-
caFactorialDesign(data=experiment,type="orthogonal")
design

X <- read.csv("shityou2.csv", header=T)
tprefm<- X
tprof<-caEncodedDesign(design)
tlevn<-as.matrix(c("明確である","明確でない","フィードバックあり","フィードバックなし","経済的報酬","体験的報酬","報酬なし"))
Conjoint(tprefm,tprof,tlevn)


《視聴意図③》
library(conjoint)

experiment<-expand.grid(
Base=c("定期開催あり","定期開催なし"), 
Cheer=c("モバイル端末から参加可能","モバイル端末から参加不可能"), 
Post=c("選択","ワード","文"))

design <-
caFactorialDesign(data=experiment,type="orthogonal")
design

X <- read.csv("shityou3.csv", header=T)
tprefm<- X
tprof<-caEncodedDesign(design)
tlevn<-as.matrix(c("定期開催あり","定期開催なし","モバイル端末から参加可能","モバイル端末から参加不可能","選択","ワード","文"))
Conjoint(tprefm,tprof,tlevn)


《アーティスト好意度①》
library(conjoint)

experiment<-expand.grid(
Base=c("日本","海外"), 
Cheer=c("応援している","応援していない"), 
Post=c("投稿している","投稿していない"))

design <-
caFactorialDesign(data=experiment,type="orthogonal")
design

X <- read.csv("shityou1.csv", header=T)
tprefm<- X
tprof<-caEncodedDesign(design)
tlevn<-as.matrix(c("日本","海外","応援している","応援していない","投稿している","投稿していない"))
Conjoint(tprefm,tprof,tlevn)

《アーティスト好意度②》
library(conjoint)

experiment<-expand.grid(
Base=c("明確である","明確でない"), 
Cheer=c("フィードバックあり","フィードバックなし"), 
Post=c("経済的報酬","体験的報酬","報酬なし"))

design <-
caFactorialDesign(data=experiment,type="orthogonal")
design

X <- read.csv("kouido2.csv", header=T)
tprefm<- X
tprof<-caEncodedDesign(design)
tlevn<-as.matrix(c("明確である","明確でない","フィードバックあり","フィードバックなし","経済的報酬","体験的報酬","報酬なし"))
Conjoint(tprefm,tprof,tlevn)

《アーティスト好意度③》
library(conjoint)

experiment<-expand.grid(
Base=c("定期開催あり","定期開催なし"), 
Cheer=c("モバイル端末から参加可能","モバイル端末から参加不可能"), 
Post=c("選択","ワード","文"))

design <-
caFactorialDesign(data=experiment,type="orthogonal")
design

X <- read.csv("kouido3.csv", header=T)
tprefm<- X
tprof<-caEncodedDesign(design)
tlevn<-as.matrix(c("定期開催あり","定期開催なし","モバイル端末から参加可能","モバイル端末から参加不可能","選択","ワード","文"))
Conjoint(tprefm,tprof,tlevn)