{"action":"opened","issue":{"url":"https://api.github.com/repos/kunstler/FUNDIV_sAPROPOS/issues/43","repository_url":"https://api.github.com/repos/kunstler/FUNDIV_sAPROPOS","labels_url":"https://api.github.com/repos/kunstler/FUNDIV_sAPROPOS/issues/43/labels{/name}","comments_url":"https://api.github.com/repos/kunstler/FUNDIV_sAPROPOS/issues/43/comments","events_url":"https://api.github.com/repos/kunstler/FUNDIV_sAPROPOS/issues/43/events","html_url":"https://github.com/kunstler/FUNDIV_sAPROPOS/issues/43","id":307740366,"number":43,"title":"Evaluate IPM integration in extreme climate","user":{"login":"kunstler","id":4505188,"avatar_url":"https://avatars1.githubusercontent.com/u/4505188?v=4","gravatar_id":"","url":"https://api.github.com/users/kunstler","html_url":"https://github.com/kunstler","followers_url":"https://api.github.com/users/kunstler/followers","following_url":"https://api.github.com/users/kunstler/following{/other_user}","gists_url":"https://api.github.com/users/kunstler/gists{/gist_id}","starred_url":"https://api.github.com/users/kunstler/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/kunstler/subscriptions","organizations_url":"https://api.github.com/users/kunstler/orgs","repos_url":"https://api.github.com/users/kunstler/repos","events_url":"https://api.github.com/users/kunstler/events{/privacy}","received_events_url":"https://api.github.com/users/kunstler/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":0,"created_at":"2018-03-22T17:22:38Z","updated_at":"2018-03-22T17:22:38Z","closed_at":null,"author_association":"OWNER","body":"To check that my implementation of Dylan trick to integrate the IPM work well even when we predict growth and survival in extreme climate (min and max of SGDD sum of growing degree days of the species) I compared the sum of one size at t of the P matrix with the survival (should match if we integrate well growth).\r\n\r\nTo summarise what the function `mk_P_intMulti_Diag` does. It simply does a two-dimensional integration with a Gauss-Legendre integration as described in the IPM book (level 3 in one direction and 17 on the other). This done only close to the diagonal as the rest of the matrix is ~ 0 for slow-growing trees. For the rest, the mid bin method is used. (so between the black lines below)\r\n\r\n![image](https://user-images.githubusercontent.com/4505188/37786546-445187be-2dfd-11e8-9969-defe69643b72.png)\r\n\r\n\r\nHere is in an example for *Fagus sylvatica* that seems to work well (red lines are the probability of survival and black points sum of P matrix for that size). \r\n\r\n![image](https://user-images.githubusercontent.com/4505188/37786540-404cb9ea-2dfd-11e8-9ac2-7bf390f870b4.png)\r\n\r\nThere is just some eviction left for size close to the max (just because the growth kernel predicts growth outside the range). I'm not sure if this is a problem or not and if we should correct for that (by adding the missing probability to the last size as it is done in IPMpack). @dzchilds and @robcito What do you think of that? Is it ok should I try to solve that further?"}}
{"actor":{"display_login":"kunstler"}}