求解释语句

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

不太能理解这段语句，求大神解释一下
for (i in 1:10){
        mu <- Beta[1] + Beta[2] * md[i]
        yi <- rnorm(100, mean = mu, sd = summary(M0)$sigma)
        points(jitter(rep(md[i], 100)), jitter(yi), col = grey(0.5), pch = 16, cex = 1)
}
全部程序如下：
Fish<-read.table(file.choose(),sep=',',header=T,nrows=150)
Fish$MeanDepth<-Fish$MeanDepth / 1000
par(mfrow = c(1,1), mar = c(5,5,3,2))
plot(x = Fish$MeanDepth,
     y = Fish$TotAbund,
     xlab = "Mean Depth (km)",
     ylab = "Total abundance",
     cex.lab = 1.5,
     pch = 1,
     ylim = c(-300, 1200))
abline(M0, lwd = 5)
abline(h = 0, lty = 2)

range(Fish$MeanDepth)
md <- seq(0.804, 4.865, length = 10)

Beta <- coef(M0)
for (i in 1:10){
        mu <- Beta[1] + Beta[2] * md[i]
        yi <- rnorm(100, mean = mu, sd = summary(M0)$sigma)
        points(jitter(rep(md[i], 100)), jitter(yi), col = grey(0.5), pch = 16, cex = 1)
}

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：For points sigma 程序

前面没调好。。。重新把全部的代码发一下

1. par(mfrow = c(1,1), mar = c(5,5,3,2))
   
2. plot(x = Fish$MeanDepth,
   
3.      y = Fish$TotAbund,
   
4.      xlab = "Mean Depth (km)",
   
5.      ylab = "Total abundance",
   
6.      cex.lab = 1.5,
   
7.      pch = 1,
   
8.      ylim = c(-300, 1200))
   
9. abline(M0, lwd = 5)
   
10. abline(h = 0, lty = 2)
    
11. 
    
12. range(Fish$MeanDepth)
    
13. md <- seq(0.804, 4.865, length = 10)
    
14. 
    
15. Beta <- coef(M0)
    
16. for (i in 1:10){
    
17.         mu <- Beta[1] + Beta[2] * md[i]
    
18.         yi <- rnorm(100, mean = mu, sd = summary(M0)$sigma)
    
19.         points(jitter(rep(md[i], 100)), jitter(yi), col = grey(0.5), pch = 16, cex = 1)
    
20. }

复制代码