Parameter apa start, etastart, mustartdi GLM () function ? Saya telah mencari di dokumen dan internet tetapi saya belum menemukan penjelasan yang jelas apa artinya ini. Ini menyerupai "nilai awal" bayesian untuk rantai, tapi saya ragu ini terkait, karena fungsi glm () dalam R adalah statistik...
Bagaimana matriks kesalahan var / cov dihitung dengan paket analisis statistik dalam praktik? Ide ini jelas bagi saya dalam teori. Tetapi tidak dalam praktik. Maksud saya, jika saya memiliki vektor variabel acak X =( X1, X2, ... , Xn)⊤X=(X1,X2,...,Xn)⊤\textbf{X}=(X_{1}, X_{2}, \ldots, X_{n})^\top...
Saya melakukan regresi logistik multivariat dengan variabel dependen Yadalah kematian di panti jompo dalam periode masuk tertentu dan mendapatkan hasil berikut (perhatikan jika variabel yang dimulai di Adalamnya adalah nilai yang berkelanjutan sedangkan yang mulai dalam Bkategori): Call: glm(Y ~...
Saya memiliki kumpulan data deret waktu yang saya coba paskan dengan Hidden Markov Model (HMM) untuk memperkirakan jumlah status laten dalam data. Kode pseudo saya untuk melakukan ini adalah sebagai berikut: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ......
Saya menyesuaikan model efek acak dengan glmerbeberapa data bisnis. Tujuannya adalah untuk menganalisis kinerja penjualan oleh distributor, dengan mempertimbangkan variasi regional. Saya memiliki variabel berikut: distcode: ID distributor, dengan sekitar 800 level region: ID geografis tingkat...
Tutup. Pertanyaan ini di luar topik . Saat ini tidak menerima jawaban. Ingin meningkatkan pertanyaan ini? Perbarui pertanyaan sehingga sesuai topik untuk Cross Validated. Ditutup 6 tahun yang lalu . Saya telah menggunakan file teks untuk menyimpan data...
Apakah wikipedia salah ... atau saya tidak memahaminya? Wikipedia: Kotak putih dan hitam ("pola catur") tersebar dengan sempurna sehingga Moran's saya akan −1. Jika kotak putih ditumpuk ke setengah papan dan kotak hitam ke yang lain, Moran's saya akan mendekati +1. Susunan acak warna kotak akan...
Sunting utama: Saya ingin mengucapkan terima kasih yang sebesar-besarnya kepada Dave & Nick sejauh ini atas tanggapan mereka. Kabar baiknya adalah bahwa saya mendapatkan loop untuk bekerja (prinsip dipinjam dari posting Prof Hydnman pada peramalan batch). Untuk mengkonsolidasikan kueri yang...
Saya punya beberapa kode dan keluaran, dan saya ingin membuat model. Saya tidak tahu bagaimana membangun model menggunakan output ini: require("splines") x <- c(0.2, 0.23, 0.26, 0.29, 0.33, 0.46, 0.53 ) y <- c(0.211, 0.2026, 0.2034, 0.2167, 0.2177, 0.19225, 0.182) fit <- lm(y ~...
Saya mendapat kesan bahwa fungsi lmer()dalam lme4paket tidak menghasilkan nilai-p (lihat lmer, nilai-p dan semua itu ). Saya telah menggunakan MCMC menghasilkan nilai p sebagai per pertanyaan ini: Efek signifikan dalam lme4model campuran dan pertanyaan ini: Tidak dapat menemukan nilai-p dalam...
Saya bukan ahli matematika. Saya telah mencari di internet tentang KL Divergence. Apa yang saya pelajari adalah divergensi KL mengukur informasi yang hilang ketika kami memperkirakan distribusi model sehubungan dengan distribusi input. Saya telah melihat ini di antara dua distribusi kontinu atau...
Saya kesulitan menemukan cara untuk menghitung nilai-p untuk area di bawah karakteristik operator penerima (ROC). Saya memiliki variabel kontinu dan hasil tes diagnostik. Saya ingin melihat apakah AUROC signifikan secara statistik. Saya menemukan banyak paket yang berhubungan dengan kurva ROC:...
Dari hasil saya, tampak bahwa GLM Gamma memenuhi sebagian besar asumsi, tetapi apakah ini merupakan peningkatan yang berharga atas LM yang ditransformasikan log? Kebanyakan literatur yang saya temukan berhubungan dengan Poisson atau Binomial GLMs. Saya menemukan artikel EVALUASI ASUMSI MODEL LINEAR...
Saya punya dataset yang mewakili 1000 dokumen dan semua kata yang muncul di dalamnya. Jadi baris mewakili dokumen dan kolom mewakili kata-kata. Jadi misalnya, nilai dalam sel berarti waktu kata j muncul dalam dokumen i . Sekarang, saya harus menemukan 'bobot' dari kata-kata itu, menggunakan metode...
Saya menggunakan paket caret untuk melatih objek randomForest dengan 10x10CV. library(caret) tc <- trainControl("repeatedcv", number=10, repeats=10, classProbs=TRUE, savePred=T) RFFit <- train(Defect ~., data=trainingSet, method="rf", trControl=tc, preProc=c("center", "scale")) Setelah...
Pertanyaan Skor tes dari tiga kelompok orang disimpan sebagai vektor terpisah dalam R. set.seed(1) group1 <- rnorm(100, mean = 75, sd = 10) group2 <- rnorm(100, mean = 85, sd = 10) group3 <- rnorm(100, mean = 95, sd = 10) Saya ingin tahu apakah ada perbedaan signifikan dalam median...
Untuk model linier Gaussian mana diasumsikan terletak pada ruang vektor dan memiliki distribusi normal standar pada , statistik untuk mana adalah ruang vektor, adalah fungsi satu-ke-satu yang meningkat dari statistik penyimpangan : Bagaimana kita bisa tahu bahwa statistik ini menyediakan tes...
Saya ingin memahami uji pasti fisher dengan lebih baik, jadi saya menyusun contoh mainan berikut ini, di mana f dan m berhubungan dengan pria dan wanita, dan n dan y berhubungan dengan "konsumsi soda" seperti ini: > soda_gender f m n 0 5 y 5 0 Jelas, ini adalah penyederhanaan yang...
Saya mencoba menghitung log-kemungkinan untuk regresi kuadrat terkecil nonlinear umum untuk fungsi dioptimalkan olehfungsi dalam paket R, menggunakan matriks varians kovarians yang dihasilkan oleh jarak pada pohon filogenetik dengan asumsi gerakan Brown (daripaket). Kode R yang dapat direproduksi...
Saya sedang memeriksa bagian dari dataset saya yang berisi 46840 nilai ganda mulai dari 1 hingga 1690 yang dikelompokkan dalam dua grup. Untuk menganalisis perbedaan antara kelompok-kelompok ini saya mulai dengan memeriksa distribusi nilai-nilai untuk memilih tes yang tepat. Mengikuti panduan...