phu but not phu

R Markdown

x <- rnorm(100)
x_na <- ifelse(x > 0, x, NA)
result <- mean(x_na, na.rm = TRUE)
cat("Avg = ", result)

## Avg =  0.6998367

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("\nRMSE = ", rmse)

## 
## RMSE =  1.52732