Error in Concentration Dependence

[robertwbuchkowski]

I am running a mixing model with a single isotope (15N) and concentration dependence. When I load the source data in using concentration dependence I get the error message:

Error in FUN(X[[i]], ...) :
'x' must be an array of at least two dimensions

That arises from the function lapply trying to calculate column means for a list with one dimensional components.

I edited the function by changing the lapply function from colMeans to mean and it appears to work well.

Does this make sense and is this a good fix?

Thanks,

Rob

CODE
if (conc_dep) {
CONC_names <- paste("Conc", mix$iso_names, sep = "")
if (sum(is.na(match(CONC_names, colnames(SOURCE)))) >
0) {
stop(paste("*** Error: Concentration dependence column names mislabeled.\n Should be 'Conc' + iso_names, e.g. 'Concd13C' if iso_names = 'd13C'.\n Please ensure Conc headings in source data file match iso_names\n in mix data file and try again. Alternatively, you may have set conc_dep=T\n by mistake.",
sep = ""))
}
CONC_iso_cols <- match(CONC_names, colnames(SOURCE))

if (length(CONC_iso_cols)>1){
conc <- do.call(rbind, lapply(split(SOURCE[, CONC_iso_cols], 
                                    list(SOURCE[, 1])), colMeans))
}else{
  
  conc <- do.call(rbind, lapply(split(SOURCE[, CONC_iso_cols], 
                                      list(SOURCE[, 1])), mean))
  colnames(conc) = "d15N"
}

}else conc <- NULL

[Heidi550]

I'm having the same issue as robertwbuchkowski and getting the error:

Error in FUN(X[[i]], ...) :
'x' must be an array of at least two dimensions

How do I access the code that robertwbuchkowski changed? Is there a filename with extension I could search for? Or is changing the code a challenging thing to do? If so, is there another workaround?