-
Notifications
You must be signed in to change notification settings - Fork 4
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
revise model_read_lib to output visit-level covariates for PC data #106
Comments
I can work on this over the next few days. I don't see temp and wind in the PC file, so like you said, I think we need to go to Michael Park's database to access them. |
@mkclapp I think that I was able to create indices for the observer in read_model_lib. It's currently in the It makes a matrix similar to that of y.ind, but instead of the number of counts it has the observer name. I can work on adding the time and date to the "data list," but I just want to be sure that this formatting would work when indexing it in the jags model. The date and time would have similar formats (at least in the way that I was planning on implementing them). |
I need to export the relevant data T and wind and we will need to modify the readPC script to accommodate. Finally we can tidy this up and have a separate table of visits with covariates in order to do data imputation etc. A separate issue is to get and clean the metadata from the ARU temperatures. |
(after a) Happy Thanksgiving y'all! |
I just pushed my changes to I updated model_read_lib to include the observer, date, and time for PC visits. Date and Time are standardized like they were in @mkclapp 's manual script. I'll work on doing a similar thing for the ARUs, but I'm not sure if I should put them in a vector similar to that of the |
Automatically generate data list with visit-level data for PC data, for use in adding visit-level covariates to the detection level of the model (ie, covariates on p11).
(Right now, we don't have an analogous parameter for the ML detections. We could add this, and use the start time of the file, the mean of the SongMeter temperature log during that time (or morning, if aggregating by day), etc... but I think that should be another issue).
I have code to do this manually, but it would be nice to have it all generated along with "data".
Creating a standalone list might be good so we can have flexibility to include various combinations of covariates in the models and compare them.
List of candidate covariates:
-date
-time of morning
-observer (relevant for point count data/p11 only, not ARU data)
-temperature*
-wind speed*
*I am not sure whether these two variables are available to us at the moment; we may need to extract them from our point count database.
Covariates should generally be standardized to mean = 0. I am going to branch off this issue and add a script that includes my manual way of generating these.
(TODO: I need to remind myself how categorical covariates (such as observer identity) need to be treated in a jags data object; for now, just use numerical covariates)
The text was updated successfully, but these errors were encountered: