new file to run the script and more designs for drive

SE_DRIVE/Designs/bayeffdesignconstr.ngd

+Design	Choice situation	alt1.x1	alt1.x2	alt1.x3	alt2.x1	alt2.x2	alt2.x3	Block	
+1	1	20	200	0	40	25	0	7	
+1	2	40	100	0	80	50	25	7	
+1	3	80	100	50	40	100	100	10	
+1	4	20	25	50	40	50	25	6	
+1	5	60	50	50	20	100	50	8	
+1	6	60	25	25	40	200	0	7	
+1	7	80	25	100	60	200	100	1	
+1	8	80	50	0	20	25	100	3	
+1	9	60	100	0	80	25	50	5	
+1	10	60	100	50	20	50	100	9	
+1	11	60	200	50	80	50	25	3	
+1	12	80	50	25	20	50	50	2	
+1	13	20	200	25	60	25	0	8	
+1	14	20	100	50	80	25	100	10	
+1	15	40	25	100	60	100	0	3	
+1	16	80	25	25	40	25	50	2	
+1	17	20	100	50	40	25	0	4	
+1	18	20	100	100	60	50	25	4	
+1	19	20	25	50	60	25	0	1	
+1	20	80	50	50	80	100	0	4	
+1	21	80	25	0	40	200	100	8	
+1	22	40	200	0	20	100	100	7	
+1	23	40	100	50	20	50	100	2	
+1	24	20	200	50	80	50	25	7	
+1	25	80	100	0	40	200	100	7	
+1	26	20	50	25	40	50	0	6	
+1	27	40	50	50	20	100	25	4	
+1	28	40	50	25	60	25	0	5	
+1	29	80	50	0	20	200	0	10	
+1	30	20	100	50	80	50	100	1	
+1	31	60	25	50	20	200	50	8	
+1	32	40	200	25	80	25	0	1	
+1	33	40	200	0	60	100	0	5	
+1	34	80	50	0	60	25	100	8	
+1	35	60	25	0	20	100	25	3	
+1	36	20	100	100	40	25	25	6	
+1	37	60	200	100	80	25	100	3	
+1	38	40	25	25	20	100	25	2	
+1	39	40	200	0	80	100	100	1	
+1	40	20	100	50	80	25	50	2	
+1	41	40	100	25	80	25	50	10	
+1	42	60	100	0	20	200	50	10	
+1	43	60	200	0	80	100	25	9	
+1	44	80	50	0	40	25	25	8	
+1	45	20	100	100	80	100	50	4	
+1	46	60	200	50	80	50	100	5	
+1	47	60	50	0	60	25	25	5	
+1	48	60	50	25	20	50	100	1	
+1	49	60	100	50	20	200	25	8	
+1	50	60	200	100	80	200	25	4	
+1	51	40	50	0	20	200	100	2	
+1	52	60	25	100	80	50	0	3	
+1	53	40	25	50	60	50	25	9	
+1	54	40	200	50	60	25	0	10	
+1	55	40	50	100	80	25	100	5	
+1	56	40	100	50	60	25	50	9	
+1	57	80	50	25	40	50	50	2	
+1	58	80	200	0	60	200	100	5	
+1	59	80	50	0	20	100	50	1	
+1	60	20	100	100	60	25	100	10	
+1	61	60	25	25	20	100	0	2	
+1	62	80	25	0	20	25	25	6	
+1	63	60	50	100	20	200	0	2	
+1	64	20	200	0	60	25	50	7	
+1	65	20	200	50	80	100	100	5	
+1	66	80	25	25	60	25	100	9	
+1	67	60	50	25	20	200	50	8	
+1	68	20	100	100	40	25	0	7	
+1	69	40	25	100	80	100	0	9	
+1	70	20	25	100	40	25	25	6	
+1	71	80	25	50	40	200	25	6	
+1	72	40	50	50	20	100	50	1	
+1	73	80	50	0	40	200	50	5	
+1	74	20	100	25	40	50	0	6	
+1	75	80	25	25	40	200	50	9	
+1	76	80	50	25	20	200	25	10	
+1	77	40	200	0	60	50	100	4	
+1	78	80	25	25	20	100	0	10	
+1	79	20	200	100	80	200	50	7	
+1	80	60	50	100	40	200	0	1	
+1	81	40	100	50	80	50	25	9	
+1	82	40	100	50	60	50	50	9	
+1	83	20	100	50	80	50	25	4	
+1	84	60	25	0	20	200	100	3	
+1	85	20	25	100	80	25	0	4	
+1	86	60	200	25	80	100	100	3	
+1	87	60	100	50	80	25	100	7	
+1	88	60	25	100	40	200	100	8	
+1	89	60	200	100	80	200	50	6	
+1	90	40	50	25	20	100	50	1	
+1	91	80	200	0	60	100	100	6	
+1	92	80	50	50	40	100	25	6	
+1	93	80	25	25	60	200	25	4	
+1	94	60	200	0	80	50	50	6	
+1	95	20	50	25	60	50	0	3	
+1	96	40	25	25	80	25	0	10	
+1	97	40	100	0	20	100	25	10	
+1	98	80	100	0	40	200	25	1	
+1	99	60	100	100	80	25	25	5	
+1	100	40	25	100	20	200	50	3	
+1	101	20	200	50	60	50	50	8	
+1	102	40	200	100	60	50	0	3	
+1	103	40	50	25	20	100	0	2	
+1	104	20	200	25	40	50	0	7	
+1	105	80	25	25	60	100	0	2	
+1	106	40	50	100	40	100	50	9	
+1	107	20	200	50	60	50	0	4	
+1	108	40	100	25	20	200	100	7	
+1	109	80	25	0	60	200	50	6	
+1	110	20	200	0	80	25	0	3	
+1	111	80	50	0	60	50	25	8	
+1	112	60	100	25	40	200	0	4	
+1	113	80	25	100	20	50	100	5	
+1	114	40	200	0	60	50	50	5	
+1	115	40	200	100	80	100	25	2	
+1	116	20	50	100	40	100	25	8	
+1	117	60	100	25	20	200	50	1	
+1	118	20	100	100	40	100	0	9	
+1	119	40	200	100	60	100	50	9	
+1	120	80	25	100	40	200	0	10	
+||||||||||
+design
+;alts = alt1, alt2
+;block = 10
+
+;eff = (mnl,d,mean)
+;rep = 1000
+;bdraws = halton(1000)
+;bseed = 2333344
+;rseed = 2333344
+;alg = swap
+;rows = 120
+
+
+
+;cond:
+if(alt1.x2+alt1.x3 >alt2.x2+alt2.x3   , alt1.x1<alt2.x1),
+if(alt1.x2+alt1.x3 <alt2.x2+alt2.x3   , alt1.x1>alt2.x1),
+if(alt1.x2=alt2.x2  , alt1.x3<>alt2.x3 )
+
+;model:
+U(alt1) = b1[(n,-0.036,0.02)]   * x1[20,40,60,80]
+        + b2[(n,-0.0034,0.002)] * x2[25,50,100,200]
+        + b3[(n,-0.0049,0.003)] * x3[0,25,50,100]
+/
+U(alt2) = b1 * x1
+        + b2 * x2
+        + b3 * x3
+$
\ No newline at end of file

