楼主: cbw1243
44312 504

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

371
2000918918 发表于 2016-11-29 14:17:54
建立空间滞后模型SLM后用predict做预测,遇到以下问题:
1.如何设定predict的参数使预测的值与SLM$fitted.values对应值相等?
data(oldcol)
lw <- nb2listw(COL.nb)
SLM <- lagsarlm(CRIME ~ INC + HOVAL, data=COL.OLD, lw,  type="mixed")
2.我想用predict函数预测空间上一个点的值,怎么预测?

372
cbw1243 发表于 2016-11-29 22:45:21
fanfuqiu 发表于 2016-11-25 11:24
大神,R语言可以做空间面板的联立方程吗? gs3sls
参考systemfit包

373
cbw1243 发表于 2016-11-29 22:48:44
zyx1991125 发表于 2016-11-22 22:15
楼主,这好像是做DEA和SFA的,我意思是做投入产出分析,列昂惕夫矩阵什么的。。求指导
可以自己写代码进行Leontief 矩阵运算啊。很简单

374
chenliangqiang 在职认证  发表于 2016-11-29 23:05:47
请教一下,我用vegan包做的CCA分析。如何改变图中字体大小、颜色等,线条颜色粗细等,以及在拉丁文前面标上特定符号呢?

CCA分析结果.jpeg (26.14 KB)

CCA分析结果.jpeg

375
3210sq 发表于 2016-11-30 00:11:18
小问题,小骚扰。

例如c(4 3 6 5 1 3 5 1 8 4 3 1 8 1 1 8 4 7 5 7 2 2 3 2),属于1至8随机出现的数字,
那么1前面的数字,1到8中,1出现1次,5出现2次,3出现1次,8出现1次,其余各数出现0次,
2前面的数字,1到8中,7出现1次,2出现1次,3出现1次,其余各数出现0次,
以此类推。。。
希望返回的形式和table一样,但是我希望得到的是一个矩阵,拜托各位大神了!!!


结果的格式如:
       1   2   3   4   5   6   7   8
1     1   0   1   0   2   0   0   1
2     0   1   1   0   0   0   1   0
3
4
5
6
7
8

376
myohmy88 发表于 2016-11-30 10:45:12
楼主你好~有个问题想请教以下。现有一个参数的分布,取值为0`1,但是不是连续的,函数形式是对于z in (1:500),要对dbinom(c*z,z,0.1)循环求和,c是所求参数。请问如何得到该分布的95%临界值?谢谢!

377
myohmy88 发表于 2016-11-30 10:45:22
楼主你好~有个问题想请教以下。现有一个参数的分布,取值为0`1,但是不是连续的,函数形式是对于z in (1:500),要对dbinom(c*z,z,0.1)循环求和,c是所求参数。请问如何得到该分布的95%临界值?谢谢!

378
jamesgd 发表于 2016-11-30 10:48:29
果真是 大神 呀, 这个贴要保存好

379
cbw1243 发表于 2016-11-30 12:19:18
myohmy88 发表于 2016-11-30 10:45
楼主你好~有个问题想请教以下。现有一个参数的分布,取值为0`1,但是不是连续的,函数形式是对于z in (1:50 ...
用最大似然估计进行参数估计,求出标准差,然后得到置信区间。

380
cbw1243 发表于 2016-11-30 12:22:33
3210sq 发表于 2016-11-30 00:11
小问题,小骚扰。

例如c(4 3 6 5 1 3 5 1 8 4 3 1 8 1 1 8 4 7 5 7 2 2 3 2),属于1至8随机出现的数字, ...
先构建n*8矩阵,如果n不大用for loop;如果很大就用 sapply。
至于随机的话用sample函数

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-31 15:46