ich versuche eine Präferenzprognose mittels Choice-Based-Conjoint-Analyse mit R hinzukriegen.
Ich habe den Leitfaden aus von Gansser dafür genutzt, der das hier vorsieht:
Code: Alles auswählen
library(AlgDesign)
vd<- gen.factorial( c(4,3,3,3), varNames=c("Ort", "SL", "ML", "Unt"), factors="all")
set.seed(54321)
fd<- optFederov(~.,vd,16)
alt1 <- fd$design
alt2 <- alt1
alt1 <- transform(alt1,z1=runif(16))
alt2 <- transform(alt2,z2=runif(16))
alt1_sort <- alt1[order(alt1$z1),]
alt2_sort <- alt2[order(alt2$z2),]
alt1_sort
alt2_sort
library(survival)
clogout1 <- clogit(E~K+Ort+SL+ML+Unt+strata(APN), data=Datenmatrix)
clogout1
clogout1$loglik
Beim letzten Schritt, also
Code: Alles auswählen
library(survival)
clogout1 <- clogit(E~K+Ort+SL+ML+Unt+strata(APN), data=Datenmatrix)
clogout1
clogout1$loglik
Wie komme ich an die gesamte Tabelle?
Danke und schöne Grüße
Tom
Edit: Die Namen der Variablen habe ich natürlich angepasst.