[问题] 想用R写一个统计回归分析的程式

楼主: eegame (yu)   2015-11-20 14:24:02
初学者(两个月) R是第一个程式语言
不知道各位前辈有没做过类似的程式
!!!假设!!!
想用R统计10条绳子的长度
再根据统计特性的率随机产生10组数据
程式如下
G <- 10 #10组数据
rowname <- c(1:G) #Y轴名称
colname <- c("c1","c2","c3","c4") #X轴名称
F <- matrix(,nrow=G,ncol=4,dimnames=list(rowname,colname))?
#开一个10*4的矩阵并以rowname和colname
F[1,"c1"] <- 19
F[2,"c1"] <- 59
F[3,"c1"] <- 48
F[4,"c1"] <- 23
F[5,"c1"] <- 31
F[6,"c1"] <- 48
F[7,"c1"] <- 56
F[8,"c1"] <- 60
F[9,"c1"] <- 74
F[10,"c1"] <- 83
#在矩阵c1的位置放入这10条绳子的调查数据
#假设1符合常态分布
F[,"c2"] <- rnorm(10,10,90)?
#在矩阵c2的位置放入根据常态分布的机率随机产生一组新的数据介于10~90之间的10笔数

#假设2符合均匀分布
F[,"c3"] <- dunif(10,10,90)
#在矩阵c3的位置放入根据均匀分布的机率随机产生一组新的数据介于10~90之间的10笔数

#假设3符合Poisson distribution分布
F[,"c4"] <- dpois(10, lambda=1)
#在矩阵c4的位置放入根据Poisson distribution分布的机率随机产生一组新的数据介于1
0~90之间的10笔数据
#结果如下
F
c1 c2 c3 c4
1? 19 -70.255568 0.0125 1.013777e-07
2? 59 9.346476 0.0125 1.013777e-07
3? 48? -2.139539 0.0125 1.013777e-07
4? 23 -33.231206 0.0125 1.013777e-07
5? 31? 72.630937 0.0125 1.013777e-07
6? 48? -6.032507 0.0125 1.013777e-07
7? 56? -8.420308 0.0125 1.013777e-07
8? 60 -69.542557 0.0125 1.013777e-07
9? 74? 12.044775 0.0125 1.013777e-07
10 83 -13.474063 0.0125 1.013777e-07
问题1.我函数的指令是不是给的不对导致c2 c3 c4的结果错误?
问题2.关于随机分布的应用 如何做到根据统计特性产生结果?
问题3.关于数据统计特性如何检定比较符合 常态分布 均匀分布 负指数分布等等的分布
特性? 是用卡方检定 T-test检定? 还是用数值回归???
这两个月没有看到类似的文章
麻烦各位前辈指点迷津
拜托 拜托
作者: celestialgod (天)   2015-11-20 14:28:00
1. 你知道各分配参数意义吗?2. 你知道第一个问题就知道怎样根据你观察的资料去产生适当的随机数3. 你的目标是检定什么?十条绳子长度符合某分配?4. 你学过适合度检定吗?另外,我很好奇这跟回归(标题)有什么关系啊?5. 你产生随机数的目的是什么?6. 会不会你等一下又要删文啦?
作者: allen1985 (我要低调 拯救形象)   2015-11-21 17:03:00
所以你的问题跟回归没有关系?最合适的分布特性? 你是说资料最符合哪种分配吗?你问题中的自变量 应变量 是?建议你简化一下你的题目 一步一步问 从头开始大概知道你的问题是什么了 你有听过qq-plot吗对我来说 你用ggplot2 还是什么function画图 不重要重要的是你想要画什么图 你的图想表达什么qq-plot 是用来检测资料是否为常态分配的一种方法而你现在最欠缺的是 你得先把统计方法学好其次才是用R把结果写出来
作者: celestialgod (天)   2015-11-22 14:25:00
推楼上,重点不是R,他只是工具,你怎么呈现你的想法才是重点

Links booklink

Contact Us: admin [ a t ] ucptt.com