D
​ear community,

This is to kindly request for your help.
 I have an error from regression of values in one stack with another
# s1 and s2 have 720 layers; coefficients[2] is the slope


### script
rstack1 <- stack(s1,s2)
s <- stack('D:/Correlation/rstack.tif')
fun <- function(x) { lm(x[1:360] ~ x[361:720])$coefficients[2] }
x1 <- calc(s, fun)

I get the following error:

Error in .calcTest(x[1:5], fun, na.rm, forcefun, forceapply) : cannot
use this function

Many thanks for your help in sloving this problem

John

​

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.

Reply via email to