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. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
