I'm absolutely not brilliant with R yet and working with package SIAR.
1) SIAR reads the first column of a table.txt as "int" instead of
"num". Who
can help?
2) I therefore made tables with a data.frame() > that seems to work. But at
the end of the loading proces of the SIARmodel, I get a warning screen with
"Session Aborted". Everything crashes and I can start all over again.
This is my script:
/# loading SIAR:
library(SIAR)
#get CONUSMERtable (NB: I do not use location now, since 'help' told me
a
single group is more easy):> d15N<-c(13.3175, 13.8675, 12.8775, 15.0785, 15.351, 14.6755,
> 14.94002083, 14.2072, 13.4987, 13.6522, 14.8447, 13.4257, 13.5632,
> 15.4282, 11.4187, 14.4442, 13.5187, 12.7762, 14.8627, 13.1272)
> d13C<-c(-23.4585, -23.078, -22.386, -17.8665, -16.1325, -20.914,
> -20.2144375, -15.52133333, -18.86783333, -18.54733333, -19.02683333,
> -15.12533333, -19.65833333, -17.56533333, -19.11833333, -19.17783333,
> -24.32233333, -24.03133333, -16.01333333, -23.81783333)
> location<-c(3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3)
> CONSUMER<- data.frame(d13C, d15N)
#get SOURCEtable:> mndN<-c(12.92019907, 12.42179762, 7.164152778)
> sdN<-c(0.924014709, 0.459319215, 2.106391697)
> mndC<-c(-14.43802315, -15.70665476, -31.72316667)
> sdC<-c(0.607857866, 0.429523631, 2.840500346)
> SOURCES<-data.frame(sp,mndN,sdN,mndC,sdC)
#get TEF-table:> sp<-c(1,2,3)
> meandN<-c(3.54,3.54,3.54)
> sddN<-c(0.001,0.001,0.001)
> meandC<-c(1.63,1.63,1.63)
> sddC<-c(0.01,0.01,0.01)
> TEF<-data.frame(sp,meandN,sddN,meandC,sddC)
#reading-in MODEL-SIAR:> model<-siarmcmcdirichletv4(CONSUMER, SOURCES, TEF,
> concdep=0,500000,50000)
##### after loading a lot, SESSION ABORTS here ####/
Who can help me?
--
View this message in context:
http://r.789695.n4.nabble.com/SIAR-problem-session-aborted-tp4649481.html
Sent from the R help mailing list archive at Nabble.com.
