[R] SIAR: problem> session aborted

Jeltje ecoqueen at live.nl
Wed Nov 14 14:25:31 CET 2012


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.



More information about the R-help mailing list