SE_DRIVE/Designs/effconstrsmall.ngd

+Design	Choice situation	alt1.x1	alt1.x2	alt1.x3	alt2.x1	alt2.x2	alt2.x3	Block	
+1	1	20	100	25	60	100	0	3	
+1	2	80	25	25	40	200	100	5	
+1	3	60	200	50	80	50	25	2	
+1	4	60	25	25	40	50	50	3	
+1	5	20	25	25	40	25	0	2	
+1	6	40	100	100	80	100	25	5	
+1	7	40	200	50	80	200	25	3	
+1	8	60	25	50	20	200	50	3	
+1	9	80	25	100	40	200	50	3	
+1	10	20	50	50	80	50	0	5	
+1	11	60	50	0	20	200	100	2	
+1	12	60	25	0	20	200	100	1	
+1	13	80	100	50	60	200	0	5	
+1	14	40	50	100	60	100	0	3	
+1	15	80	100	100	60	200	0	5	
+1	16	20	100	50	40	50	0	1	
+1	17	20	200	0	80	100	100	1	
+1	18	60	200	25	80	25	50	5	
+1	19	20	200	50	80	50	100	3	
+1	20	40	100	0	20	100	50	2	
+1	21	60	100	25	80	50	25	4	
+1	22	60	100	50	80	50	25	2	
+1	23	80	100	0	40	100	100	2	
+1	24	80	50	50	60	200	100	5	
+1	25	60	50	100	20	200	50	1	
+1	26	80	50	25	40	100	0	4	
+1	27	80	50	100	40	100	50	2	
+1	28	20	25	100	80	25	50	5	
+1	29	20	50	100	80	50	0	1	
+1	30	20	200	50	60	100	50	4	
+1	31	20	200	50	40	100	50	3	
+1	32	40	100	25	80	25	50	1	
+1	33	80	50	100	60	200	25	1	
+1	34	40	50	100	80	25	50	4	
+1	35	40	25	100	20	200	25	3	
+1	36	20	200	25	60	100	0	2	
+1	37	60	50	50	20	200	100	4	
+1	38	80	50	25	60	25	50	5	
+1	39	60	100	0	20	200	0	5	
+1	40	40	200	100	80	25	25	2	
+1	41	20	200	0	60	25	100	5	
+1	42	40	200	50	20	200	100	4	
+1	43	20	50	0	60	25	0	4	
+1	44	20	25	50	80	25	25	4	
+1	45	60	25	25	20	100	100	2	
+1	46	80	25	25	40	200	25	1	
+1	47	60	50	0	40	100	100	1	
+1	48	80	25	25	40	50	25	2	
+1	49	60	25	25	20	50	25	3	
+1	50	40	200	0	60	50	25	4	
+1	51	40	100	0	20	25	100	4	
+1	52	80	50	100	40	100	100	2	
+1	53	20	100	50	40	100	0	5	
+1	54	40	100	0	20	50	50	4	
+1	55	20	50	100	80	100	0	1	
+1	56	40	200	100	60	25	0	3	
+1	57	60	200	25	80	25	0	1	
+1	58	60	25	0	20	50	100	3	
+1	59	60	200	0	80	25	50	1	
+1	60	80	50	50	20	25	100	4	
+||||||||||
+design
+;alts = alt1, alt2
+;block = 5
+;eff = (mnl, d)
+;alg = swap
+;rows = 60
+
+
+
+;cond:
+if(alt1.x2+alt1.x3 >alt2.x2+alt2.x3   , alt1.x1<alt2.x1),
+if(alt1.x2+alt1.x3 <alt2.x2+alt2.x3   , alt1.x1>alt2.x1),
+if(alt1.x2=alt2.x2  , alt1.x3<>alt2.x3 )
+
+;model:
+U(alt1) = b1[-0.036]  * x1[20,40,60,80]
+        + b2[-0.0034] * x2[25,50,100,200]
+        + b3[-0.0049] * x3[0,25,50,100]
+/
+U(alt2) = b1 * x1
+        + b2 * x2
+        + b3 * x3
+$
\ No newline at end of file

