Diberikan matriks kovarians , bagaimana cara menghasilkan data sedemikian rupa sehingga memiliki sampel matriks kovarians \ hat {\ boldsymbol \ Sigma} = \ boldsymbol \ Sigma_s ?Σ = Σ sΣsΣs\boldsymbol \Sigma_sΣ^=ΣsΣ^=Σs\hat{\boldsymbol \Sigma} = \boldsymbol \Sigma_s Lebih umum: kita sering...
Bagaimana cara kerja metode inversi? Katakanlah saya memiliki sampel acak X1, X2, . . . , XnX1,X2,...,XnX_1,X_2,...,X_n dengan kepadatan lebih dari dan karenanya dengan cdf pada . Kemudian dengan metode inversi saya mendapatkan distribusi sebagai . 0<x<1FX(x)=x1/θ(0,1)XF - 1 X(u)=uθf( x ; θ...
Misalnya, dalam R, MASS::mvrnorm()fungsi ini berguna untuk menghasilkan data untuk menunjukkan berbagai hal dalam statistik. Dibutuhkan Sigmaargumen wajib yang merupakan matriks simetris yang menentukan matriks kovarians dari variabel. Bagaimana saya membuat simetris n×nn×nn\times n matriks dengan...
Saya bertanya-tanya apakah mungkin untuk menghasilkan variabel binomial acak berkorelasi mengikuti pendekatan transformasi linier? Di bawah ini, saya mencoba sesuatu yang sederhana dalam R dan menghasilkan beberapa korelasi. Tapi saya bertanya-tanya apakah ada cara berprinsip untuk melakukan...
Saya mencoba beberapa pencarian google biasa dll tetapi sebagian besar jawaban yang saya temukan agak ambigu atau bahasa / perpustakaan tertentu seperti Python atau C ++ stdlib.hdll. Saya mencari bahasa agnostik, jawaban matematis, bukan spesifik perpustakaan. Sebagai contoh, banyak yang...
Bagaimana saya bisa mengambil sampel dari distribusi campuran, dan khususnya campuran distribusi normal R? Misalnya, jika saya ingin mengambil sampel dari: 0.3×N(0,1)+0.5×N(10,1)+0.2×N( 3 , .1 )0,3×N(0,1)+0,5×N(10,1)+0,2×N(3,.1) 0.3\!\times\mathcal{N}(0,1)\; + \;0.5\!\times\mathcal{N}(10,1)\; +...
Saya membaca tentang MCMC adaptif (lihat misalnya, Bab 4 dari Handbook of Markov Chain Monte Carlo , ed. Brooks et al., 2011; dan juga Andrieu & Thoms, 2008 ). Hasil utama dari Roberts dan Rosenthal (2007) adalah bahwa jika skema adaptasi memenuhi kondisi adaptasi menghilang (ditambah beberapa...
Saya mencoba untuk menulis sebuah program dalam R yang mensimulasikan angka acak semu dari distribusi dengan fungsi distribusi kumulatif: F(x)=1−exp(−ax−bp+1xp+1),x≥0F(x)=1−exp(−ax−bp+1xp+1),x≥0F(x)= 1-\exp \left(-ax-\frac{b}{p+1}x^{p+1}\right), \quad x \geq 0 di...
Apa cara terbaik untuk mengambil sampel dari distribusi Cantor ? Ini hanya memiliki cdf dan kami tidak dapat
Saya menggunakan dekomposisi Cholesky untuk mensimulasikan variabel acak berkorelasi diberi matriks korelasi. Masalahnya, hasilnya tidak pernah mereproduksi struktur korelasi seperti yang diberikan. Berikut ini adalah contoh kecil dalam Python untuk menggambarkan situasi. import numpy as np...
Saya mencoba untuk menghasilkan urutan acak berkorelasi dengan mean = 000 , varians = 111 , koefisien korelasi = 0.80.80.8 . Dalam kode di bawah ini, saya menggunakan s1& s2sebagai standar deviasi, dan m1& m2sebagai sarana. p = 0.8 u = randn(1, n) v = randn(1, n) x = s1 * u + m1 y = s2 *...
Di sini 'pemotongan' menyiratkan pengurangan presisi angka acak dan tidak memotong rangkaian angka acak. Sebagai contoh, jika saya memiliki angka yang benar-benar acak (diambil dari distribusi apa pun, misalnya, normal, seragam, dll.) Dengan presisi sewenang-wenang dan saya memotong semua angka...
Saya ingin menggambar bilangan bulat dari 1 ke beberapa tertentu dengan menggulirkan sejumlah dadu enam sisi yang adil (d6). Jawaban yang baik akan menjelaskan mengapa metodenya menghasilkan bilangan bulat seragam dan independen .NNN Sebagai contoh ilustrasi, akan sangat membantu untuk menjelaskan...
Saya telah menulis sebuah program yang menghasilkan data acak. Jika program bekerja dengan benar, data tersebut harus mengikuti distribusi probabilitas khusus yang diketahui. Saya ingin menjalankan program, melakukan perhitungan pada hasilnya, dan menghasilkan nilai-p. Sebelum orang lain...
Saya ingin menghasilkan sampel dari wilayah biru yang didefinisikan di sini: Solusi naif adalah menggunakan sampel penolakan di unit square, tetapi ini hanya memberikan efisiensi (~ 21,4%).1−π/41−π/41-\pi/4 Apakah ada cara agar saya dapat mencicipi lebih efisien?
Saya perlu membuat angka acak mengikuti distribusi normal dalam interval . (Saya bekerja di R.)(a,b)(a,b)(a,b) Saya tahu fungsi ini rnorm(n,mean,sd)akan menghasilkan angka acak mengikuti distribusi normal, tetapi bagaimana cara menetapkan batas interval di dalamnya? Apakah ada fungsi R tertentu...
Saya ingin menggambar sampel . Wikipedia menyarankan untuk menggunakan komposisi Cholesky atau Eigend , yaitu atau x∼N(0,Σ)x∼N(0,Σ)\mathbf{x} \sim N\left(\mathbf{0}, \mathbf{\Sigma} \right)Σ=D1DT1Σ=D1D1T \mathbf{\Sigma} = \mathbf{D}_1\mathbf{D}_1^T Σ = QΛQTΣ=QΛQT \mathbf{\Sigma} =...
Dari randoness statistik Wikipedia : Keacakan global dan keacakan lokal berbeda. Kebanyakan konsepsi filosofis tentang keacakan bersifat global — karena mereka didasarkan pada gagasan bahwa "dalam jangka panjang" suatu urutan tampak benar-benar acak, bahkan jika sub-urutan tertentu tidak akan...
Misalkan kita punya X 2 ∼ unif ( n , 0 , 1 ) ,X1∼unif(n,0,1),X1∼unif(n,0,1),X_1 \sim \textrm{unif}(n,0,1), X2∼unif(n,0,1),X2∼unif(n,0,1),X_2 \sim \textrm{unif}(n,0,1), di mana adalah sampel acak ukuran seragam,
Saya memperkirakan sampel matriks kovarian CCC sampel dan mendapatkan matriks simetris. Dengan CCC , saya ingin membuat nnn -variate rn didistribusikan normal, tetapi karena itu saya membutuhkan Cholesky dekomposisi . Apa yang harus saya lakukan jika tidak pasti