SpaceVacation <- function(n=100) {T1 = 90 + rnorm(2*n)*rexp(2*n,0.5) T2 = 50 + rnorm(4*n,0,8) TempY = c(T1,T2) TempX =rnorm(6*n,63.5,7) Vacation = data.frame(TempX,TempY) return(Vacation)}