generatemd.R

+rmarkdown::render("simulation_output.rmd",
+                  output_file = "SE_DRIVE/simulation_DRIVE_wrongpriors.html")
\ No newline at end of file

parameters_SE Design-Agri.R

-
-
-designpath<- "designs-seagri/"
-
-resps =360  # number of respondents
-nosim=2 # number of simulations to run (about 500 is minimum)
-
-#betacoefficients should not include "-"
-basc = 4.2 ## very high asc
-bprof = 0.3
-bexp = 0.3
-bdomestic = 0.3
-bforeign = 0.3
-bdamage = 0.6
-bprice = 0.2
-
-manipulations = list(alt1.professional=     expr(alt1.initiator==1), 
-                     alt2.professional=     expr(alt2.initiator==1),
-                     alt1.expert      =     expr(alt1.initiator==2), 
-                     alt2.expert      =     expr(alt2.initiator==2),
-                     alt1.domestic    =     expr(alt1.funding==1), 
-                     alt2.domestic    =     expr(alt2.funding==1),
-                     alt1.foreign     =     expr(alt1.funding==2), 
-                     alt2.foreign     =     expr(alt2.funding==2))
-
-
-
-
-
-
-#place your utility functions here
-u<-list(
-  v1 =V.1 ~  bprof*alt1.professional+ bexp * alt1.expert + bdomestic * alt1.domestic + bforeign * alt1.foreign + bdamage*alt1.damage + bprice * alt1.compensation,
-  v2 =V.2 ~  bprof*alt2.professional + bexp * alt2.expert + bdomestic * alt2.domestic + bforeign * alt2.foreign + bdamage*alt2.damage + bprice * alt2.compensation,
-  v3 =V.3 ~ basc)
-
-source("simulationcore_purrr.R")
+
+
+designpath<- "designs-seagri/"
+
+resps =360  # number of respondents
+nosim=2 # number of simulations to run (about 500 is minimum)
+
+#betacoefficients should not include "-"
+basc = 4.2 ## very high asc
+bprof = 0.3
+bexp = 0.3
+bdomestic = 0.3
+bforeign = 0.3
+bdamage = 0.6
+bprice = 0.2
+
+manipulations = list(alt1.professional=     expr(alt1.initiator==1), 
+                     alt2.professional=     expr(alt2.initiator==1),
+                     alt1.expert      =     expr(alt1.initiator==2), 
+                     alt2.expert      =     expr(alt2.initiator==2),
+                     alt1.domestic    =     expr(alt1.funding==1), 
+                     alt2.domestic    =     expr(alt2.funding==1),
+                     alt1.foreign     =     expr(alt1.funding==2), 
+                     alt2.foreign     =     expr(alt2.funding==2))
+
+
+
+
+
+
+#place your utility functions here
+u<-list(
+  v1 =V.1 ~  bprof*alt1.professional+ bexp * alt1.expert + bdomestic * alt1.domestic + bforeign * alt1.foreign + bdamage*alt1.damage + bprice * alt1.compensation,
+  v2 =V.2 ~  bprof*alt2.professional + bexp * alt2.expert + bdomestic * alt2.domestic + bforeign * alt2.foreign + bdamage*alt2.damage + bprice * alt2.compensation,
+  v3 =V.3 ~ basc)
+
\ No newline at end of file

