我正在解決重建(或恢復)概率分佈函數時只知道分佈的矩的問題。我已經爲它編寫了代碼,雖然邏輯看起來是正確的,但我沒有得到我想要的輸出。在R中嵌套for循環 - 問題與最終結果
我試圖用作近似(或重建或恢復)CDF的公式是您在下圖中看到的。我正在編寫等式右邊的代碼,並將其與我在代碼中稱爲F的向量等同起來。
包含原始公式的紙張鏈接可以在此處找到。
它被標記爲(2)在紙方程。
這裏是我寫的:
#R Codes:
alpha <- 50
T <- 1
x <- seq(0, T, by = 0.1)
# Original CDF equation
Ft <- (1-log(x^3))*(x^3)
plot(x, Ft, type = "l", ylab = "", xlab = "")
# Approximated CDF equation using Moment type reconstruction
k<- floor(alpha*y/T)
for(i in 1:length(k))
{
for(j in k[i]:alpha)
{
F[x+1] <- (factorial(alpha)/(factorial(alpha-j)*factorial(j-k)*factorial(k)))*(((-1)^(j-k))/(T^j))*((9/(j+3))^2)
}
}
plot(x[1:7], F, type = "l", ylab = "", xlab = "")
任何幫助將這裏讚賞,因爲使用我的代碼獲得的逼近和圖形是從原來的曲線非常不同的代碼。
僅供參考,我刪除了您的內聯代碼減價塊的報價替代。這可以通過用四個空格縮進代碼塊來實現,或者在突出顯示要縮進的代碼後觸擊'{}'按鈕。你可以仔細檢查一下,我沒有破壞代碼的初始意圖嗎? – Chase
這幾乎不可能幫助你,因爲你的鏈接不公開。因此,我既沒有你的等式,也沒有你想要的結果。 – Andrie
當我運行你的代碼時,我得到一個錯誤,找不到y。我猜y是一個向量,因爲否則不是k = 1的長度? –