我正在运行一个简单的比赛时间与温度的回归,只是为了发展一些基本的直觉。我的数据集非常大,每个观察值都是给定比赛中一个单位在给定年份中的比赛完成时间。
对于起跑者,我正在运行一个非常简单的关于温度箱的比赛时间回归。
temp变量摘要:
|
Variable | Obs Mean Std. Dev Min Max
------------+--------------------------------------------
avg_temp_scc| 8309434 54.3 9.4
我看过很多相关的博主,他们解释了如何产生边际效应,但不知怎么的,他们中的大多数人并没有以交互效应为样本。
当我用Stata来做它的时候,边缘效应很容易计算,但是R似乎以另一种方式出现,并且显示了不同的结果。
我想知道如何用适当的相互作用计算适当的边际效应,
为什么函数margins()不显示交互因子?
如果我用手工coef(glm)*mean(dlogis(predict(glm, type = "link")),na.rm=T)计算它,它看起来更接近Stata结果,而不是使用函数margins()。
更重要的是,我如何手工计算standard error of m
我在Stata 15的拟合logistic回归模型中检索边缘效应时遇到了问题。结果变量mathtsbv是二进制变量,性别变量sex也是假的,记录的种族eth变量是绝对的,数值范围在0到5之间。所有缺失的值都被排除了。
下面是我的do文件的摘录:
logit mathtsbv sex eth sex##i.eth if (mathtsbv>=0&mathtsbv<.)&(sex>=0&sex<.)&(eth>=0ð<.)
margins, dydx(sex eth sex##i.eth) atmeans
这是我在S
我试图在Stata中编写一个foreach循环,它将自动生成数据集中所有连续变量(不包括字符串、二进制变量)的日志转换。我尝试过的代码(似乎不起作用)如下:
qui foreach v of varlist _all {
gen 'v' = log(v)
}
此外,我还想编写一个循环,它将生成所有可能的交互项排列,但不确定如何处理。
我使用xtivreg,xtabond实现了动态面板数据模型。在样本中,我有305组:
.xtsum id
+--------------+
| Observations |
+--------------+
| N=13588 |
| n=305 |
| T-bar=44.5508|
+--------------+
Stata只使用我示例中的292组作为xtabond。
Number of obs=11636
Number of groups=292
对于xtivreg,它只使用213组
Number of obs=8199
Number of groups =21
在Stata中,我最近发现,当我在一个回归模型中的多个交互项中使用相同的变量时,Stata会将该变量标记为共线性。例如,运行: regress dep i.gender##c.age i.ethnicity##c.age 标记以下消息: note: age omitted because of collinearity 年龄仍然包括在随后的回归表中,但有两次,第一次带有系数,SE等,正如人们所预期的,但第二次说的是omitted。 我以前做过很多次类似的分析,但从来没有(或者至少从来没有注意到)。这是不和谐的,因为不言而喻,年龄与其自身是共线的。但这并不重要,因为我试图同时输入的并不是两个名
我想要将一个分隔的文本文件导入Stata。有些字段是数字,其中数字用逗号(即2,144.20)格式化。当我在infix命令中为这些列指定数值数据类型时,这些值将被归为缺失。
infix 2 first str id 2-15 double amount 16-25 using "{datasetname}"
是否有一种方法来指定数字格式(例如%20.2fc),以便Stata不将它们视为非数字格式?另一种方法是将其作为字符串导入,并在以后将其转换为数字。但是,我想知道是否有一种方法可以在infix命令本身中指定格式。
根据此的答案,如何将共线列名打印到控制台?
脚本:
library(corrplot)
library(caret)
x <- seq(0, 100, 1)
# colinear with x
y <- x + 2.3
# almost colinear with x / some small gaussian noise
z <- x + rnorm(mean = 0, sd = 5, n = 101)
# uncorrrelated gaussian
w <- rnorm(mean = 0, sd = 1, n = 101)
a <- z+seq(10
我查了几个帖子,但没有得到我想要的答案。
我想运行一个Stata脚本作为Makefile的一部分。所以在我的目标中,我定义了
live_data:
@echo "Executing Stata code"
StataIC -e 'stata_code.do'
@exit_status=$?
@echo "Finished execution of Stata code."
@echo Code finished with exit code of $(exit_status)
但是,每当我查看make al