Parameter pi in HetClosed model

[mspinola]

I am fitting a Hetrogeneity model in RMark. See below. But I dont understand the output. If pi is set as pi=formula=~1. I need to have 2 p values, right? Or how I need to set pi to get 2 p values? Or how I need to set p to get 2 p values?

riqmark <- process.data(data=riq_m, model = "HetClosed", mixtures = 2, groups=c("manejo"))
riq.ddl=make.design.data(riqmark)

N.manejo=list(formula=~manejo)
p.dot=list(formula=~1)
pi.dot=list(formula=~1)

rich.N.manejo.p.dot.pi.dot= mark(riqmark, riq.ddl, model.parameters=list(N=N.manejo,p=p.dot, pi=pi.dot))

> rich.N.manejo.p.dot.pi.dot= mark(riqmark, riq.ddl, model.parameters=list(N=N.manejo,p=p.dot, pi=pi.dot))
STOP NORMAL EXIT

Output summary for HetClosed model
Name : pi(~1)p(~1)N(~manejo)

Npar : 3
-2lnL: 263.3825
AICc : 269.3942

Beta
estimate se lcl ucl
pi:(Intercept) 0.0000000 0.0000000 0.0000000 0.0000000
p:(Intercept) -1.2865473 0.0622208 -1.4085002 -1.1645945
N:(Intercept) 3.1505685 0.2597246 2.6415081 3.6596288
N:manejo1 -0.2303438 0.3399449 -0.8966358 0.4359482


Real Parameter pi
Group:manejo0

0.5

Group:manejo1

0.5


Real Parameter p
Group:manejo0
1
mixture:1 0.2164378
mixture:2 0.2164378

Group:manejo1
1
mixture:1 0.2164378
mixture:2 0.2164378


Real Parameter N
Group:manejo0

167.3493

Group:manejo1

133.5454

[jlaufenb]

You must specify mixtures in the submodel formula for p:

p.mix =list(formula=~mixture)
pi.dot =list(formula=~1)

See the Edwards and Eberhardt example in the RMark documentation:

?edwards.eberhardt

[mspinola]

Thank you very much. That answered my question.