search for: mu_x

...ons and correlation sig_x <- 1 sig_y <- 1 rho_xy <- 0.0 # Covariance between X and Y sig_xy <- rho_xy * sig_x *sig_y # Covariance matrix Sigma_xy <- matrix(c(sig_x ^ 2, sig_xy, sig_xy, sig_y ^ 2), nrow = 2, ncol = 2) # Load the mvtnorm package library("mvtnorm") # Means mu_x <- 0 mu_y <- 0 # Simulate 1000 observations set.seed(12345) # for reproducibility xy_vals <- rmvnorm(1000, mean = c(mu_x, mu_y), sigma = Sigma_xy) # Have a look at the first observations head(xy_vals) # Create scatterplot plot(xy_vals[, 1], xy_vals[, 2], pch = 16, cex = 2, col = &quot...

...Island Pattern in India". To achieve the results I am applying a* two-dimensional Gaussian fit* on an LST raster of 1 km spatial resolution but I am facing two errors in the following code. library(raster) LST <- raster("D:/Celsius_Day/MOD_01.tif") gaussian2d <- function(x, y, mu_x, mu_y, sigma_x, sigma_y, amp) { exp(-((x - mu_x)^2/(2*sigma_x^2) + (y - mu_y)^2/(2*sigma_y^2)))*amp } #define a function for the sum of squared errors between the data and the Gaussian sse <- function(p) { mu_x <- p mu_y <- p sigma_x <- p sigma_y <- p amp <- p[5] fitted <- g...

I am trying to understand the meaning of the eps parameter of the equiv.test parameter of the equiv.test function (package equivUMP) The help file for equiv.test states that the parameter eps is "a single strictly positive number giving the equivalence limits" What is the scale of measurement of eps? It does not appear to be the same scale as the scale used for the two vectors that are