similar to: Errors in SIAR

Help! model1 <- siarmcmcdirichletv4(data, source, tef, concdep=0, 500000, 50000) Error in matrix(1, ncol = (numsources + numiso) * numgroups, nrow = (siardata$iterations - : invalid 'ncol' value (too large or NA) In addition: Warning message: In Ops.ordered((numsources + numiso), numgroups) : '*' is not meaningful for ordered factors Previous message, but more problems:

Hi, my name is Julio, I'm trying to run a script, specifically with the siar package on R (I have attached the script) but the program collapses on during the run. The error message says: *** caught segfault *** address 0xb31941cc, cause 'memory not mapped' Traceback: 1: .C("siarmcmcv4", as.integer(numdata), as.integer(numsources), as.integer(numiso),

Hi there, I am trying to run an MCMC on stable isotope data from certain organisms to determine their dietary habits in the package SIAR. I have prepared my data according to Inger, R., Jackson, A., Parnell, A., Bearhop, S. : SIAR V4 (Stable Isotope Analysis in R) an Ecologist's Guide (also better known as 'SIAR for dummies'), which is as follows: - consumer data table (columns

Hello, Beginner, sorry if this is wasting anyone's time, but have been working on this for a couple of days now, think it should have take a few hours! The Problem: Error in Summary.factor(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 4L, : max not meaningful for factors I have tried to re-arrange and check formula. I am working in SIAR, and cannot get model1 to run because of the error

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

Hello, I am trying to run the demoscript for the package "siar". However, using the example script and example data supplied I cannot get 1) the model to run unless I specify iterations =500 and burnin = 500 and 2) if the model does work without crashing R then the function "siarhistograms" does not work . Can anyone point out where I am going wrong? >model1 <-

Hi, I have a problem while loading the "siar" program in R. When I am loading siar, system does not load convexhull. On the screen I have seen such writings. The following object(s) are masked from ‘package:spatstat’: convexhull How can I load the convexhull, how can I unmask from this package? I will be appreciated if you give advice about this. Best Sukran

Hello, I recently bought a new Macbook Pro and installed R via a disk image from an older computer - it seems to work fine. I then installed the package "siar", tried to load it via the library command and received the following messages (note the "mvtnorm" package did not appear to load correctly): trying URL

Hi. I have the following code in a loop. It splits a vector into subvectors of equal size. But if the size of the original vector is not an exact multiple of the desired subvector size, then the first few subvectors have one more element than the last few. I know that the cut function could be used to determine where to break up the vector, but it doesn't seem to provide control over

On Nov 12, 2011 at 8:29pm Alex wrote: > Has anyone used SIAR package add on? I posted a reply to an earlier question from you on this subject. See http://r.789695.n4.nabble.com/Errors-in-SIAR-td4029804.html. In it I note that there are problems with the function from siar (not SIAR) you are using, but that this may not be "your" problem, that the function calls for matrices (you

Hello! I have larger and a smaller data frame with 1 factor in each - it's the same factor: large.frame<-data.frame(myfactor=LETTERS[1:10]) small.frame<-data.frame(myfactor=LETTERS[c(9,7,5,3,1)]) levels(large.frame$myfactor) levels(small.frame$myfactor) table(large.frame$myfactor) table(small.frame$myfactor) myfactor has 10 levels in large.frame and 5 levels in small.frame. All 5

Hi, Does tapply change the order when applied on a factor? Below is the code I tried. > mylevels<-c("IN0020020155","IN0019800021","IN0020020064") >

I have a data frame with some factor columns. I want to drop the rows with rare factor values (and remove the factor values from the factors). E.g., frame$MyFactor takes values A 1,000 times, B 2,000 times, C 30 times and D 4 times. I want to remove all rows which assume rare values (<1%), i.e., C and D. i.e., frame <- frame[[! (frame$MyFactor %in% c("A","B"))]] except

Dear R help group, I am teaching myself linear mixed models with missing data since I would like to analyze a stats design with these kind of models. The textbook example is for the procedure "proc MIXED" in SAS, but I would like to know if there is an equivalent in R. This example only includes two time-measurements across subjects (a t-test "with missing values"), but I

Saludos: Hace unos días me encontré con un problema en la función readSIAR del paquete solaR; pensé en actualizar pero no lo hice en el momento. El caso es que acabo de actualizar el paquete solaR -y el resto de paquetes también- y veo que el error no se va. ------ Corte desde la consola -------- R version 2.15.1 (2012-06-22) -- "Roasted Marshmallows" Copyright (C) 2012 The R

Hi folks, I am working with repeated measures data and I ran into issues where the paired t-test results did not match those obtained by employing glht() contrasts on a lme model. While the lme model itself appears to be fine, there seems to be some discrepancy with using glht() on the lme model (unless I am missing something here). I was wondering if someone could help identify the issue. On

Hello, I have a list of data with multiple elements, and each element in the list has multiple variables in it. Here's an example: ### Make the fake data dv <- c(1,3,4,2,2,3,2,5,6,3,4,4,3,5,6) subject <- factor(c("s1","s1","s1","s2","s2","s2","s3","s3","s3",

Hi, I searched for this in the mailing list, but found no results. I have a large dataframe ( dim(mydata)= 1297059 16, object.size(mydata= 145280576) ) , and I want to perform some calculations which can be done by a factor's levels, say, mydata$myfactor. So what I want is to split this dataframe into nlevels(mydata$myfactor) = 80 levels. But I must do this efficiently, that is, I

Hi I'm new to R and would like to implement a SAS-like macro variable in R. What I'd like to do is take the simple R code below and change the "=TEF" to different letters to refer to different companies' data for download. # DOWNLOADS FILES FROM YAHOO INTERNET

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20071208/5409f1a7/attachment.pl