Saya mencoba menganalisis lead-lag antara time series dua harga saham. Dalam analisis deret waktu reguler, kita dapat melakukan Cross Correlaton, VECM (Granger Causality). Namun bagaimana seseorang menangani hal yang sama dalam rangkaian waktu yang tidak beraturan. Hipotesisnya adalah bahwa salah...
Saya ingin memodelkan dua variabel waktu yang berbeda, beberapa di antaranya sangat collinear dalam data saya (usia + kohort = periode). Melakukan ini saya mengalami beberapa masalah dengan lmerdan dan interaksi poly(), tetapi mungkin tidak terbatas lmer, saya mendapat hasil yang sama dengan...
Gambar di bawah ini dari artikel ini di Ilmu Psikologi . Seorang kolega menunjukkan dua hal yang tidak biasa tentang hal itu: Menurut keterangan, bilah kesalahan menunjukkan "± 2,04 kesalahan standar, interval kepercayaan 95%." Saya hanya pernah melihat ± 1,96 SE digunakan untuk CI 95%, dan saya...
Saya melakukan ANOVA satu arah (per spesies) dengan kontras khusus. [,1] [,2] [,3] [,4] 0.5 -1 0 0 0 5 1 -1 0 0 12.5 0 1 -1 0 25 0 0 1 -1 50 0 0 0 1 di mana saya membandingkan intensitas 0,5 melawan 5, 5 melawan 12,5 dan seterusnya. Ini adalah data yang saya kerjakan dengan hasil sebagai...
Saya perhatikan bahwa ketika saya mencoba untuk menemukan nilai kritis untuk Mann-Whitney U menggunakan R, nilainya selalu 1 + nilai kritis. Misalnya, untuk , nilai kritis (berekor dua) adalah 8, sedangkan untuk , the (berekor dua) ) nilai kritis adalah 22 (periksa tabel ),
Dalam contoh berikut ini saya memiliki kerangka data yang terdiri dari serangkaian waktu pengukuran suhu air yang direkam pada 5 kedalaman di lautan di mana setiap nilai Tempsesuai dengan tanggal DateTimedan kedalaman Depth. set.seed(1) Temp <- rnorm(43800,sd=20) AirT <-
Model koin bias biasanya memiliki satu parameter . Salah satu cara untuk memperkirakan dari serangkaian undian adalah dengan menggunakan beta sebelum dan menghitung distribusi posterior dengan kemungkinan binomial.θ=P(Head|θ)θ=P(Head|θ)\theta = P(\text{Head} | \theta)θθ\theta Dalam pengaturan...
Saya telah memasang model lognormal menggunakan R dengan satu set data. Parameter yang dihasilkan adalah: meanlog = 4.2991610 sdlog = 0.5511349 Saya ingin mentransfer model ini ke Scipy, yang belum pernah saya gunakan sebelumnya. Menggunakan Scipy, saya bisa mendapatkan bentuk dan skala 1 dan...
Struktur data > str(data) 'data.frame': 6138 obs. of 10 variables: $ RT : int 484 391 422 516 563 531 406 500 516 578 ... $ ASCORE : num 5.1 4 3.8 2.6 2.7 6.5 4.9 2.9 2.6 7.2 ... $ HSCORE : num 6 2.1 7.9 1 6.9 8.9 8.2 3.6 1.7 8.6 ... $ MVMNT : Factor w/ 2 levels "_Withd","Appr": 2 2 1 1 2...
Saya tidak mendapatkan perbedaan antara rfobject$importancedan importance(rfobject)di kolom MeanDecreaseAccuracy. Contoh: > data("iris") > fit <- randomForest(Species~., data=iris, importance=TRUE) > fit$importance setosa versicolor virginica MeanDecreaseAccuracy
Pertanyaan saya berkaitan dengan GAM dalam paket mgcv R. Karena ukuran sampel yang kecil saya ingin menentukan kesalahan prediksi menggunakan validasi silang leave-one-out. Apakah ini masuk akal? Apakah ada paket atau kode bagaimana saya bisa melakukan ini? The errorest()fungsi dalam ipred paket...
Saya memiliki percobaan yang dijalankan pada ratusan komputer yang didistribusikan di seluruh dunia yang mengukur kejadian peristiwa tertentu. Peristiwa masing-masing bergantung satu sama lain sehingga saya dapat memesannya dalam urutan yang meningkat dan kemudian menghitung perbedaan...
Saya punya pertanyaan tentang perbandingan model menggunakan faktor Bayes. Dalam banyak kasus, ahli statistik tertarik untuk menggunakan pendekatan Bayesian dengan prior yang tidak tepat (misalnya beberapa prior Jeffreys dan prior reference). Pertanyaan saya adalah, dalam kasus-kasus di mana...
Saya memiliki GLMM formulir: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Ketika saya menggunakan drop1(model, test="Chi"), saya mendapatkan hasil yang berbeda daripada jika saya menggunakan Anova(model, type="III")dari paket mobil atau...
Saya ingin melakukan analisis kekuatan untuk sampel tunggal dari data binomial, dengan H0:p=0H0:p=0H_0: p = 0 , vs H1:p=0.001H1:p=0.001H_1: p = 0.001 , di mana ppp adalah proporsi keberhasilan dalam populasi. Jika 0<p<10<p<10 < p
Apa yang terjadi ketika Anda tidak memiliki gagasan tentang distribusi parameter? Pendekatan apa yang harus kita gunakan? Sebagian besar waktu kami bertujuan untuk menggarisbawahi jika variabel tertentu memiliki pengaruh terhadap ada / tidaknya spesies tertentu, dan variabel tersebut diterima atau...
Saya merasa kewalahan setelah mencoba menggali literatur tentang bagaimana menjalankan analisis model campuran saya dengan menggunakan AIC untuk memilih model atau model terbaik. Saya tidak berpikir bahwa data saya serumit itu, tetapi saya mencari konfirmasi bahwa apa yang telah saya lakukan adalah...
Saya mencoba menyesuaikan model campuran multivarian (yaitu, beberapa respons) di R. Selain dari paket ASReml-rdan SabreR(yang memerlukan perangkat lunak eksternal), sepertinya ini hanya mungkin di MCMCglmm. Dalam makalah yang menyertai MCMCglmmpaket (hal.6), Jarrod Hadfield menggambarkan proses...
Utama masalah adalah: Saya tidak dapat memperoleh estimasi parameter yang sama dengan EViews dan R. Untuk alasan yang saya sendiri tidak tahu, saya perlu memperkirakan parameter untuk data tertentu menggunakan EViews. Ini dilakukan dengan memilih opsi NLS (nonlinear least square) dan menggunakan...
Salah satu metode yang disarankan kepada saya adalah dengan melihat plot scree dan memeriksa "siku" untuk menentukan jumlah PC yang tepat untuk digunakan. Tetapi jika plotnya tidak jelas, apakah R memiliki perhitungan untuk menentukan angka? fit <- princomp(mydata,