我試圖運行一個GMM估計,其中包括正常CDF的積分,其參數同時存在於函數和積分區間中。代碼爲的要點如下: g1 <- function(b,x){
e <- b[1] +b[2]*x$x1
r <- e + b[3]*x$x2
n <- b[4]+b[5]*x$x3
ncdf <- function(z){
return((pnorm((log(z
我使用MATLAB來計算包括自然指數在內的複數函數的數值積分。 我得到一個警告: 無限或沒有非數字值遇到 ,如果我使用的功能integral,而另一個則會引發錯誤: 輸出的功能必須與輸入尺寸相同 如果我使用功能quadgk。 我認爲其原因可能是當變量ep接近零時,被積函數是無限的。 代碼如下所示。希望你們能幫我弄明白。 close all
clear
clc
%%
N = 10^5;
e
我想計算下面的積分,即積分內有積分。 integrate(function(v){
integrate(function(s){
pnorm(qnorm(v)-qnorm(s))
},0,1)$value
},0,1)
不幸的是我所有的時間得到兩種: Error in integrate(function(v) { :
evaluation of fu