Hello,

I am running a single-season two-species model in RPresence using PsiBA parameterization. I am primarily interested in determining the effect of an apex predator on smaller carnivores. I am little bit caught up in formalizing the codes. For species-specific occupancy without any interaction where psiA, psiBA=psiBa. the code is :

mod1<-occMod(model=list(psi~SP,p~SP),data=data,type="so.2sp.1",parameterization="PsiBA")

The matrix of the model is:

a1 a2

psiA "1" "0"

psiBA "1" "1"

psiBa "1" "1"

Is it correct for estimating species-specific occupancy without any interactions? or will it should be like this:

a1 a2

psiA 1 0

psiBA 0 1

psiBa 0 1

How should I enter my codes for getting the above matrix where the a1 parameter of psiBA=psiBa=0. Does both the matrices signify the same thing that is occupancy without any interactions? if not what is the difference between the two.

Moreover, if anyone can tell me how will I predict my beta estimates generate through RPresence for plotting purposes.

Thank You