Commit 700ac282 authored by matbuoro's avatar matbuoro
Browse files

replace NA by average lambda

parent 029586ad
......@@ -60,6 +60,7 @@ inits_fix <- list(
# 180 # 2016
# )
lambda_tot = (data$Cm_Eu[,1] + data$Cm_Eu[,2]) * ((rowSums(data$Cm_B) + rowSums(data$Cum_B))/rowSums(data$Cm_B))
lambda_tot = ifelse(is.na(lambda_tot),mean(lambda_tot,na.rm=TRUE),lambda_tot)
# METTRE A JOUR /!\ TAILLE MATRICE
# logit_pi_B = structure(.Data = c(
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment