individual_patents_t81 <- read_dta("D:/Rdata/data/individual_patents_t81.dta")
table8<-list(
result1 <- lm(patented ~ treated_post+treated+post,
data = individual_patents_t81[individual_patents_t81$appyr>=1995,]),
result2 <- lm(patented ~ treated_post+treated+appyr96+ appyr97+ appyr98+ appyr99+ appyr00+appyr01+ appyr02+ appyr03+appyr04+ appyr05+ appyr06+ appyr07+appyr08+ appyr09+ appyr10,
data = individual_patents_t81[individual_patents_t81$appyr>=1995,]),
result3 <- plm(patented ~ treated_post+treated+appyr96+ appyr97+ appyr98+ appyr99+ appyr00+appyr01+ appyr02+ appyr03+appyr04+ appyr05+ appyr06+ appyr07+appyr08+ appyr09+ appyr10|id_npo,
data = individual_patents_t81[individual_patents_t81$appyr>=1995,],
effect = 'individual',
model = 'within'),
)
数据本身的不大,但是观测值有10多万,没有加“|id_npo”可以跑出来,加了跑不出来。


雷达卡




京公网安备 11010802022788号







