normal-distribution

0熱度

1回答

你好，我希望得到一些幫助，以整合基於隨機變量的函數。該功能是獲得連續分配的預期值。這是我到目前爲止的代碼。 montecarlo = function(r,v,t,x,k) { y = rnorm(1) e = (y*(x*exp((-v*sqrt(t)*y)+((r-(.5*v^2))*t))-k)) MCOP=exp(-r*t)*integrate(e, lower = -Inf, u

3熱度

1回答

確定正態分佈給予其位數信息

我不知道我怎麼會具有R告訴我SD（如qnorm一個參數（）內建R中）正態分佈，其95％的極限值已知？作爲一個例子，我知道我的法線的兩個95％極限值分別是158和168。因此，在下面的R代碼SD顯示爲「x」。如果「Y」（在這個簡單的qnorm（）功能的回答）需要爲（158，168），然後可以告訴[R我應該是什麼X？ y <- qnorm(c(.025,.975), 163, x)

2熱度

1回答

dnorm和dmvnorm結果不匹配

我注意到從dmvnorm（mvtnorm包）返回到dnorm（baseR）的值之間的實質區別。哪一個是正確的結果，爲什麼？在獨立正態隨機變量的情況下，我應該使用哪一個？ > mvtnorm::dmvnorm(rep(1,2),rep(0,2), diag(2,2)) [1] 0.04826618 > prod(dnorm(rep(1,2),0,2^2)) [1] 0.009344515

2熱度

2回答

ARIMA模型與非正常的錯誤

我適合ARIMA(0,0,1)模型在R與一個外生變量。性擬合後，我測試的誤差項，它的高度非正常（這就像t-distributed錯誤）：我的問題是：是否有任何包R可以容納ARIMA模型t-distributed錯誤？或者有沒有其他補救措施來解決這個問題？數據已經過對數轉換的數據，所以我想我不能執行另一個數據轉換。非常感謝您的幫助！下面是數據： dput(x) c(1.098612289

1熱度

1回答

在python中用scipy截斷正態分佈

我想在Python3中用scipy來截斷正態分佈。我希望做一些簡單的：情節截斷正以0.5爲中心，範圍從0到1，我有以下的代碼行 from scipy import truncnorm import matplotlib.pyplot as plt plt.plot([truncnorm.pdf(p,0,1, loc=0.5) for p in np.arange(0,1.1,0.1)])

2熱度

1回答

用拒絕方法創建正態分佈產生錯誤的前因子

我試圖用拒絕方法在Fortran中的均勻分佈值中創建正態分佈。它實際上或多或少地運作良好，但我並不完全符合我想要的結果。我生成與該段的代碼 function generator result(c) implicit none integer, dimension(2) :: clock double precision :: c,d cal

1熱度

1回答

Python：計算N個多元正態分佈值的可能性

所以我有一組N個多元正態分佈，它們都具有相同的協方差。對於這些分佈中的每一個，我想計算得到值x的可能性。對於單個分佈情況，以及「X」值倍數，這是微不足道現在 from scipy.stats import multivariate_normal import numpy as np cov = [[1 ,0.1],[0.1 ,1]] mean = [0,0] Values = np.r

0熱度

2回答

如何將長python列表格式化爲Fortran多行數組？

我有一個問題。我在Fortran編程一個數學計算程序。在那裏，我必須初始化一個隨機值數組。這些值必須正態分佈，平均值爲0，標準偏差爲1. 我在Python 2.7中做了以下腳本來生成900個這樣的值。 import numpy as np mu, sigma = 0, 1.0 list = [] i = 0 while i < 901: s = np.random.norma

0熱度

1回答

R 1000 rnorm samples

x < - rorm（25）將從標準正態分佈產生一個大小爲25的單個樣本。我該如何從標準正態分佈中同時抽取1000個大小爲25的樣本？我希望能夠有效地做到這一點，以便我能夠計算1000個樣本中每個樣本的平均值和標準偏差等事情，並通過直方圖進行比較。 [另外：我將隨後要uniformally和隨機地選擇這些1000個樣品中的一個和它引導]

1熱度

1回答

爲夏皮羅Wilk檢驗臨界值

我試圖得到一個夏皮羅Wilk檢驗臨界W值在R. Shapiro-Wilk normality test data: samplematrix[, 1] W = 0.69661, p-value = 7.198e-09 其中n = 50和α= 0.05，我知道臨界值W = .947，通過執行臨界值表。但是，如何使用R來獲得這個臨界值？