Pertanyaan yang diberi tag «r»

Inilah konteks saya untuk pertanyaan ini: Dari apa yang bisa saya katakan, kita tidak bisa menjalankan regresi kuadrat terkecil biasa dalam R ketika menggunakan data tertimbang dan surveypaket. Di sini, kita harus menggunakan svyglm(), yang sebaliknya menjalankan model linear umum (yang mungkin hal yang sama? Saya bingung di sini dalam …

13 r  generalized-linear-model  least-squares  r-squared  deviance 

Prediksi saya yang berasal dari model regresi logistik (glm dalam R) tidak dibatasi antara 0 dan 1 seperti yang saya harapkan. Pemahaman saya tentang regresi logistik adalah bahwa parameter input dan model Anda digabungkan secara linear dan responsnya diubah menjadi probabilitas menggunakan fungsi tautan logit. Karena fungsi logit dibatasi antara …

13 r  regression  logistic  generalized-linear-model 

Jika saya membuat matriks 2-D yang seluruhnya terdiri dari data acak, saya berharap komponen PCA dan SVD pada dasarnya tidak menjelaskan apa-apa. Sebaliknya, sepertinya kolom SVD pertama muncul untuk menjelaskan 75% dari data. Bagaimana ini bisa terjadi? Apa yang saya lakukan salah? Ini plotnya: Berikut adalah kode R: set.seed(1) rm(list=ls()) …

13 r  pca  svd 

The bildpaket tampaknya menjadi paket yang sangat baik untuk respon biner serial. Tapi ini untuk waktu diskrit. Saya ingin menentukan fungsi waktu yang mulus untuk koneksi rasio peluang dari respons Y saat ini dengan respons biner yang diukur pada waktu sebelumnya, atau setidaknya versi Markov orde pertama dari ini. Saya …

13 r  repeated-measures  binary-data  panel-data 

Saya mencari paket R umum, bersih dan cepat (yaitu menggunakan C ++ rutin) untuk mensimulasikan jalur dari difusi nonlinear non-homogen seperti (1) menggunakan skema Euler-Maruyama, skema Milstein (atau lainnya). Ini ditakdirkan untuk tertanam ke dalam kode estimasi yang lebih besar dan karenanya layak dioptimalkan. dXt= f( θ , t , …

13 r  simulation  stochastic-processes  markov-process 

Saya ingin memprediksi Tree Heights di area tertentu menggunakan beberapa variabel yang diperoleh melalui penginderaan jauh. Seperti perkiraan Biomassa, dll. Saya ingin menggunakan regresi linier terlebih dahulu (saya tahu itu bukan ide terbaik tetapi merupakan langkah wajib untuk proyek saya). Saya ingin tahu bagaimana autokorelasi spasial dapat memengaruhinya dan apa …

13 r  multiple-regression  spatial  autocorrelation 

"Machine Learning For Hackers" O'Reilly mengatakan bahwa setiap komponen utama mewakili persentase perbedaan. Saya mengutip bagian yang relevan dari halaman di bawah ini (bab 8, hal.207). Berbicara dengan ahli lain, mereka sepakat itu persentase. Namun 24 komponen berjumlah 133,2095%. Bagaimana itu bisa terjadi? Setelah meyakinkan diri sendiri bahwa kita dapat …

13 r  pca 

Saya mencoba menjalankan Bayesian logit pada data di sini . Saya menggunakan bayesglm()dalam armpaket di R. Pengodean cukup mudah: df = read.csv("http://dl.dropbox.com/u/1791181/bayesglm.csv", header=T) library(arm) model = bayesglm(PASS ~ SEX + HIGH, family=binomial(link="logit"), data=df) summary(model) memberikan hasil sebagai berikut: Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.10381 0.10240 1.014 0.311 …

13 r  bayesian  p-value 

Saya akan senang untuk melakukan tes post-hoc TukeyHSD setelah Anova dua arah saya dengan R, memperoleh tabel berisi pasangan yang disortir yang dikelompokkan berdasarkan perbedaan yang signifikan. (Maaf tentang kata-katanya, saya masih baru dengan statistik.) Saya ingin memiliki sesuatu seperti ini: Jadi, dikelompokkan dengan bintang atau huruf. Ada ide? Saya …

13 r  anova  multiple-comparisons  post-hoc  tukey-hsd 

Dua fungsi ini ada di R tapi saya tidak tahu perbedaannya. Tampaknya mereka hanya mengembalikan nilai-p yang sama saat menelepon wilcox.testdengan correct=FALSE, dan wilcox_test(dalam paket koin) dengan distribution="aymptotic". Untuk nilai lain mereka mengembalikan nilai p yang berbeda. Jugawilcox.test selalu mengembalikan W = 0 untuk dataset saya, terlepas dari pengaturan parameternya: …

13 r  hypothesis-testing  p-value  wilcoxon-mann-whitney 

Saya mencoba memodelkan beberapa data menggunakan glmnetpaket dalam R. Katakanlah saya memiliki data berikut training_x <- data.frame(variable1 = c(1, 2, 3, 2, 3), variable2 = c(1, 2, 3, 4, 5)) y <- c(1, 2, 3, 4, 5) (Ini adalah penyederhanaan; data saya jauh lebih rumit.) Lalu saya menggunakan kode berikut …

13 r  glmnet 

Saya mencoba membuat prediksi menggunakan model hutan acak di R. Namun saya mendapatkan kesalahan karena beberapa faktor memiliki nilai yang berbeda di set tes daripada di set pelatihan. Misalnya, suatu faktor Cat_2memiliki nilai 34, 68, 76, dll., Di set tes yang tidak muncul di set pelatihan. Sayangnya, saya tidak memiliki …

13 r  categorical-data  random-forest 

Saya bertanya-tanya bagaimana Anda akan menghasilkan data dari persamaan regresi Poisson di R? Saya agak bingung bagaimana mendekati masalah. Jadi jika saya berasumsi kita memiliki dua prediktor dan X 2 yang didistribusikan N ( 0 , 1 ) . Dan intersepnya adalah 0 dan kedua koefisiennya sama dengan 1. Maka …

13 r  regression  poisson-distribution  simulation 

Sebagai jawaban atas pertanyaan tentang pemilihan model di hadapan multikolinieritas , Frank Harrell menyarankan : Masukkan semua variabel dalam model tetapi jangan menguji efek dari satu variabel yang disesuaikan untuk efek variabel yang bersaing ... Tes chunk dari variabel yang bersaing kuat karena variabel collinear bergabung dalam keseluruhan uji derajat …

13 r  model-selection  multicollinearity 

Ini juga merupakan pertanyaan yang sangat terkait dengan Python sebagai meja kerja statistik dan unggul sebagai meja kerja statistik . Saya tahu ada diskusi besar tentang Ruby versus Python tetapi ini bukan poin dari pertanyaan ini. Saya berpikir bahwa Ruby lebih cepat daripada Python dan memiliki sintaksis yang sangat alami …

13 r  python  software  ruby