小编bri*_*ian的帖子

据我所知,Stata和R都有"预测"功能.我试图复制使用R在Stata中执行的结果,结果涉及计算预测值的标准偏差.R中是否有功能,可能使用其"预测"功能,这将允许我这样做？我似乎无法完美地复制结果.如果它有帮助,Stata代码执行以下操作:

reg Y X1 X2 if condition
predict resid, r
predict stdf, stdf

stdf论证的定义是:

stdf计算预测的标准误差,这是1次观测的点预测的标准误差.它通常被称为未来或预测值的标准误差.通过构造,产生的标准误差stdf 总是大于由此产生的误差stdp; 参见[R]预测中的方法和公式

我写的R代码是:

fit <- lm(Y ~ X1 + X2, data=df)
new.df <- data.frame(...) # This is a new data frame with my new data period I want to predict in
predict(fit, new.df, se.fit = TRUE)

但是,当我将标准误差转换为标准偏差时,它们与Stata输出不匹配.

提前致谢!

我正在写python中的csv,我想知道为什么我的输出文件在程序完成运行时是空白的.首先,这是我的代码:

reader = csv.reader(open(location, 'rU'), delimiter = ',', quotechar = '"')
writer = csv.writer(open('temp.csv', 'wb'), delimiter = ',', quotechar = '"')

writer.writerow(["test", "test", "test", "test", "test", "test", "test"])

count = 0
for row in reader:
    if count != 0 and count < len(split):
        writer.writerow([row[0], row[1], row[2], row[3], row[4], row[5], split[count-1]])
    count = count + 1

我知道一些事情:

• 当程序完成时,写入'temp.csv'的所有内容都是"test"行.文档的其余部分是空的.
• 如果我强制在运行过程中停止程序,实际上文本正在写入文件:它只是因某种原因被擦除.
• 我做了一些研究,很多人发现"关闭"csv文件解决了他们的问题,但这对我来说似乎没有用,除非我遗漏了一些明显的东西!

在此先感谢您的帮助!

编辑:这是我正在使用的一些更新代码:

infile = open(location, 'rU')
outfile = open('temp.csv', 'wb')
reader = csv.reader(infile, delimiter = ',', quotechar = '"')
writer = …