parameters_SE_DRIVE.R

@@ -4,17 +4,25 @@
 
 designpath<- "SE_DRIVE/Designs/"
 
-resps =360  # number of respondents
-nosim= 100 # number of simulations to run (about 500 is minimum)
+resps =120  # number of respondents
+nosim= 500 # number of simulations to run (about 500 is minimum)
 
 
 
 
-## If beta values come from other source
 
-b1 = -0.036
-b2 = -0.034
-b3 = -0.049
+
+# bpreis = -0.036
+# blade  = -0.034
+# bwarte = -0.049
+
+
+# wrong parameters
+
+# 
+bpreis = -0.01
+blade = -0.07
+bwarte = 0.02
 
 manipulations = list(alt1.x2=     expr(alt1.x2/10), 
                      alt1.x3=     expr(alt1.x3/10),
@@ -25,6 +33,6 @@ manipulations = list(alt1.x2=     expr(alt1.x2/10),
 
 #place your utility functions here
 u<-list(
-  v1 =V.1~  b1 * alt1.x1 + b2*alt1.x2 + b3*alt1.x3   ,
-  v2 =V.2~  b1 * alt2.x1 + b2*alt2.x2 + b3*alt2.x3  
+  v1 =V.1~  bpreis * alt1.x1 + blade*alt1.x2 + bwarte*alt1.x3   ,
+  v2 =V.2~  bpreis * alt2.x1 + blade*alt2.x2 + bwarte*alt2.x3  
 )

simulationcore_purrr.R

@@ -3,9 +3,10 @@ rm(list=ls())
 source("functions.R")
 
 library("formula.tools")
+library(ggplot2)
+library(stringr)
 
-
-  #source("parameters_SE Design-Agri.R")
+#  source("parameters_SE Design-Agri.R")
 source("parameters_SE_DRIVE.R")