下面的手绘图显示了对泊松回归拟合值绘制的残差。
如何创建两个向量,在绘制时,显示与附图中相同的模式?
请注意,拟合值无关紧要(因此没有比例)。
答案 0 :(得分:5)
使y的标准偏差成为x的增加函数:
x <- runif(100, 0, 10)
y <- rnorm(100, 0, x)
plot(x,y)
答案 1 :(得分:2)
这是我的尝试(将Stata代码从http://www.econometricsbysimulation.com/2012/11/modeling-heteroskedasticity.html翻译为R):
set.seed(123)
# Let's generate a sample data set that has heteroskedasticity in z.
z <- runif(1000)*5
# This means the standard deviation on u is 18 plus 16*z (which has a mean of 2.5)
u <- ((18+16*z))*rnorm(1000)
# Plot
plot(z, u)
答案 2 :(得分:0)