全部版块 我的主页
› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛
接受各种R语言难题
fanfuqiu
2016-12-7 15:53:57
cbw1243 发表于 2016-12-6 22:41
systemfit 包是我用过比较好的包。
我是做空间面板  一直提醒错误：Error in mat2listw(mydata1) : x is not a matrix。
说我的矩阵不是对称的，但我检查了，矩阵没问题？
这是怎磨回事？
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

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

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

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)
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

fanfuqiu
2016-12-8 10:31:20
cbw1243 发表于 2016-12-7 22:49
is.matrix(x)
谢啦  
还有 楼主 做动态空间面板模型，R的程序包是什么？
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

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 为其余两个虚拟变量
      }
    }
  }
}

二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-8 12:39:30
fanfuqiu 发表于 2016-12-8 10:31
谢啦  
还有 楼主 做动态空间面板模型，R的程序包是什么？
plm 是经常用的针对panel的包
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-8 12:42:01
myohmy88 发表于 2016-12-7 00:15
这个函数不连续，所以不可导，所以海塞矩阵不存在吧
如果不连续的话用区间内的MLE. 然后进行区间比较
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-8 12:43:26
fanghang716 发表于 2016-12-8 12:39
你好，再回归时需要对一个two-level CES函数进行回归
用R回归时，涉及到四个参数的设置，a,b区间在（0,1） ...
在R中尽量不要用循环。
可以试试tapply
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

fanfuqiu
2016-12-8 15:43:47
cbw1243 发表于 2016-12-8 12:39
plm 是经常用的针对panel的包
plm是针对普通面板的程序包
splm是针对空间静态面板的程序包
我想问的是 针对空间动态面板的程序包是什么？
谢谢啦亲
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

Licko123456
2016-12-9 13:21:28
今天上矩阵代数，我朋友电脑的Rstudio运行程序只运行第一行代码，也没有报错，就是无论运行什么程序、函数都只运行第一句，我就帮她弄，没搞定，不知道原因是什么，后来我们为此吵架了，特来此求助！希望有人能够解答一下怎么处理解决
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-9 23:45:45
Licko123456 发表于 2016-12-9 13:21
今天上矩阵代数，我朋友电脑的Rstudio运行程序只运行第一行代码，也没有报错，就是无论运行什么程序、函数都 ...
这种情况我还是第一次听说...你确定你运行代码的方式是对的吗？
试试写两行代码，然后把结果截过来看一下
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

zyx1991125
2016-12-10 12:07:22
谢谢哈~~
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

风之将尽
2016-12-10 21:05:56
用wordcloud画词云图结果好几个词没输出出来，显示
Warning messages:
1: In strwidth(words[i], cex = size[i], ...) :
  Windows字体数据库里没有这样的字体系列
2: In wordcloud(t词词,t词频, random.order = F, random.color = F, colors = mycolors,  :
  供给侧改革 could not be fit on page. It will not be plotted.
请问这是怎么回事呢?很奇怪为啥报错了，因为单独输出供给侧、改革都是没问题的，就是一起输出时出问题了。
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ccj3333
2016-12-10 23:07:17
我怎么样来进行大数据的操作（包括导入的速度，操作的函数与包的选择）
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

冰小二
2016-12-11 20:48:43
已经过了暑假，不知大神还登不登论坛，想问下大神空间贝叶斯生存模型的R程序实现难吗？作为程序小白感受不到难易（捂脸），要自己编代码吗？有没有现成的程序，只需要代入数据直接就跑出来了。我看不懂它的程序包，不知道怎么代入空间生存数据，还有模型的前奏DIC和PD的计算~~小硕毕业论文所用，大神能回复的话感激不尽
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

feeloob
2016-12-11 23:32:49
大神，求求你，我的R语言 包安装不起，帮助语言？mean() 没法用啊！！
打开rstudio 先出现这个 Error in tools::startDynamicHelp() : unable to create socket
输入？mean() 会出现错误代码  starting httpd help server ...Error in startDynamicHelp(TRUE) : unable to create socket
R现在根本用不了，请大神教教怎么办啊！！！
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-12 04:34:32
ccj3333 发表于 2016-12-10 23:07
我怎么样来进行大数据的操作（包括导入的速度，操作的函数与包的选择）
大数据处理并不是简单能总结的。但是，基础的处理方式用dplyr这一类型的包。在函数的选择上尽量选择 base里的函数，或者外接Rcpp的函数
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-12 04:36:45
冰小二 发表于 2016-12-11 20:48
已经过了暑假，不知大神还登不登论坛，想问下大神空间贝叶斯生存模型的R程序实现难吗？作为程序小白感受不到 ...
贝叶斯模型在R中的应用还不太广泛，这就意味着在大部分的情况下，都需要自己coding
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-12 04:39:04
feeloob 发表于 2016-12-11 23:32
大神，求求你，我的R语言 包安装不起，帮助语言？mean() 没法用啊！！
打开rstudio 先出现这个 Error in t ...
在Rstudio中改一下设置就行。参考
http://stackoverflow.com/questio ... le-to-create-socket
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

冰小二
2016-12-12 10:52:42
cbw1243 发表于 2016-12-12 04:36
贝叶斯模型在R中的应用还不太广泛，这就意味着在大部分的情况下，都需要自己coding
唉，蓝瘦香菇~~还是谢谢大神~
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

松石空月
2016-12-13 20:41:45
上次咨询的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)




二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

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

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

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

原图尺寸 11.7 KB

误判概率图.jpg

二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-14 14:09:24
gtliu2016 发表于 2016-12-14 09:47
大神，我想画个两条正态密度曲线图表示误判概率的，如附件图，求助无人回应，请大神指点，非常感谢！
用常用的plot或者ggplot画出来就行。
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-14 14:10:59
gtliu2016 发表于 2016-12-14 09:47
大神，我想画个两条正态密度曲线图表示误判概率的，如附件图，求助无人回应，请大神指点，非常感谢！
可以参考：
https://www.r-bloggers.com/creating-shaded-areas-in-r/
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-14 14:11:25
gtliu2016 发表于 2016-12-14 09:47
大神，我想画个两条正态密度曲线图表示误判概率的，如附件图，求助无人回应，请大神指点，非常感谢！
可以参考：
https://www.r-bloggers.com/creating-shaded-areas-in-r/
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

cbw1243
2016-12-14 14:13:36
tju_gyj 发表于 2016-12-14 01:08
大神，您好，我在用R软件进行面板进行“固定效应变截距”模型的估计时，R并不给出截距项，但给出了每个虚拟 ...
回归结果中，对某一个虚拟变量值的估计应该是没有的，那个就是基准，默认为0
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

tju_gyj
2016-12-14 17:43:11
cbw1243 发表于 2016-12-14 14:13
回归结果中，对某一个虚拟变量值的估计应该是没有的，那个就是基准，默认为0
谢谢，明白了。下午我自己应该也是确定了，R在固定效应模型中，并没有公共截距项，它所给出的就是每个部门的固定效应。和我参考的书出现冲突，估计是因为R的版本不一样吧。另又有一个问题请教一下大神，R在随机效应模型面板估计时，R用什么函数命令能提取该随机效应估计值呢？多谢~
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

阿森纳冠军
2016-12-15 10:40:52
程序最后结果提示：non-numeric argument to mathematical function什么意思
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

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

求大神回复！
二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

首页上一页下一页跳至第页
栏目导航
热门文章
推荐文章

说点什么

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群