Double Mating
Anastasia Shavrova
2022-03-22
Fecundity Over Time
- trial = experimental trial (1-8)
- mate_order = male morph mating order in either of 4 treatments: “FF”, “FS”, “SS”, “SF”
- day = the day of egg laying starting at the mating day relates to female longevity
- egg_count = the number of eggs the female laid
- times_mated = how many times a female mated either 1 or 2 times
- replicate = replicate vial within a treatment (1-5)
- indv = unique identification number for female
- morph = the morph of the male mite (scrambler or fighter)
- order = mating order of the mite (1st or 2nd)
Only using females mated twice
Did females who mated twice have a different rate of fecundity depending on treatment?
log(egg_count+1)~day+I(day^2)x(mate order)+(1|indv)+(1|trial)
Differences arise between treatments: slowest to fastest rate: FF>SS>SF>FS
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## day 5.354 5.354 1 2454.31 2.6162 0.1059
## I(day^2) 82.814 82.814 1 2456.62 40.4668 2.377e-10 ***
## mate_order 5.243 1.748 3 125.22 0.8540 0.4670
## I(day^2):mate_order 147.611 49.204 3 2488.69 24.0434 2.458e-15 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Note: Use 'contrast(regrid(object), ...)' to obtain contrasts of back-transformed estimates## $emtrends
## mate_order (day)^2.trend SE df lower.CL upper.CL
## FF 1.97e-05 0.000147 2476 -0.000268 3.08e-04
## FS -4.92e-04 0.000235 2490 -0.000954 -3.12e-05
## SF -1.50e-03 0.000202 2490 -0.001899 -1.11e-03
## SS -1.35e-03 0.000203 2489 -0.001750 -9.56e-04
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## FF - FS 0.000512 0.000265 2488 1.932 0.2147
## FF - SF 0.001523 0.000217 2489 7.014 <.0001
## FF - SS 0.001373 0.000224 2489 6.124 <.0001
## FS - SF 0.001011 0.000306 2490 3.306 0.0053
## FS - SS 0.000861 0.000307 2489 2.804 0.0261
## SF - SS -0.000150 0.000277 2491 -0.542 0.9486
##
## Note: contrasts are still on the 1(mu + 1) scale
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 4 estimatesFirst Day Fecundity
Does time of first egg lay change due to treatment?
day~mate_order + egg_count [covariate] + (1|trial)
Female first egg lay day does not differ between treatments
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## mate_order 62.790 20.930 3 234.21 1.2700 0.28534
## egg_count 71.887 71.887 1 237.98 4.3622 0.03781 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## $emmeans
## mate_order emmean SE df lower.CL upper.CL
## FF 5.89 1.63 8.25 2.154 9.63
## FS 5.53 1.63 8.20 1.797 9.27
## SF 4.45 1.63 8.34 0.709 8.19
## SS 5.23 1.62 8.10 1.497 8.96
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## FF - FS 0.359 0.754 234 0.476 0.9643
## FF - SF 1.441 0.775 235 1.858 0.2490
## FF - SS 0.663 0.741 234 0.895 0.8073
## FS - SF 1.082 0.746 234 1.450 0.4698
## FS - SS 0.304 0.720 234 0.423 0.9746
## SF - SS -0.778 0.737 234 -1.056 0.7168
##
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 4 estimatesDoes amount of eggs laid on the first day change due to treatment?
egg_count~mate_order + day [covariate] + (1|trial)
There are differences between how many eggs are laid on the first day of egg lay between treatments, from greatest to least SF>FS>SS>FF
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## mate_order 7413.6 2471.2 3 234.15 4.4617 0.004542 **
## day 1973.7 1973.7 1 232.94 3.5635 0.060307 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## $emmeans
## mate_order emmean SE df lower.CL upper.CL
## FF 25.6 6.48 9.98 11.1 40.0
## FS 37.5 6.46 9.84 23.1 51.9
## SF 40.8 6.50 10.11 26.3 55.2
## SS 35.7 6.41 9.59 21.3 50.0
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## FF - FS -11.94 4.30 234 -2.774 0.0303
## FF - SF -15.21 4.42 235 -3.443 0.0038
## FF - SS -10.08 4.25 234 -2.373 0.0853
## FS - SF -3.27 4.34 235 -0.754 0.8749
## FS - SS 1.85 4.17 234 0.444 0.9707
## SF - SS 5.13 4.27 235 1.201 0.6270
##
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 4 estimatesPeak Fecundity
Does time of peak fecundity change due to treatment?
day~mate_order + egg_count [covariate] + (1|trial)
Yes timing of peak fecundity differs between treatments, from earliest day to latest day: SS<SF<FS<FF
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## mate_order 308.112 102.704 3 242.25 2.7418 0.04388 *
## egg_count 12.937 12.937 1 244.52 0.3454 0.55729
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## $emmeans
## mate_order emmean SE df lower.CL upper.CL
## FF 12.15 2.01 8.96 7.60 16.7
## FS 10.87 2.00 8.73 6.33 15.4
## SF 9.56 2.02 9.09 5.00 14.1
## SS 9.34 2.00 8.72 4.81 13.9
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## FF - FS 1.277 1.10 243 1.159 0.6531
## FF - SF 2.593 1.14 242 2.284 0.1047
## FF - SS 2.805 1.10 243 2.554 0.0544
## FS - SF 1.316 1.11 243 1.187 0.6358
## FS - SS 1.528 1.06 242 1.446 0.4721
## SF - SS 0.212 1.10 242 0.193 0.9974
##
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 4 estimatesDoes peak fecundity change due to treatment?
egg_count~mate_order + day[covariate] + (1|trial)
Peak fecundity does not differ between treatments
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## mate_order 4886.2 1628.74 3 243.85 0.8605 0.4622
## day 248.6 248.55 1 208.15 0.1313 0.7174## $emmeans
## mate_order emmean SE df lower.CL upper.CL
## FF 83.4 8.12 17.2 66.3 100.6
## FS 72.4 7.89 15.2 55.6 89.2
## SF 81.2 8.16 17.6 64.0 98.4
## SS 75.3 7.89 15.4 58.5 92.1
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## FF - FS 11.07 7.81 244 1.417 0.4897
## FF - SF 2.27 8.15 244 0.278 0.9925
## FF - SS 8.18 7.88 243 1.037 0.7280
## FS - SF -8.81 7.87 245 -1.118 0.6786
## FS - SS -2.90 7.54 243 -0.385 0.9806
## SF - SS 5.91 7.80 243 0.757 0.8735
##
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 4 estimatesLast Day Fecund
Does day of last lay differ between treatments?
day~mate_order + egg_count [covariate] + (1|trial)
Day of last egg lay does differ between treatments, from earliest to latest SS<FS<SF<FF
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## mate_order 1315.8 438.59 3 241.19 4.470 0.004472 **
## egg_count 2547.4 2547.37 1 246.61 25.962 6.929e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## $emmeans
## mate_order emmean SE df lower.CL upper.CL
## FF 23.1 2.11 12.3 18.5 27.7
## FS 18.8 2.14 12.8 14.2 23.4
## SF 19.1 2.16 13.5 14.5 23.8
## SS 16.8 2.12 12.5 12.2 21.4
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## FF - FS 4.30 1.79 242 2.405 0.0789
## FF - SF 3.96 1.82 242 2.172 0.1340
## FF - SS 6.29 1.76 242 3.576 0.0024
## FS - SF -0.34 1.82 240 -0.187 0.9977
## FS - SS 1.99 1.77 241 1.125 0.6745
## SF - SS 2.33 1.82 241 1.281 0.5761
##
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 4 estimatesDoes fecundity differ on the last day between treatments?
Last egg lay does change between morphs and mating order, from greatest to least SF>FS>FF>SS
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## mate_order 9129.7 3043.2 3 241.11 3.5184 0.0158 *
## day 22361.2 22361.2 1 247.00 25.8528 7.286e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## $emmeans
## mate_order emmean SE df lower.CL upper.CL
## FF 28.2 6.07 13.4 15.11 41.2
## FS 30.8 6.10 13.6 17.68 43.9
## SF 35.7 6.17 14.3 22.50 48.9
## SS 19.1 6.04 13.1 6.03 32.1
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## FF - FS -2.61 5.37 241 -0.486 0.9621
## FF - SF -7.53 5.44 241 -1.384 0.5108
## FF - SS 9.11 5.32 241 1.712 0.3197
## FS - SF -4.92 5.39 241 -0.912 0.7983
## FS - SS 11.72 5.21 240 2.248 0.1134
## SF - SS 16.65 5.32 241 3.131 0.0105
##
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 4 estimatesTotal Fecundity
- trial = experimental trial (1-8)
- mate_order = male morph mating order in either of 4 treatments: “FF”, “FS”, “SS”, “SF”
- day = the day of egg laying starting at the mating day relates to female longevity
- egg_count = the number of eggs the female laid
- times_mated = how many times a female mated either 1 or 2 times
- replicate = replicate vial within a treatment (1-5)
- indv = unique identification number for female
- total = total fecundity per female
- day = longevity of the female
- morph = the morph of the male mite (scrambler or fighter)
- order = mating order of the mite (1st or 2nd)
Total fecundity does not differ between treatments
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## mate_order 254798 84933 3 112.48 0.8189 0.486## $emmeans
## mate_order emmean SE df lower.CL upper.CL
## FF 585 70.8 27.9 440 730
## FS 490 70.5 26.8 345 634
## SF 529 72.0 29.1 381 676
## SS 465 69.2 24.9 322 607
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## FF - FS 95.8 83.0 113 1.154 0.6572
## FF - SF 56.7 84.4 113 0.672 0.9076
## FF - SS 120.9 81.9 113 1.476 0.4554
## FS - SF -39.1 83.5 113 -0.468 0.9658
## FS - SS 25.1 80.7 112 0.311 0.9895
## SF - SS 64.2 82.3 113 0.780 0.8634
##
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 4 estimatesSurvival
Survival does not differ between treatments
## Analysis of Deviance Table (Type II tests)
##
## Response: Surv(day)
## Df Chisq Pr(>Chisq)
## mate_order 3 7.6117 0.054757 .
## egg_count 1 7.1182 0.007631 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Mating Behaviour
- trial = experimental trial (1-8)
- mate.order = First mating with first morph, second mating with second morph
- male.morph = male morph as fighter (F) or scrambler (S) in either of 4 treatments: “FF”, “FS”, “SS”, “SF”
- time = time to initiate mating (in seconds)
- replicate = replicate vial within a treatment (1-8)
- mate.type = the type of mating pair either the same (FF/SS) or different (SF/FS)
- male = male morph presented to female: “Fighter” or “Scrambler”
Mating Rejection
Do females reject one morph over another, depending on mating order?
Females do not reject one morph more than another
Mating Latency
Using data where females mated twice
Does mating latency differ between morph type and mating order?
No mating latency does not differ but this is marginally not significant
## Type III Analysis of Variance Table with Satterthwaite's method
## Sum Sq Mean Sq NumDF DenDF F value Pr(>F)
## male 2.5 2.5 1 295.64 0.0047 0.9453306
## as.factor(mate.order) 5877.0 5877.0 1 295.14 11.1678 0.0009395 ***
## male:as.factor(mate.order) 1770.3 1770.3 1 295.55 3.3639 0.0676437 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## $emmeans
## male mate.order emmean SE df lower.CL upper.CL
## Fighter 1 37.6 3.94 16.0 29.2 45.9
## Scrambler 1 32.6 3.91 15.5 24.2 40.9
## Fighter 2 41.5 3.93 15.8 33.2 49.8
## Scrambler 2 46.1 3.92 15.6 37.8 54.5
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## Fighter 1 - Scrambler 1 5.00 3.72 296 1.344 0.5357
## Fighter 1 - Fighter 2 -3.95 3.74 295 -1.056 0.7167
## Fighter 1 - Scrambler 2 -8.59 3.73 295 -2.304 0.0994
## Scrambler 1 - Fighter 2 -8.95 3.70 295 -2.417 0.0761
## Scrambler 1 - Scrambler 2 -13.59 3.69 295 -3.683 0.0016
## Fighter 2 - Scrambler 2 -4.64 3.71 295 -1.249 0.5962
##
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 4 estimates