1条回答 默认 最新
- CSDN专家-HGJ 2021-10-27 16:24关注
参考下如下代码:
est<-function(x,y){ mux=mean(x) muy=mean(y) sx=0 sy=0 sxy=0 for (i in 1:length(x)){ sx=sx+(x[i]-mux)**2 sy=sy+(y[i]-muy)**2 sxy=sxy+(x[i]-mux)*(y[i]-muy) } res<-list(mu=c(mux,muy),s=matrix(data=list(sx,sxy,sxy,sy),dimnames=list(c('sx','sxy','sxy','sy')))) return(res) } x=c(150,149,131,143,133,140) y=c(40,42,35,38,30,43) est(x,y)
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 无用评论 打赏 举报
悬赏问题
- ¥20 ML307A在使用AT命令连接EMQX平台的MQTT时被拒绝
- ¥20 腾讯企业邮箱邮件可以恢复么
- ¥15 有人知道怎么将自己的迁移策略布到edgecloudsim上使用吗?
- ¥15 错误 LNK2001 无法解析的外部符号
- ¥50 安装pyaudiokits失败
- ¥15 计组这些题应该咋做呀
- ¥60 更换迈创SOL6M4AE卡的时候,驱动要重新装才能使用,怎么解决?
- ¥15 让node服务器有自动加载文件的功能
- ¥15 jmeter脚本回放有的是对的有的是错的
- ¥15 r语言蛋白组学相关问题