Day 4 added

Author: psolymos

day3-2-detection.Rmd

@@ -714,7 +714,6 @@ AIC(Mmulti1, Mmulti2)
 
 tau1 <- exp(coef(Mmulti1))
 names(tau1) <- gsub("log.tau_spp", "", names(tau1))
-tau1
 
 X <- model.matrix(~ logmass + MaxFreqkHz + Mig2 + Hab2,
   lhreg_data[match(names(tau1), lhreg_data$spp),])
@@ -723,6 +722,8 @@ tau2 <- exp(X %*% coef(Mmulti2))
 res <- data.frame(count=sapply(yall_dis[names(tau1)], sum), tau1=tau1, tau2=tau2)
 res$diff <- res$tau1 - res$tau2
 
+res
+
 plot(diff ~ count, res)
 abline(h=0, lty=2)
 ```

day4-1-intro.Rmd

@@ -162,6 +162,7 @@ include_graphics("images/variable-selection-2.png")
 # Disadvantages
 
 - Taking care of uncertainty needs more work
+- Model validation gets a bit tricky
 
 ***