mtrxnames <- c("L1" , "L2") akyuz2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/akyuz2.csv", header=FALSE)) rownames(akyuz2) <-mtrxnames colnames(akyuz2) <-mtrxnames akkul2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/akkul2.csv", header=FALSE)) rownames(akkul2) <-mtrxnames colnames(akkul2) <-mtrxnames altundag2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/altundag2.csv", header=FALSE)) rownames(altundag2) <-mtrxnames colnames(altundag2) <-mtrxnames arcan2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/arcan2.csv", header=FALSE)) rownames(arcan2) <-mtrxnames colnames(arcan2) <-mtrxnames armagan2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/armagan2.csv", header=FALSE)) rownames(armagan2) <-mtrxnames colnames(armagan2) <-mtrxnames arslantas2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/arslantas2.csv", header=FALSE)) rownames(arslantas2) <-mtrxnames colnames(arslantas2) <-mtrxnames atasayar2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/atasayar2.csv", header=FALSE)) rownames(atasayar2) <-mtrxnames colnames(atasayar2) <-mtrxnames aydemir2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/aydemir2.csv", header=FALSE)) rownames(aydemir2) <-mtrxnames colnames(aydemir2) <-mtrxnames bahadir2 <-as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/bahadir2.csv", header=FALSE)) rownames(bahadir2) <-mtrxnames colnames(bahadir2) <-mtrxnames baltaci2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/baltaci2.csv", header=FALSE)) rownames(baltaci2) <-mtrxnames colnames(baltaci2) <-mtrxnames baran2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/baran2.csv", header=FALSE)) rownames(baran2) <-mtrxnames colnames(baran2) <-mtrxnames bilge2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/bilge2.csv", header=FALSE)) rownames(bilge2) <-mtrxnames colnames(bilge2) <-mtrxnames bilgili2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/bilgili2.csv", header=FALSE)) rownames(bilgili2) <-mtrxnames colnames(bilgili2) <-mtrxnames budakoglu2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/budakoglu2.csv", header=FALSE)) rownames(budakoglu2) <-mtrxnames colnames(budakoglu2) <-mtrxnames cakir2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/cakir2.csv", header=FALSE)) rownames(cakir2) <-mtrxnames colnames(cakir2) <-mtrxnames certel2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/certel2.csv", header=FALSE)) rownames(certel2) <-mtrxnames colnames(certel2) <-mtrxnames ciraci2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/ciraci2.csv", header=FALSE)) rownames(ciraci2) <-mtrxnames colnames(ciraci2) <-mtrxnames civitci2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/civitci2.csv", header=FALSE)) rownames(civitci2) <-mtrxnames colnames(civitci2) <-mtrxnames civitci22 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/civitci22.csv", header=FALSE)) rownames(civitci22) <-mtrxnames colnames(civitci22) <-mtrxnames ciyin2 <-as.matrix (read.csv("C:/Users/vilda/Desktop/higherorder/ciyin2.csv", header=FALSE)) rownames(ciyin2) <-mtrxnames colnames(ciyin2) <-mtrxnames durak2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/durak2.csv", header=FALSE)) rownames(durak2) <-mtrxnames colnames(durak2) <-mtrxnames duru2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/duru2.csv", header=FALSE)) rownames(duru2) <-mtrxnames colnames(duru2) <-mtrxnames ekinci2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/ekinci2.csv", header=FALSE)) rownames(ekinci2) <-mtrxnames colnames(ekinci2) <-mtrxnames erdem2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/erdem2.csv", header=FALSE)) rownames(erdem2) <-mtrxnames colnames(erdem2) <-mtrxnames fikri2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/fikri2.csv", header=FALSE)) rownames(fikri2) <-mtrxnames colnames(fikri2) <-mtrxnames gamsiz2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/gamsiz2.csv", header=FALSE)) rownames(gamsiz2) <-mtrxnames colnames(gamsiz2) <-mtrxnames genc2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/genc2.csv", header=FALSE)) rownames(genc2) <-mtrxnames colnames(genc2) <-mtrxnames golcu2 <-as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/golcu2.csv", header=FALSE)) rownames(golcu2) <-mtrxnames colnames(golcu2) <-mtrxnames guven2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/guven2.csv", header=FALSE)) rownames(guven2) <-mtrxnames colnames(guven2) <-mtrxnames kapkiran2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/kapkiran2.csv", header=FALSE)) rownames(kapkiran2) <-mtrxnames colnames(kapkiran2) <-mtrxnames kapkiran22 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/kapkiran22.csv", header=FALSE)) rownames(kapkiran22) <-mtrxnames colnames(kapkiran22) <-mtrxnames karababa2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/karababa2.csv", header=FALSE)) rownames(karababa2) <-mtrxnames colnames(karababa2) <-mtrxnames karaoglu2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/karaoglu2.csv", header=FALSE)) rownames(karaoglu2) <-mtrxnames colnames(karaoglu2) <-mtrxnames meralbahar2<- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/meralbahar2.csv", header=FALSE)) rownames(meralbahar2) <-mtrxnames colnames(meralbahar2) <-mtrxnames mercan2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/mercan2.csv", header=FALSE)) rownames(mercan2) <-mtrxnames colnames(mercan2) <-mtrxnames mert2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/mert2.csv", header=FALSE)) rownames(mert2) <-mtrxnames colnames(mert2) <-mtrxnames meryem2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/meryem2.csv", header=FALSE)) rownames(meryem2) <-mtrxnames colnames(meryem2) <-mtrxnames nurtenkaya2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/nurtenkaya2.csv", header=FALSE)) rownames(nurtenkaya2) <-mtrxnames colnames(nurtenkaya2) <-mtrxnames oner2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/oner2.csv", header=FALSE)) rownames(oner2) <-mtrxnames colnames(oner2) <-mtrxnames pamukcu2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/pamukcu2.csv", header=FALSE)) rownames(pamukcu2) <-mtrxnames colnames(pamukcu2) <-mtrxnames saricam2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/saricam2.csv", header=FALSE)) rownames(saricam2) <-mtrxnames colnames(saricam2) <-mtrxnames selma2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/selma2.csv", header=FALSE)) rownames(selma2) <-mtrxnames colnames(selma2) <-mtrxnames sinem2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/sinem2.csv", header=FALSE)) rownames(sinem2) <-mtrxnames colnames(sinem2) <-mtrxnames tumkaya2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/tumkaya2.csv", header=FALSE)) rownames(tumkaya2) <-mtrxnames colnames(tumkaya2) <-mtrxnames tutgun2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/tutgun2.csv", header=FALSE)) rownames(tutgun2) <-mtrxnames colnames(tutgun2) <-mtrxnames tutgun22 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/tutgun22.csv", header=FALSE)) rownames(tutgun22) <-mtrxnames colnames(tutgun22) <-mtrxnames ucar2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/ucar2.csv", header=FALSE)) rownames(ucar2) <-mtrxnames colnames(ucar2) <-mtrxnames uluyalcinkaya2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/uluyalcinkaya2.csv", header=FALSE)) rownames(uluyalcinkaya2) <-mtrxnames colnames(uluyalcinkaya2) <-mtrxnames ummet2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/ummet2.csv", header=FALSE)) rownames(ummet2) <-mtrxnames colnames(ummet2) <-mtrxnames yazicilar2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/yazicilar2.csv", header=FALSE)) rownames(yazicilar2) <-mtrxnames colnames(yazicilar2) <-mtrxnames kostak2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/kostak2.csv", header=FALSE)) rownames(kostak2) <-mtrxnames colnames(kostak2) <-mtrxnames yakut2 <- as.matrix(read.csv("C:/Users/vilda/Desktop/higherorder/yakut2.csv", header=FALSE)) rownames(yakut2) <-mtrxnames colnames(yakut2) <-mtrxnames #####52 çalışma ################3 dataucla=list(akkul2, akyuz2 , altundag2 , arcan2 , armagan2 , arslantas2 , atasayar2 , aydemir2 , bahadir2 , baltaci2 , baran2 , bilge2 , bilgili2 , budakoglu2 , cakir2, certel2 , ciraci2 , civitci2 , civitci22 , ciyin2 , durak2 , duru2 , ekinci2 , erdem2 , fikri2 , gamsiz2 , genc2 , golcu2 , guven2 , kapkiran2 , kapkiran22 , karababa2 , karaoglu2 , kostak2 , meralbahar2 , mercan2 , mert2 , meryem2 , nurtenkaya2 , oner2 , pamukcu2 , saricam2 , selma2, sinem2 , tumkaya2 , tutgun2 , tutgun22 , ucar2, uluyalcinkaya2 , ummet2 , yakut2, yazicilar2) dataucla nucla <- c(117, 464, 144, 381, 93, 669, 300, 242, 161, 873, 426, 200, 428, 225, 556, 144, 55, 439, 398, 267, 786, 190, 392, 204, 665, 325, 202, 411, 182, 336, 110, 683, 446, 881, 630, 343, 281, 160, 667, 657, 278, 326, 144, 942, 422, 686, 595, 888, 1451, 117, 429, 319) is.pd(dataucla) dataucla[!is.pd(dataucla)] #TSSEM 1. adım analizleri fixed effect model için fixed1 <- tssem1(dataucla, nucla, method="FEM") summary(fixed1) coef(fixed1) #TSSEM 2. adım analizleri #kestirimlerden önce RAM yaklasımı ile matrisleri hesaplarız#### ## Define the S matrix that includes both Phi and Psi ## Phi matrix: 1x1 correlation matrix between the latent factors ( Phi <- matrix(c(1), nrow=1, ncol=1) ) ## Psi matrix: 2x2 diagonal matrix of the error variances ( Psi <- Diag(c(".2*e1",".2*e2")) ) ## Create a block diagonal matrix as the S matrix S1 <- bdiagMat(list(Psi, Phi)) #BU ADIM GEREKLİ OLMASADA YORUMLAMA KOLAYLIĞI AÇISINDAN ÖNEMLİ dimnames(S1)[[1]] <- c("L1","L2", "L") dimnames(S1)[[2]] <- dimnames(S1)[[1]] S1 ## Convert it into a MxMatrix class S1 <- as.mxMatrix(S1) ## Lambda matrix: 2x1 factor loadings ## Arrange the data by row ( Lambda <- matrix(c(".3*Loading1", ".3*Loading2"), nrow=2, ncol=1, byrow=TRUE) ) ## Create a 2x2 of zeros ( Zero2x2 <- matrix(0, nrow=2, ncol=2) ) ## Create a 1x3 of zeros ( Zero1x3 <- matrix(0, nrow=1, ncol=3) ) ## Define the A matrix A1 <- rbind(cbind(Zero2x2, Lambda), Zero1x3) ## bu adım yine modeli incelemek için kullanışlı dimnames(A1)[[1]] <- c("L1","L2", "L") dimnames(A1)[[2]] <- dimnames(A1)[[1]] A1 ## Convert it into a MxMatrix class A1 <- as.mxMatrix(A1) ## F matrix to select the observed variables ## First 20 elements are observed variables ## Last 1 elements are latent variables F1 <- create.Fmatrix(c(1,1,0), name="F1", as.mxMatrix=FALSE) ## bu adım yine modeli incelemek için kullanışlı dimnames(F1)[[1]] <- c("L1","L2") dimnames(F1)[[2]] <- c("L1","L2", "L") ## Convert it into a MxMatrix class F1 <- as.mxMatrix(F1) ##matrisler hesaplandıkan sonra TSSEM ikinci asamama kestirimlerine gecilir### ##tssem ikinci asaması icin fixed1 den ögeleri alarak rem ya da seme göre ###birleÅŸtirilmiÅŸ kor matrisi ve asmptotk örneklem kor mat oluÅŸturur. fixed2 <- tssem2(fixed1, Amatrix=A1, Smatrix=S1, Fmatrix=F1, diag.constraints=FALSE) summary(fixed2) coef(fixed2) plot(fixed2)