Displaying 3 results from an estimated 3 matches for "thinby".
2010 Sep 07
1
problems with siar package
...pped'
Traceback:
1: .C("siarmcmcv4", as.integer(numdata), as.integer(numsources),
as.integer(numiso), as.integer(numgroups), as.integer(startgroup),
as.integer(endgroup), as.integer(siardata$iterations),
as.integer(siardata$burnin), as.integer(siardata$howmany),
as.integer(siardata$thinby), as.double(prior), as.data.frame(data2),
as.data.frame(concdepdata), as.data.frame(sourcedata),
as.data.frame(correctionsdata), as.data.frame(parameters))
2: siarmcmcdirichletv4(target, source, fraction, iterations = 5e+05, burnin
= 1e+05, howmany = 2000, thinby = 10)
Possible actions:
1: abor...
2013 Oct 31
0
SIAR problem with model running
...5,
50000): NAs in foreign function call (arg 5)'
there aren't any NAs in my data set, so I assume that NAs are being created
by argument 5 in my code? I've tried to change the number for the
'burning', but then I get another error warning:
'Error in iterations, burning or thinby: (iterations-burnin)/thinby must be
an integer.'
What can I do to get the MCMC running?
thanks a lot for your help.
Best,
Claudia
[[alternative HTML version deleted]]
2011 Nov 10
4
Error in Summary.factor(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 4L, : max not meaningful for factors
...0.4 1.3 3.4 0.99
2 Algal Paste 0.4 1.3 3.4 0.99
3 Amphipoda 0.4 1.3 3.4 0.99
> model1 <- siarmcmcdirichletv4(data, source, tef, concdep = 0,
> iterations=500000, burnin=50000, howmany=50000, thinby=15, prior = rep(1,
> nrow(sources)), siardata=list(SHOULDRUN=FALSE))
Error in Summary.factor(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 4L, :
max not meaningful for factors*
> source
sources Meand13C SDd13C Meand15N SDd15N
1 2 2 2...