x <- rnorm(100)
x_na <- ifelse(x > 0, x, NA)
result <- mean(x_na, na.rm = TRUE)
cat("Avg =", result)
## Avg = 0.8759744
t <- rnorm(100)
y <- rnorm(100)
df <- data.frame(t, y)
er <- df$t - df$y
sq <- er^2
mse <- mean(sq)
rmse <- sqrt(mse)
cat("RMSE =", rmse)
## RMSE = 1.311479