[问答] 接受各种R语言难题 [推广有奖]

421楼

fanfuqiu 发表于 2016-12-7 15:53:57

cbw1243 发表于 2016-12-6 22:41
systemfit 包是我用过比较好的包。

我是做空间面板一直提醒错误：Error in mat2listw(mydata1) : x is not a matrix。
说我的矩阵不是对称的，但我检查了，矩阵没问题？
这是怎磨回事？

422楼

wenny1124 发表于 2016-12-7 17:21:55

您好，我用因子分析来做几个城市的排名，有几个问题想咨询一下：
1、用的fa函数，我用的原始数据，然后用fa.score得出几个因子得分，用方差贡献率作为权重，来得出最终的得分来，这个思路对不？
2、我按照上面思路做的，但是发现因子载荷系数好像不太对，比如森林覆盖率越大越好，但是对应的因子载荷却是负的，这个是对的吗？我有点弄不清这个因子载荷跟因子得分的关系了。
万分感谢！

423楼

cbw1243 发表于 2016-12-7 22:49:01

fanfuqiu 发表于 2016-12-7 15:53
我是做空间面板一直提醒错误：Error in mat2listw(mydata1) : x is not a matrix。
说我的矩阵不是对称 ...

is.matrix(x)

424楼

fanfuqiu 发表于 2016-12-8 10:31:20

cbw1243 发表于 2016-12-7 22:49
is.matrix(x)

谢啦
还有楼主做动态空间面板模型，R的程序包是什么？

425楼

fanghang716

发表于 2016-12-8 12:39:11

你好，再回归时需要对一个two-level CES函数进行回归捕获.JPG

用R回归时，涉及到四个参数的设置，a,b区间在（0,1），\[\theta\] 和\[\rho\] 区间设置为（-20, 0.99）。要求出Y对双CES表达式的回归结果，问题时=是，我R只是初学者，只会嵌套循环，运行一次都要等很久。怎么样才能避免使用循环

r2 = NULL
r2 = array(rep(0, len_j1*len_i1*len_j2*len_i2),dim= c(len_j1, len_i1, len_j2, len_i2))
str(r2)
for (i1_index in 1:len_i) # j1表示theta, j2表示rho, i1表示b, i2表示a
{
  for (j1_index in 1:len_j)
  {
for (i2_index in 1:len_i)
{
   for (j2_index in 1:len_j)
   {
r2[j1_index, i1_index, j2_index, i2_index] = summary(lm(GDP.ln ~ ces.ln[,j1_index, i1_index, j2_index, i2_index] + KH1H2$d + tech1))$r.squared    #KH1H2$d 和 tech1 为其余两个虚拟变量
   }
}
  }
}

我爱人大论坛