m1 <- model.dmc(
p.map = list(a="1",v="F",z="1",d="1",sz="1",sv="1", t0="1",st0="1"),
constants = c(st0=0,d=0),
match.map = list(M=list(s1="r1",s2="r2")),
factors = list(S=c("s1","s2"), F=c("f1", "f2")),
responses = c("r1","r2"),
type = "rd")
pVec <- c(a=1, v.f1=1, v.f2=1.5, z=0.5, sz=0.25, sv=0.2,t0=.15)
p.df.dmc(m1, pVec, "s1.f1.r1", n1order=TRUE)
Run the code above in your browser using DataLab