I got the fatty acid data loaded again today, but I had some trouble and it took me a while to figure out what was up. I eventually tracked it down to these places and thought some warning messages / checks might help.
- DataIO: line 294: "stopifnot(dim(mean_c)[1]==n.preys & dim(mean_c)[2]==n.fats)": add something like, "number of fatty acids in Conv.Coeff.mean does not match number of fatty acids in the predators file" and/or "number of prey in Conv.Coeff.mean does not match number of preys in preys file".
- DataIO: line 294: "stopifnot(dim(var_c)[1]==n.preys & dim(var_c)[2]==n.fats)": "number of fatty acids in Conv.Coeff.var does not match number of fatty acids in the predators file" and/or "number of prey in Conv.Coeff.var does not match number of preys in preys file".
- DataIO: line 371: " R[,,i]=cov(alr(preys[preys.ix==unique(preys.ix)[i],]))*ni[i]". I got an error here, I think because there was one prey category for which there was only one prey item, so it couldn't do a covariance matrix? I removed that prey item and it went okay. If that is indeed the problem, maybe check to make sure there is more than one prey item in each category and throwing an error before doing the cov(alr(...)) call?
I got the fatty acid data loaded again today, but I had some trouble and it took me a while to figure out what was up. I eventually tracked it down to these places and thought some warning messages / checks might help.