Pertanyaan yang diberi tag «r»

Lee dan Lemieux (hlm. 31, 2009) menyarankan peneliti untuk menyajikan grafik saat melakukan analisis desain diskontinuitas Regresi (RDD). Mereka menyarankan prosedur berikut: "... untuk beberapa bandwidth , dan untuk beberapa jumlah nampan K 0 dan K 1 ke kiri dan kanan dari nilai cutoff, masing-masing, idenya adalah untuk membangun nampan …

10 r  regression  data-visualization  stata  regression-discontinuity 

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, cor=TRUE)

10 r  pca 

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 rumus berikut:indep_var c dep_var ar(1) EViews mengklaim bahwa …

10 r  autoregressive  software 

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 pemasangan model seperti seperti membentuk kembali beberapa variabel respon menjadi satu variabel format panjang …

10 r  mixed-model  multivariate-analysis  lme4-nlme  multivariate-regression 

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 summary(model). Dua yang terakhir ini memberikan jawaban yang sama. Menggunakan banyak data yang …

10 r  anova  glmm  r  mixed-model  bootstrap  sample-size  cross-validation  roc  auc  sampling  stratification  random-allocation  logistic  stata  interpretation  proportion  r  regression  multiple-regression  linear-model  lm  r  cross-validation  cart  rpart  logistic  generalized-linear-model  econometrics  experiment-design  causality  instrumental-variables  random-allocation  predictive-models  data-mining  estimation  contingency-tables  epidemiology  standard-deviation  mean  ancova  psychology  statistical-significance  cross-validation  synthetic-data  poisson-distribution  negative-binomial  bioinformatics  sequence-analysis  distributions  binomial  classification  k-means  distance  unsupervised-learning  euclidean  correlation  chi-squared  spearman-rho  forecasting  excel  exponential-smoothing  binomial  sample-size  r  change-point  wilcoxon-signed-rank  ranks  clustering  matlab  covariance  covariance-matrix  normal-distribution  simulation  random-generation  bivariate  standardization  confounding  z-statistic  forecasting  arima  minitab  poisson-distribution  negative-binomial  poisson-regression  overdispersion  probability  self-study  markov-process  estimation  maximum-likelihood  classification  pca  group-differences  chi-squared  survival  missing-data  contingency-tables  anova  proportion 

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 MeanDecreaseGini Sepal.Length 0.027078501 0.019418330 0.040497602 0.02898837 9.173648 Sepal.Width 0.008553449 0.001962036 0.006951771 0.00575489 2.472105 Petal.Length 0.313303381 0.291818815 0.280981959 0.29216790 41.284869 Petal.Width 0.349686983 0.318527008 0.270975757 0.31054451 46.323415 > …

10 r  random-forest  importance 

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 …

10 r  mixed-model  repeated-measures  model-selection  lme4-nlme 

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 ), tetapi:α=.05,n=10,m=5α=.05,n=10,m=5\alpha=.05, n = 10, m = …

10 r  hypothesis-testing  nonparametric  wilcoxon-mann-whitney 

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 3.1626716539637488e + 90 - angka yang …

10 r  python  numpy  scipy 

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 …

10 r  anova  contrasts  generalized-least-squares 

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 satu instrumen memimpin yang lain. Saya punya data untuk kedua …

10 r  time-series  cross-correlation  unevenly-spaced-time-series 

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 nlmeIIRC. Jelas, pemahaman saya tentang apa fungsi …

10 r  lme4-nlme  polynomial 

Saya mencoba untuk mereplikasi karya seorang kolega dan saya memindahkan analisis dari Stata ke R. Model yang dia gunakan memanggil opsi "cluster" dalam fungsi nbreg untuk mengelompokkan kesalahan standar. Lihat http://repec.org/usug2007/crse.pdf untuk deskripsi yang cukup lengkap tentang apa dan mengapa opsi ini Pertanyaan saya adalah bagaimana menjalankan opsi yang sama …

10 r  stata  negative-binomial 

Jika saya mencocokkan data saya dengan sesuatu seperti lm(y~a*b), dalam sintaks R, di mana avariabel biner dan bvariabel numerik, maka a:bistilah interaksi adalah perbedaan antara kemiringan pada y~bat a= 0 dan at a= 1. Sekarang, katakanlah hubungan antara ydan bmelengkung. Jika saya sekarang cocok lm(y~a*poly(b,2)), maka a:poly(b,2)1perubahan dalam y~bkondisi bersyarat …

10 r  multiple-regression  splines 

Ada banyak metode untuk melakukan regularisasi - , , dan regularisasi berbasis misalnya. Menurut Friedman Hastie & Tibsharani , pembuat peraturan terbaik tergantung pada masalah: yaitu sifat fungsi target yang sebenarnya, dasar tertentu yang digunakan, rasio sinyal terhadap kebisingan, dan ukuran sampel.L 1 L 2L.0L0L_0L.1L1L_1L.2L2L_2 Adakah penelitian empiris yang membandingkan …

10 r  regression  machine-learning  regularization