十五年R语言统计分析经验，开贴答疑，回答统计学问题

楼主您好。我将检验正态分布的函数和for函数连用后，为什么会报错。代码如下，谢谢您了。
c1<-1:4
c2<-1:5
c3<-1:10
c4<-1:10
for (i in c(c1,c2,c3,c4)) {
shapiro.test(i)
}

楼主你好，想问一下我用rstudio画的图不太清晰是什么原因，该怎么处理

楼主你好，我是初学R语言，下面的代码请教是有什么问题呢
legend("topleft",inset=.05,title="Drug Type",c("A","B")
lty=c(1, 2),pch=c(15, 17),col=c("red","blue"))
par(opar)
报错，错误: unexpected symbol in:
"legend("topleft",inset=.05,title="Drug Type",c("A","B")
lty"

有个略微复杂一点的问题想请教一下。
我现在有两组数据，一组是高程(dem)，一组是我们测的一个值(hq15)。现在想做两件事，一个是用高程做X轴看散点分布图，并做出来趋势线，这个基本实现了。
另外一个是找到趋势线上边的某些突变点位。也就是说目前趋势线是有起伏波动的，现在想找到的就是在某个高程下，测量值明显发生改变的值，也就是由于高程产生的突变点。并且标记在图上。目前就卡在这了，不知道用哪个统计包能做出来？

图大致就是这个样子。比如图上在600左右和900左右，还有1400左右都存在较明显的拐点，我想提取出来对应的DEM值，并在图上用直线标注在X轴。

大神，怎么根据R语言输出的coefficients: 写出ARIMA乘法模型表达式啊？