那加权最小二乘估计要在SAS中实现,是否只能自己编程实现?
- 于 概率论的一问题
不对哦,答案是1/2,问题其实我已经解决了,看到有些书上的解答很简单,而我觉得需要交代的东西很多,会不会是我想多了,简单的某点没有想到
我的解答是:由P(A并B)小于等于P(A并B并C)就可以得到x<=1/2,再找个例子说明一下1/2是可以取到的
- 于 概率论的一问题
设A,B,C两两独立,且ABC为不可能事件
如果P(A)=P(B)=P(C)=x,求x的最大值
- 刚才拜读了那两篇文章,我支持老先生!!!
- hehe
- look
- 我想应该是这样的吧:
1,方差分析的原假设是检验均值相等,这是肯定的
2,方差分析模型中是有一些基本假设的,一般会去假设方差满足齐性,也就是各水平下的方差相等
3,你说教科书上公认F检验是检验方差是否相等的,这是你认识不足,可能你刚开始学概率统计时最早看到的F检验是检验正态总体的方差是否相等,但事实上是F检验有很多地方都有。 所谓F检验是指我们用到的统计量在原假设成立的时候服从F分布而已,不是说F检验就是检验方差相等的。
事实你在教科书上最早接触到正态总体的均值相等的假设的时候,可能显示U检验,然后是T检验,事实也可以做F检验,事实你去看看t分布和F分布之间的关系,你就明白了 - 数学中不是同样有很多概念不是客观存在的吗,比如说四维的向量空间难道是客观存在的?
很多概念不是因为实际应用中的需要才给出来的吗,那么总体价格指数是不是也是为了说明某一时期的价格水平而产生的呢,您的意思是现在的价格指数不能表示价格水平还是它根本不该存在? - 标准化回归系数应该是指2楼说的那样吧
- 于 统计人的悲哀人家用平均数有什么问题,国民都有了基本统计素质以后,他那个数据碍得了任何事情?
你或者说你从那个数据里看不到你真正想看的东西,这也是可以理解的吧,哪个国家不会藏点掖点啥的,该
公布的时候自然会告诉你,有时候要从对方的角度换位想想啊 - 那么综合国力是个啥?请老先生指教一下后辈
- 看看
- 应该是个高等统计的问题了吧,最优检验记得在本科教学的时候都是带星号的内容,现在的本科教材上很多干脆就没有这部分的内容,到了研究生的时候这部分内容好像才详细的拿出来讨论
可能是我没有清楚你的意思,用p值的频率性质来说明检验的好坏应该是可以的,但是从一个p值肯定是无法说明检验的好坏的,所谓好的检验就是犯两种错误的概率要低,一般都是控制犯第一类错的概率的前提下,尽可能地使犯第二类错误的概率低,
我在猜想你讲到的论文中首先是有很多样本,其次样本的来源是从已知分布中模拟得出的。
总之我还是不明白你的意思 - 你这个问题似乎很难回答哦,看到你的问题似乎你数理统计知识已经到高等统计的那块了哦
可是你说通过p值来判断检验的好坏,我似乎从没有听说过,p值一般是拿来下结论的,也就是是该拒绝原假设还是接受原假设
讲到检验的好坏一般是用势(power),具体理论建议你看看相关的高等统计的书,这里很难一下自讲清楚吧,呵呵,反正我讲不清 - 于 R的数据导入x<-matrix(scan("e:\\exp1data.txt"),ncol=1,byrow=T)
y<-x[31:60]/*将x的后30个数据赋给变量y,前30个给变量x,由于这里的矩阵x和变量x同名,所以要注意先给y赋值*/
x<-x[1:30] - 于 R的数据导入[quote]引用楼主ypchen于2009-07-23 14:17发表的 R的数据导入 :
逗号分隔的文本文件
<br /> # first row contains variable names, comma is separator <br /> # assign the variable id to row names <br /> # note the / instead of \ on mswindows systems <br /> ....... [url=http://cos.name/bbs/job.php?action=topost&tid=15900&pid=tpc][/url]<br /> [/quote]<br /> 弱弱地问一下,逗号分隔的文本文件是不是就是将excel文件另存为扩展名为csv的那个文件?
- 觉得自己在这个论坛得到不少帮助,自己似乎没有啥贡献,有些惭愧。
最近自己在学习R,主要是线性模型方面的,学了点皮毛,在这里把自己学到的一点点东西在这里向大家报告一下,或许对象我这样的初学者有些帮助,不对的地方也好请高手们指正:
//*从excel表中读取数据,先将数据文件存成csv文件,再用read.csv去读取*/
data<-read.csv(file="e:\\work\\2009下半年\\modern statistics-2009new\\ding learning note\\linear model\\例2.2.csv")
/*将data的第二列数值赋给变量x,第三列赋给变量y*/
x<-data[,2]
y<-data[,3]
/*做散点图*/
plot(x,y)
/*线性回归用的是函数lm()*/
z<-lm(y~x)
summary(z)
/*以下代码是没有截距的回归*/
z<-lm(y~x-1)
z<-lm(y~0+x)
library(MASS)
z<-lm(y~x)
summary(z)
plot(x,y)
plot(fitted(z),studres(z))/*残差分析,拟合值与学生化残差图*/
qqnorm(studres(z))/*对残差进行正态性检验*/
qqline(studres(z))
这里调用library(mass)我还是不太明白,估计是将回归的一些结果拿出来用吧,因为之前我没有这个的时候,做不了plot(fitted(z),studres(z))
这是一元线性回归的一点点内容,多元的我正在学,有机会在向大家报告 - [quote]引用第5楼yuan^2于2009-07-24 20:52发表的 :
? read.csv
便能看到所有相关详细信息
header 的默认值是FALSE
header=TRUE就是保留你设置的列名,例如x,y
[/quote]
好像不是噢,我没有设置这个参数值,第一行的文字也出来了哦
data<-read.csv(file="e:\\work\\例2.2.csv")
data
这是我的代码,得出来的是
年份 人均国民收入 人均消费金额
1 1980 460 234.75
2 1981 489 259.26
3 1982 525 280.58
4 1983 580 305.97
5 1984 692 347.15
6 1985 853 433.53
7 1986 956 481.36
8 1987 1104 545.40
9 1988 1355 687.51
10 1989 1512 756.27
11 1990 1634 797.08
12 1991 1879 890.66
13 1992 2287 1063.39
14 1993 2939 1323.22
15 1994 3923 1736.32
16 1995 4854 2224.59
17 1996 5576 2627.06
18 1997 6053 2819.36
19 1998 6392 2958.18 - 太感谢了
按照你给的方法真的就可以了噢,为什么呢?
如果我再excel表格的第一行写上文字描述每列的数据,读出数据后,R在作图时会忽略第一行吗