-
Notifications
You must be signed in to change notification settings - Fork 0
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
Pollinators-bees system data #12
Comments
Yes, I'll upload a template with a few rows this coming week and updated it
when perfectly clean. (still travelling, this is never ending!)
On Sun, Dec 4, 2016 at 5:45 PM Steve Vissault ***@***.***> wrote:
Assigned #12
<#12> to
@ibartomeus <https://github.com/ibartomeus>.
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#12 (comment)>,
or mute the thread
<https://github.com/notifications/unsubscribe-auth/AB8IgmI6W8ZksbQ91XqAcDxq053-xRDRks5rEu4vgaJpZM4LDoPd>
.
--
--
Ignasi Bartomeus PhD
www.bartomeuslab.com
@ibartomeus
Skype: nachobartomeus
Dpto. Ecología Integrativa
Estación Biológica de Doñana (EBD-CSIC)
Avda. Américo Vespucio s/n
Isla de la Cartuja
41092, Sevilla (Spain)
|
Perfect, let me know if you need help (by writing on this issue). That was nice to meet you! Have a good trip! |
Data available at : https://github.com/bw4sz/Montreal_Workshop/tree/master/Data Interactions, plant and pollinator datasets can be merged, should I upload just one unique dataset? Or the three tables separately? Other issues:
|
Hi @ibartomeus, I started to work on your issue before leaving for vacation. I'm still working on it. What I'm doing is this: I'm writing an update to your markdown document (montreal_workshop repo) to add a new section which reshape your data in order to get the data input format needed by the function |
Great, I'll have updated and hopefully final data (same format) next week. Happy to help! |
@SteveViss I don't see any update in the M_w repo. If you tell me the data format needed I can try to do it this week and get data part done! |
@ibartomeus, sorry for the delay. I restart to work on this today, I'll push before the end of the day. Can you tell me what means |
Intertegular distance. A proxy of body size. Just let me know what do you need next. |
library(devtools)
install_github("TheoreticalEcosystemEcology/alienR", ref = "alienData")
library("alienR")
load("data/argsAlienData.RData") # How I created this file is explained in the explnation.Rmd. File is in attachment
dataBartomeus <- as.alienData(idObs=idObs,interactPair=interactPair,traitInd=traitInd,traitSp=traitSp)
str(dataBartomeus)
head(dataBartomeus$idObs)
head(dataBartomeus$interactSp) #being this a bipartite netw, you don't want to have a square matrix here
dim(dataBartomeus$interactSp)
head(dataBartomeus$interactInd)
dim(dataBartomeus$interactInd) #same here, this is quite a waste of space, because most indiv won't interact with others, and certainly not with themselves
head(dataBartomeus$coOcc) #this I think is not well calculated. SHould be from idObs
head(dataBartomeus$traitSp)
head(dataBartomeus$traitInd) |
@ibartomeus, you're the specialist. So I'll include of these comments. |
Sorry I should have used the issue and not the mail. some details: I don't get how cOcc is created from the interaction matrix (interactPair) as this do not has any information on where and when. I expect co-ocurrence ignores real interactions and only uses if two species are present at the same time (Round in my case) and place (site in my case). Right? This would need to be amended if I am right. as.alienData also create potential interaction matrices that are squared. This lead to a lot of structural NA's in bipartite networks, which IMO are suboptimal. |
* Added the as.alienData function with help file written in roxygen2 * CamelCase functions * restore Bartomeus data file in R/ * Fixing issue with weight_sd and weight_mean functions * removin weight_mean and weigthed_sd, already present in stats - close #17 * Rebuild package and add new clean target in makefile * fixinx syntax issues * code review; fix minor issues/syntax errors * Set mandatory idObs object in as.alienData() #21 * cleaning up data format * removing ressourceConsumer matrix * Details on roxygen documentation * Add data consistency among interacPair and idObs #21 * set build interac matrix #21 * Add interacSp and interacInd from interacPair * finish to reference interacSp and interacInd with idObs #21 * check consistency among traitSp and traitInd with idObs #21 * imp. coOcc methods: inferred or not; last checks on out objects; make sure coOcc rownames and colnames reffered to idObs * make sure siteEnv are referenced in idObs #21 * build doc of as.alienData function * Push last edits on as.alienData function * set traitSp and traitInd as long format * fixing functions errors * Fixing last issue with res * building docs; and close #21 * Add data and test it * Add comments on as.alienData posible enhancements * set Nacho's comments in issue #12 * set bartomeus in data folder * init as.alienData test * init test on integrity #18 * cover unit tests as.alienData() - issue #18 * add unit tests to data structure issue #18
This is much easier to make tests on matrices, than a long format Thanks, |
According to your comment @ibartomeus, the
I'm fine with that. |
Good for the coOcc! Regarding the interactSp, whatever is wise programmatically is ok with me. I think in long formats, but It0s not important if it works fine in the package. |
This part is done, have a look here: https://github.com/TheoreticalEcosystemEcology/alien/blob/alienData/R/as.alienData.R#L360-L372 |
Thanks, now looks good to me! |
@MoralesCastilla, Do you still working on the phylo? I'm not familiar with it but I might help you to include it in the |
Hi Steve, I'm working into assembling the pollinator phylogeny to run the different models with Nacho's data. The plant phylogeny is ready. I guess that for coding purposes we could use the plant phylogeny and a dummy pollinator phylogeny (until I manage to build a final version). I have seen that you coded the phylo param as a square matrix (I'm assuming you mean the patristic phylogenetic distance matrix). I wonder if we could input the phylogeny as a "phylo" object instead, which is the most common format for phylogenies in R (https://cran.rstudio.com/web/packages/ape/index.html). By doing so we would avoid loosing information and we can always transform the phylo object into a square phylogenetic distance matrix using the function cophenetic.phylo(). Let me know if you need me to upload the phylogenies as they are to the data folder (I'm happy to upload them as either phylo objects or square matrices). |
Sure, I'm fine with the Thanks @MoralesCastilla ! |
Hi @MoralesCastilla, any progress on the phylo object? Let me know if you need help. |
Hi @SteveViss, I have pushed a bit of code to include phylo objects in |
Hi @MoralesCastilla, Nice, I'll will have a look at it. I will open a new issue if any things are missings or have to be improve. Travel safe. |
@ibartomeus, when you have a moment, can you take a look at the |
Looks ok to me. @MoralesCastilla Where is the phylo for plants? Can we add it already? And which is the best way to construct the phylo for pollinators (using a distance metric of its taxonomy?) |
@ibartomeus in the |
Hi @ibartomeus,
It will be nice to include your data as part of the package (
data
folder).Hence, I can use them to do do some unit testing on package functionalities.
Let me know what do you think,
Cheers.
The text was updated successfully, but these errors were encountered: