Advertisement
Guest User

Untitled

a guest
Dec 14th, 2019
116
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.31 KB | None | 0 0
  1. df <- read.delim('ais.txt')
  2. sapply(df,class)
  3. summary(df)
  4. cor(df[-c(1,2)])
  5. plot(df[-c(1,2)], col = df$Sex)
  6. x_v = df$Wt
  7. y_v = df$Ht
  8. B1 = sum((y_v-mean(y_v))*(x_v-mean(x_v)))/sum((x_v-mean(x_v))^2)
  9. B0 = mean(y_v)-B1*mean(x_v)
  10. plot(df$Wt, df$Ht)
  11. abline(B0,B1)
  12.  
  13. y_pred <- B0 + B1x
  14. points(x_v, y_pred, col = "red")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement