DRIVE project added

SE_DRIVE/Designs/effdesignconstr.ngd

+Design	Choice situation	alt1.x1	alt1.x2	alt1.x3	alt2.x1	alt2.x2	alt2.x3	Block	
+1	1	40	50	50	80	25	50	5	
+1	2	20	200	25	40	50	100	6	
+1	3	80	50	0	20	50	100	1	
+1	4	40	50	100	60	100	0	3	
+1	5	60	25	50	20	200	50	1	
+1	6	40	200	100	60	25	100	2	
+1	7	60	100	0	40	200	100	8	
+1	8	60	50	50	40	200	25	10	
+1	9	60	100	100	80	50	0	5	
+1	10	80	25	25	20	200	0	1	
+1	11	20	200	100	60	200	50	7	
+1	12	60	100	100	80	100	50	4	
+1	13	40	50	50	80	25	25	3	
+1	14	80	25	0	40	200	0	2	
+1	15	20	25	50	80	25	25	8	
+1	16	80	25	25	40	50	100	4	
+1	17	40	100	50	60	25	0	2	
+1	18	60	50	25	40	25	50	3	
+1	19	60	50	0	20	200	0	8	
+1	20	60	25	100	80	25	0	6	
+1	21	80	25	50	20	200	25	2	
+1	22	80	25	50	60	100	100	8	
+1	23	40	25	100	60	50	25	7	
+1	24	40	100	100	80	100	50	9	
+1	25	60	25	25	20	200	0	3	
+1	26	80	50	25	20	50	100	10	
+1	27	60	100	0	20	25	100	7	
+1	28	40	50	100	20	200	50	7	
+1	29	20	200	25	80	25	50	8	
+1	30	20	50	100	80	25	0	8	
+1	31	60	50	100	20	100	50	4	
+1	32	40	25	50	60	50	0	2	
+1	33	80	100	0	40	100	25	8	
+1	34	80	25	50	60	50	50	1	
+1	35	60	200	100	80	25	50	10	
+1	36	80	25	0	20	200	25	10	
+1	37	40	50	100	60	50	25	1	
+1	38	60	200	100	80	25	25	7	
+1	39	20	200	25	60	25	25	1	
+1	40	20	200	0	40	50	50	6	
+1	41	20	100	50	60	25	100	7	
+1	42	80	25	25	60	200	50	7	
+1	43	20	50	50	60	50	0	2	
+1	44	20	100	50	40	50	0	5	
+1	45	40	100	0	80	50	50	3	
+1	46	40	200	100	60	25	0	10	
+1	47	40	25	25	20	200	100	1	
+1	48	20	200	100	60	25	100	4	
+1	49	20	100	50	60	50	50	8	
+1	50	60	100	50	40	200	25	7	
+1	51	80	50	0	40	50	100	3	
+1	52	20	100	25	80	25	100	5	
+1	53	40	200	0	60	50	25	1	
+1	54	40	50	50	80	50	25	9	
+1	55	80	100	50	20	200	50	1	
+1	56	60	200	0	80	50	100	3	
+1	57	20	200	0	20	100	100	3	
+1	58	40	25	0	20	100	100	2	
+1	59	20	200	100	40	200	25	8	
+1	60	40	50	25	20	100	0	7	
+1	61	60	25	100	20	100	25	9	
+1	62	40	200	0	80	100	100	5	
+1	63	80	50	25	40	100	100	6	
+1	64	60	50	25	80	50	0	1	
+1	65	80	50	25	40	25	50	10	
+1	66	60	25	25	40	200	50	6	
+1	67	60	25	25	20	25	50	5	
+1	68	20	25	100	60	100	25	4	
+1	69	20	200	25	60	50	0	6	
+1	70	60	100	50	20	100	100	2	
+1	71	40	200	100	80	200	25	4	
+1	72	80	25	100	20	200	100	2	
+1	73	20	25	25	80	25	0	5	
+1	74	60	25	0	40	100	50	9	
+1	75	40	50	50	20	200	100	4	
+1	76	60	25	100	60	100	25	7	
+1	77	80	50	50	20	25	100	9	
+1	78	60	100	50	20	200	50	6	
+1	79	80	25	100	40	100	25	9	
+1	80	40	50	50	20	200	25	8	
+1	81	40	200	0	60	25	100	10	
+1	82	60	25	50	80	50	0	4	
+1	83	20	100	25	80	100	0	2	
+1	84	60	50	0	20	25	100	9	
+1	85	80	50	100	60	100	50	3	
+1	86	80	25	50	20	50	100	4	
+1	87	40	25	100	20	50	100	6	
+1	88	20	100	25	80	50	0	4	
+1	89	40	50	0	60	25	0	5	
+1	90	60	25	50	40	100	50	7	
+1	91	20	25	100	60	50	25	3	
+1	92	80	25	0	60	200	50	2	
+1	93	80	100	0	40	200	25	5	
+1	94	20	100	25	40	50	50	10	
+1	95	20	200	100	80	25	0	5	
+1	96	20	200	100	80	50	0	5	
+1	97	80	200	0	20	100	100	3	
+1	98	60	100	25	40	200	0	1	
+1	99	20	25	100	60	50	50	10	
+1	100	40	200	50	60	200	0	9	
+1	101	60	50	0	20	50	100	8	
+1	102	60	200	50	40	200	100	8	
+1	103	20	200	25	40	25	0	4	
+1	104	40	100	0	20	100	50	6	
+1	105	40	50	0	20	100	25	9	
+1	106	60	100	50	80	100	0	6	
+1	107	60	25	25	20	200	25	10	
+1	108	60	100	50	40	50	100	2	
+1	109	80	100	0	20	100	50	9	
+1	110	20	100	0	80	25	50	10	
+1	111	80	50	0	60	200	25	9	
+1	112	60	200	100	80	50	25	6	
+1	113	20	200	0	40	25	100	1	
+1	114	20	200	0	80	100	50	9	
+1	115	20	200	50	40	50	100	4	
+1	116	80	100	25	20	100	100	7	
+1	117	40	200	25	80	25	25	3	
+1	118	80	50	25	40	200	50	5	
+1	119	20	100	100	80	50	0	6	
+1	120	80	25	25	40	100	25	10	
+||||||||||
+design
+;alts = alt1, alt2
+;block = 10
+;eff = (mnl, d)
+;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[-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

SE_DRIVE/effdesign_norestr.ngd → SE_DRIVE/Designs/effdesignnorestr.ngd

-Design	Choice situation	opt1.x1	opt1.x2	opt1.x3	opt2.x1	opt2.x2	opt2.x3	Block	
+Design	Choice situation	alt1.x1	alt1.x2	alt1.x3	alt2.x1	alt2.x2	alt2.x3	Block	
 1	1	60	200	100	80	25	0	10	
 1	2	40	25	0	20	200	100	8	
 1	3	40	200	100	60	25	0	9	
@@ -121,18 +121,18 @@ Design	Choice situation	opt1.x1	opt1.x2	opt1.x3	opt2.x1	opt2.x2	opt2.x3	Block
 1	120	20	25	50	80	200	0	4	
 ||||||||||
 design
-;alts = opt1*, opt2*
+;alts = alt1*, alt2*
 ;block = 10
 ;eff = (mnl, d)
 ;alg = mfederov
 ;rows = 120
 
 ;model:
-U(opt1) = b1[-0.036]  * x1[20,40,60,80]
+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(opt2) = b1 * x1
+U(alt2) = b1 * x1
         + b2 * x2
         + b3 * x3
 $
\ No newline at end of file

SE_DRIVE/Designs/olddesign.ngd

+"Choice situation"	"Block"	"alt1.x1"	"alt2.x1"	"alt1.x2"	"alt2.x2"	"alt1.x3"	"alt2.x3"	"Design"
+97	9	20	40	200	25	100	0	1
+98	9	40	60	200	25	100	0	1
+99	9	20	80	200	25	100	25	1
+100	9	20	40	200	50	25	50	1
+101	9	80	60	25	100	25	50	1
+102	9	40	80	100	25	100	0	1
+106	9	60	20	25	100	100	50	1
+105	9	20	40	100	100	100	50	1
+104	9	20	40	50	25	100	100	1
+103	9	80	20	25	200	100	100	1
+108	9	80	40	50	200	50	25	1
+107	9	40	20	50	25	25	100	1
+49	5	80	40	50	25	50	100	1
+50	5	60	40	50	100	50	25	1
+51	5	60	20	50	100	25	50	1
+52	5	60	40	100	100	0	50	1
+53	5	20	80	100	50	50	50	1
+54	5	60	20	50	100	0	50	1
+58	5	40	80	25	100	100	0	1
+57	5	60	20	25	50	0	25	1
+56	5	80	20	100	200	0	25	1
+55	5	40	60	200	50	50	0	1
+60	5	40	20	25	200	50	25	1
+59	5	80	60	25	200	0	0	1
+37	4	20	60	200	100	100	100	1
+38	4	60	80	200	100	0	50	1
+39	4	40	80	200	25	25	0	1
+40	4	80	20	25	100	50	0	1
+41	4	80	40	25	50	100	100	1
+42	4	60	20	50	25	0	100	1
+46	4	60	20	100	200	25	0	1
+45	4	80	60	25	200	0	50	1
+44	4	20	80	50	25	25	25	1
+43	4	40	60	100	100	50	25	1
+48	4	80	20	50	200	25	0	1
+47	4	40	60	50	100	100	25	1
+109	10	80	40	25	200	100	0	1
+110	10	60	20	25	200	50	50	1
+111	10	60	40	25	200	100	25	1
+112	10	40	80	200	100	50	0	1
+113	10	40	60	200	25	100	50	1
+114	10	40	80	100	25	50	50	1
+118	10	40	20	200	200	25	50	1
+117	10	60	80	50	50	100	50	1
+116	10	80	60	25	50	50	50	1
+115	10	80	40	50	100	25	25	1
+120	10	20	80	50	50	50	0	1
+119	10	20	60	25	25	100	50	1
+85	8	20	40	100	50	25	25	1
+86	8	20	60	200	50	25	50	1
+87	8	80	20	50	200	50	25	1
+88	8	60	80	50	25	50	0	1
+89	8	40	20	25	200	50	100	1
+90	8	40	60	100	25	25	25	1
+94	8	40	80	200	25	0	0	1
+93	8	80	40	25	25	0	25	1
+92	8	20	60	200	100	0	50	1
+91	8	60	80	200	100	50	0	1
+96	8	40	20	50	200	0	100	1
+95	8	40	20	50	100	50	50	1
+61	6	40	20	25	100	0	100	1
+62	6	80	60	25	100	50	25	1
+63	6	20	40	50	50	50	25	1
+64	6	80	60	50	100	50	50	1
+65	6	60	40	25	25	0	100	1
+66	6	80	60	50	200	50	0	1
+70	6	60	20	25	50	0	0	1
+69	6	80	40	25	200	25	0	1
+68	6	60	20	25	50	0	50	1
+67	6	80	20	50	200	50	0	1
+72	6	40	60	100	50	50	0	1
+71	6	20	60	200	200	50	0	1
+73	7	80	20	25	200	50	50	1
+74	7	40	60	200	25	50	100	1
+75	7	40	60	50	50	25	0	1
+76	7	40	60	200	25	25	25	1
+77	7	20	80	100	100	50	0	1
+78	7	20	80	25	50	100	0	1
+82	7	20	80	200	25	25	0	1
+81	7	80	40	25	200	25	100	1
+80	7	20	40	50	100	100	0	1
+79	7	40	80	200	100	50	100	1
+84	7	60	20	25	50	0	100	1
+83	7	40	60	200	100	25	25	1
+25	3	20	80	200	25	50	25	1
+26	3	80	60	50	100	25	25	1
+27	3	20	40	100	25	0	50	1
+28	3	60	80	25	25	100	0	1
+29	3	80	60	25	25	0	25	1
+30	3	40	60	200	25	0	100	1
+34	3	60	20	200	200	0	25	1
+33	3	40	20	200	200	0	50	1
+32	3	80	60	25	200	25	25	1
+31	3	40	80	100	50	100	25	1
+36	3	60	80	200	200	100	50	1
+35	3	20	40	100	50	100	0	1
+13	2	40	80	200	50	25	0	1
+14	2	40	60	100	25	100	50	1
+15	2	80	40	100	100	0	25	1
+16	2	60	20	25	50	25	25	1
+17	2	20	40	200	25	50	50	1
+18	2	60	40	50	100	25	100	1
+22	2	20	80	25	25	100	50	1
+21	2	20	80	200	100	100	50	1
+20	2	20	40	50	50	100	25	1
+19	2	80	40	25	25	25	100	1
+24	2	40	20	100	200	50	50	1
+23	2	40	20	50	200	100	0	1
+1	1	80	20	25	200	100	50	1
+2	1	60	40	50	100	50	100	1
+3	1	60	20	200	200	0	100	1
+4	1	20	80	200	25	0	100	1
+5	1	40	80	100	50	100	50	1
+6	1	60	80	50	25	0	0	1
+10	1	40	60	50	25	100	0	1
+9	1	80	40	100	200	100	25	1
+8	1	40	20	25	25	25	100	1
+7	1	20	60	25	25	100	0	1
+12	1	20	60	200	25	50	100	1
+11	1	20	60	200	25	50	0	1

parameters_SE_DRIVE.R

-
-
-source("prior_est/scripts/priordetermination_parks.R")
-
-designpath<- "Designsparks/"
-
-resps =360  # number of respondents
-nosim=2000 # number of simulations to run (about 500 is minimum)
-
-
-
-
-## If beta values come from other source
-for (ano in 1:length(priors)) {
-  assign(names(priors[ano]),priors[ano]) 
-}
-
-#place your utility functions here
-u<-list(
-  v1 =V.1~  basc        +bgroesse*alt1.groesse+ bentfernung * alt1.entfernung + bgemeinschaft * alt1.gemeinschaft + bkultur * alt1.kultur + bumweltbildung*alt1.umweltbildung + btoiletten * alt1.toiletten+ bspiel*alt1.spiel + bpflegeint*alt1.pflegeint+bpflegeziele*alt1.pflegeziele+ bbeitrag*alt1.beitrag  ,
-  v2 =V.2~  basc        +bgroesse*alt2.groesse+ bentfernung * alt2.entfernung + bgemeinschaft * alt2.gemeinschaft + bkultur * alt2.kultur + bumweltbildung*alt2.umweltbildung + btoiletten * alt2.toiletten+ bspiel*alt2.spiel + bpflegeint*alt2.pflegeint+bpflegeziele*alt2.pflegeziele+ bbeitrag*alt2.beitrag,
-  v3 =V.3~ 0)
+
+
+
+
+designpath<- "SE_DRIVE/Designs/"
+
+resps =360  # number of respondents
+nosim= 100 # 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
+
+manipulations = list(alt1.x2=     expr(alt1.x2/10), 
+                     alt1.x3=     expr(alt1.x3/10),
+                     alt2.x2=     expr(alt2.x2/10), 
+                     alt2.x3=     expr(alt2.x3/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  
+)

simulation_output.rmd

 ---
 title: "Simulation experimental design"
-author: "Julian Sagebiel"
+author: ""
 date: "2023-01-05"
 format:
   html:

simulationcore_purrr.R

@@ -5,8 +5,8 @@ source("functions.R")
 library("formula.tools")
 
 
-#source("parameters_SE Design-Agri.R")
-source("parameters_ip2.R")
+  #source("parameters_SE Design-Agri.R")
+source("parameters_SE_DRIVE.R")
 
 
 designfile<-list.files(designpath,full.names = T)