接受各种R语言难题

上次咨询的Cochran-Armitage trend test已粗略解决，感谢大神  为保证可靠性，需要等距赋值。
参考文献:
  1.何春拉,潘建红,陈,靖等.Cochran—Armitage 趋势检验不同赋值的模拟研究[J].中国卫生统计,2012,29(6):843-845

#赋值观察矩阵

Age <-matrix(c(108,373,494,1283,214,387), byrow=TRUE, nrow=3,dimnames = list(age=c('0-20','20-40','40-60'), serious=c('y','n')))

#方法一：调用coin包，使用chisq_test函数

require('coin')

Age1 <- as.table(Age)

chisq_test(Age,scores=list(age=c(2,1,0)))

#方法二：调用Desctools包，使用CochranArmitageTest函数

require('DescTools')

CochranArmitageTest(Age)

大神，您好，我在用R软件进行面板进行“固定效应变截距”模型的估计时，R并不给出截距项，但给出了每个虚拟变量的的估计值；那这言外之意是不是得到的结果即是各部门的截距项呢，也不存在公共截距项，即不存在提取公共截距项这一说法呢。我在学习《计量经济学实验教材——基于EViews和R软件的应用》这本书的面板模型那章的时候，书上在讲到“固定效应变截距模型的估计”时，说summary（plm的模型）之后没有给出模型的公共截距，公共截距存在模型的alpha中，它提取公共截距的代码为”plm的模型$alpha“，但我在运行中，并不能得到该公共截距，在help(plm)中，也没有查到相关的信息。那我刚刚的认为是否正确，书上是否存在错误，抑或是R软件版本不同的原因。盼解答，多谢~

大神，我想画个两条正态密度曲线图表示误判概率的，如附件图，求助无人回应，请大神指点，非常感谢！

gtliu2016 发表于 2016-12-14 09:47
大神，我想画个两条正态密度曲线图表示误判概率的，如附件图，求助无人回应，请大神指点，非常感谢！

用常用的plot或者ggplot画出来就行。

gtliu2016 发表于 2016-12-14 09:47
大神，我想画个两条正态密度曲线图表示误判概率的，如附件图，求助无人回应，请大神指点，非常感谢！

可以参考：
https://www.r-bloggers.com/creating-shaded-areas-in-r/

gtliu2016 发表于 2016-12-14 09:47
大神，我想画个两条正态密度曲线图表示误判概率的，如附件图，求助无人回应，请大神指点，非常感谢！

可以参考：
https://www.r-bloggers.com/creating-shaded-areas-in-r/

tju_gyj 发表于 2016-12-14 01:08
大神，您好，我在用R软件进行面板进行“固定效应变截距”模型的估计时，R并不给出截距项，但给出了每个虚拟 ...

回归结果中，对某一个虚拟变量值的估计应该是没有的，那个就是基准，默认为0

cbw1243 发表于 2016-12-14 14:13
回归结果中，对某一个虚拟变量值的估计应该是没有的，那个就是基准，默认为0

谢谢，明白了。下午我自己应该也是确定了，R在固定效应模型中，并没有公共截距项，它所给出的就是每个部门的固定效应。和我参考的书出现冲突，估计是因为R的版本不一样吧。另又有一个问题请教一下大神，R在随机效应模型面板估计时，R用什么函数命令能提取该随机效应估计值呢？多谢~

程序最后结果提示：non-numeric argument to mathematical function什么意思

用R做一个模型
因变量Y是定序的，1到5分的量表
自变量，最主要的是A，A影响B和C，对因变量起作用
Y是个体层面数据，ABC是城市层面数据

求大神回复！