Derived estimates of two groups

[Rudymena88]

Hello,

I'm running a Huggins closed population model with 2 groups (Male and Female) however my output for N only produces an estimate of both groups together:

Code: Select all

> head(Capturehistoryconc)
    ch tag.list sex
1 1000     F160   1
2 1000       F3   1
3 1100     F161   1
4 1011     M176   0
5 1000     F126   1
6 1000     F170   1
> str(Capturehistoryconc)
'data.frame':   39 obs. of  3 variables:
 $ ch      : chr  "1000" "1000" "1100" "1011" ...
 $ tag.list: chr  "F160" "F3" "F161" "M176" ...
 $ sex     : Factor w/ 2 levels "0","1": 2 2 2 1 2 2 2 2 2 2 ...
> dp=process.data(Capturehistoryconc,group="sex",model="HugHet")
> dp
$data
     ch tag.list sex group
1  1000     F160   1     2
2  1000       F3   1     2
3  1100     F161   1     2
4  1011     M176   0     1
5  1000     F126   1     2
6  1000     F170   1     2
7  1000     F186   1     2
8  1100      F35   1     2
9  1000     F171   1     2
10 1000     F134   1     2
$nocc
[1] 4

$nocc.secondary
NULL

$time.intervals
[1] 1 1 1

$begin.time
[1] 1

$age.unit
[1] 1

$initial.ages
[1] 0 0

$group.covariates
  sex
1   0
2   1
>ee.huggins.Mth2.additive=mark(Capturehistoryconc,model="HugFullHet",
                              model.parameters=list(p=ptimemixtureshared),adjust=TRUE)
>ee.huggins.Mth2.additive$results$derived
$`N Population Size`
  estimate       se      lcl      ucl
1 65.49802 11.99515 50.36829 100.7635



Is there any way of viewing estimates of N for both groups?

Thanks,

Rudy

[jlaake]

Use do as first argument of mark

[jlaake]

Was supposed to be dp as first argument. Doing from new phone and it auto corrected.

[Rudymena88]

Ah perfect!

This worked:

Code: Select all

ee.huggins.m0=mark(dp,model="Huggins",
                   model.parameters=list(p=pdotshared))

Thanks!

Rudy