这条命令如何解读

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

d<-if (rdf>0) 2 *(1-pt(a,rdf))
为什么第一个括号后面就直接一个2，而没有任何连接符？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：如何

假设
rdf=2
a=0.5
很显然rdf>0  #TRUE
所以会执行 2 *(1-pt(a,rdf))
d<-if (rdf>0) 2 *(1-pt(a,rdf))
#0.6666667

反之如果设
rdf=-2
rdf>0  #FALSE
就不会执行 2 *(1-pt(a,rdf))

谢谢您，R语法真是泰简洁了，一个条件语句直接这样写了
一般都是有加括号或者换行的，呵呵，没有看过来

R语句有时候省略了赋值，d<-if (rdf>0) 2 *(1-pt(a,rdf))可以写成d<-if (rdf>0)  d<-2 *(1-pt(a,rdf))

ywh19860616 发表于 2011-8-5 11:10
谢谢您，R语法真是泰简洁了，一个条件语句直接这样写了
一般都是有加括号或者换行的，呵呵，没有看过来

呵呵，是的。不过这不算R的简介，很多语言都可以这样，包括比较基础的C++

epoh老师，您好
1、我看了您提供的systemfitFiml程序，一直运行不出来的，有一点很不解
在systemfit包中，我看了这个函数内部程序，其中也有if (method == "FIML") ，但是为什么命令Description
时只有提到这些方法：
Fits a set of linear structural equations using Ordinary Least Squares (OLS), Weighted Least Squares (WLS), Seemingly Unrelated Regression (SUR), Two-Stage Least Squares (2SLS), Weighted Two-Stage Least Squares (W2SLS) or Three-Stage Least Squares (3SLS).
我试图用method=“FIML”，但是一直报错，不能运行
2、在systemfit包的Vignettes介绍里面4.4 panel data，提到这样一句话：
The systemfit function can also be used for a (classical)"seemingly unrelated regression" analysis with panel data.
这里提到了，针对panel data只能做SUR估计，而其他估计方法没有针对panel data。且此时只能为单方程的。提到了这样一句话：
In contrast to the previously described usage of systemfit,argument formula must be a single equation(object of class formula).
不知道我理解的是否正确？

哈哈!老兄都被你说对了

page 32/40 也的确指出

full information maximum likelihood (FIML)还未完成

而且仅配合＂SUR＂

今天尝试修改,可惜交了白券.

############
library("systemfit")
data("GrunfeldGreene")
library("plm")
GGPanel <- plm.data(GrunfeldGreene, c("firm", "year"))
formulaGrunfeld <- invest ~ value + capital
fitsur <- systemfit(formulaGrunfeld, method = "SUR",data = GGPanel)
fitsur
Coefficients:
        Chrysler_(Intercept)               Chrysler_value
                   0.5043036                    0.0695456
            Chrysler_capital General.Electric_(Intercept)
                   0.3085445                  -22.4389132
      General.Electric_value     General.Electric_capital
                   0.0372914                    0.1307830
  General.Motors_(Intercept)         General.Motors_value
                -162.3641052                    0.1204930
      General.Motors_capital         US.Steel_(Intercept)
                   0.3827462                   85.4232548
              US.Steel_value             US.Steel_capital
                   0.1014782                    0.3999914
    Westinghouse_(Intercept)           Westinghouse_value
                   1.0888770                    0.0570091
        Westinghouse_capital
                   0.0415065
###############
greeneSurPooled <- systemfit(formulaGrunfeld, method = "SUR",data = GGPanel, pooled = TRUE)
greeneSurPooled
Coefficients:
        Chrysler_(Intercept)               Chrysler_value
                 -38.3612772                    0.0961894
            Chrysler_capital General.Electric_(Intercept)
                   0.3095321                  -38.3612772
      General.Electric_value     General.Electric_capital
                   0.0961894                    0.3095321
  General.Motors_(Intercept)         General.Motors_value
                 -38.3612772                    0.0961894
      General.Motors_capital         US.Steel_(Intercept)
                   0.3095321                  -38.3612772
              US.Steel_value             US.Steel_capital
                   0.0961894                    0.3095321
    Westinghouse_(Intercept)           Westinghouse_value
                 -38.3612772                    0.0961894
        Westinghouse_capital
                   0.3095321

##################

Estimation methods
In the future, we wish to include more sophisticated estimation methods
such as limited information maximum likelihood (LIML),
full information maximum likelihood (FIML), generalized
methods of moments (GMM) and spatial econometric methods
(Paelinck and Klaassen 1979; Anselin 1988).

##################

Unit_root

   Unit_root.pdf (202.01 KB)

2011-8-25 19:19:17 上传

Econometric Methods for Panel Data.pdf

Econometric Methods for Panel Data.pdf (105.2 KB)

2011-9-27 15:56:29 上传

恩，现在对simultaneous equation model for panel data的估计还不懂，目前R的systemfit包只有截面的
如果分别对单方程进行IV估计，那样不是有效的，所以看到别人用了FIML方法，才试图找FIML方法怎么实现

谢谢epoh老师，又花费您宝贵时间了
可能这个方法目前还不成熟。