Pertanyaan yang diberi tag «r»

Mari kita memiliki beberapa model linier, misalnya ANOVA hanya sederhana: # data generation set.seed(1.234) Ng <- c(41, 37, 42) data <- rnorm(sum(Ng), mean = rep(c(-1, 0, 1), Ng), sd = 1) fact <- as.factor(rep(LETTERS[1:3], Ng)) m1 = lm(data ~ 0 + fact) summary(m1) Hasilnya adalah sebagai berikut: Call: lm(formula = …

18 r  regression  confidence-interval 

Ini adalah bingkai data saya: Group <- c("G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G1","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G2","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3","G3") Subject <- c("S1","S2","S3","S4","S5","S6","S7","S8","S9","S10","S11","S12","S13","S14","S15","S1","S2","S3","S4","S5","S6","S7","S8","S9","S10","S11","S12","S13","S14","S15","S1","S2","S3","S4","S5","S6","S7","S8","S9","S10","S11","S12","S13","S14","S15") Value <- c(9.832217741,13.62390117,13.19671612,14.68552076,9.26683366,11.67886655,14.65083473,12.20969772,11.58494621,13.58474896,12.49053635,10.28208078,12.21945867,12.58276212,15.42648969,9.466436017,11.46582655,10.78725485,10.66159358,10.86701127,12.97863424,12.85276916,8.672953949,10.44587257,13.62135205,13.64038394,12.45778874,8.655142642,10.65925259,13.18336949,11.96595556,13.5552118,11.8337142,14.01763101,11.37502161,14.14801305,13.21640866,9.141392359,11.65848845,14.20350364,14.1829714,11.26202565,11.98431285,13.77216009,11.57303893) data <- data.frame(Group, Subject, Value) Lalu saya menjalankan model efek linear-campuran untuk membandingkan perbedaan 3 Grup pada "Nilai", di mana "Subjek" adalah faktor acak: library(lme4) library(lmerTest) model <- lmer (Value~Group + (1|Subject), data = data) summary(model) …

18 r  lme4-nlme  post-hoc 

Konteks : Regresi hierarkis dengan beberapa data yang hilang. Pertanyaan : Bagaimana cara saya menggunakan estimasi kemungkinan maksimum informasi penuh (FIML) untuk mengatasi data yang hilang dalam R? Apakah ada paket yang akan Anda rekomendasikan, dan apa langkah-langkah khasnya? Sumber dan contoh online akan sangat membantu juga. PS : Saya …

18 r  maximum-likelihood  missing-data 

Saya memiliki data tentang perilaku pergerakan (waktu yang dihabiskan untuk tidur, kurang gerak, dan melakukan aktivitas fisik) yang jumlahnya kira-kira 24 (seperti dalam jam per hari). Saya ingin membuat variabel yang menangkap waktu relatif yang dihabiskan dalam masing-masing perilaku ini - Saya telah diberitahu bahwa transformasi log-rasio isometrik akan mencapai …

17 r  multivariate-analysis  data-transformation  compositional-data 

Saya telah belajar tentang model bahaya proporsional Cox. Saya memiliki banyak pengalaman pas model regresi logistik, dan model sehingga untuk membangun intuisi saya telah membandingkan cocok menggunakan coxphdari R "survival" dengan model regresi logistik cocok menggunakan glmdengan family="binomial". Jika saya menjalankan kode: library(survival) s = Surv(time=lung$time, event=lung$status - 1) summary(coxph(s …

17 r  logistic  survival  cox-model  power 

Seseorang dapat melakukan regresi logit di R menggunakan kode seperti: > library(MASS) > data(menarche) > glm.out = glm(cbind(Menarche, Total-Menarche) ~ Age, + family=binomial(logit), data=menarche) > coefficients(glm.out) (Intercept) Age -21.226395 1.631968 Sepertinya algoritme pengoptimalan telah terkonvergensi - ada informasi tentang jumlah langkah algoritme penilaian skor: Call: glm(formula = cbind(Menarche, Total - …

17 r  generalized-linear-model  optimization  algorithms  logit 

Saya memiliki beberapa data yang saya pas menggunakan model LOESS di R, memberi saya ini: Data memiliki satu prediktor dan satu respons, dan bersifat heteroskedastik. Saya juga menambahkan interval kepercayaan. Masalahnya adalah intervalnya adalah interval kepercayaan untuk garis, sedangkan saya tertarik pada interval prediksi. Misalnya, panel bawah lebih bervariasi daripada …

17 r  regression  prediction-interval  loess 

Ini adalah posting pertama saya di StackExchange, tetapi saya telah menggunakannya sebagai sumber daya selama beberapa waktu, saya akan melakukan yang terbaik untuk menggunakan format yang sesuai dan melakukan pengeditan yang sesuai. Juga, ini adalah pertanyaan multi-bagian. Saya tidak yakin apakah saya harus membagi pertanyaan menjadi beberapa posting berbeda atau …

17 r  logistic  model-selection  regression-strategies 

Saya bisa menggunakan beberapa panduan tentang penyajian beberapa data. Plot pertama ini adalah perbandingan case-control untuk sitokin IL-10. Saya secara manual mengatur sumbu y untuk memasukkan 99% data. Alasan saya mengatur ini secara manual adalah karena grup kasus memiliki pencilan yang ekstrim. Kolaborator saya ragu untuk melakukan penghapusan pencilan ke …

17 r  data-visualization  outliers  boxplot  presentation 

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 yang tersedia untuk itu?

17 r  normal-distribution  matlab  simulation  random-generation 

Saya menjalankan regresi logistik ordinal ini di R: mtcars_ordinal <- polr(as.factor(carb) ~ mpg, mtcars) Saya mendapatkan ringkasan model ini: summary(mtcars_ordinal) Re-fitting to get Hessian Call: polr(formula = as.factor(carb) ~ mpg, data = mtcars) Coefficients: Value Std. Error t value mpg -0.2335 0.06855 -3.406 Intercepts: Value Std. Error t value 1|2 …

17 r  regression  logistic  interpretation  ordered-logit 

Saya tidak mengerti penjelasan dalam Rfile bantuan untuk efek () : Untuk model linier yang dipasang oleh lmatau aov, pengaruhnya adalah nilai derajat kebebasan tunggal yang tidak berkorelasi yang diperoleh dengan memproyeksikan data ke subruang ortogonal yang dihasilkan oleh dekomposisi QR selama proses pemasangan. Adakah yang bisa menjelaskan apa artinya …

17 r  regression 

Saya memiliki variabel independen yang disebut "kualitas"; variabel ini memiliki 3 modalitas respons (kualitas buruk; kualitas sedang; kualitas tinggi). Saya ingin memperkenalkan variabel independen ini ke dalam regresi linier berganda. Ketika saya memiliki variabel independen biner (variabel dummy, saya dapat kode 0/ 1) mudah untuk memperkenalkannya ke dalam model regresi …

17 r  regression  categorical-data  multicollinearity  categorical-encoding 

Saya memiliki histogram data jumlah berikut. Dan saya ingin mencocokkan distribusi diskrit dengannya. Saya tidak yakin bagaimana saya harus melakukan ini. Haruskah saya pertama menempatkan superimpose distribusi diskrit, katakanlah distribusi Binomial Negatif, pada histogram sehingga saya akan mendapatkan parameter distribusi diskrit dan kemudian menjalankan tes Kolmogorov-Smirnov untuk memeriksa nilai-p? Saya …

17 r  poisson-distribution  discrete-data  computational-statistics  negative-binomial 

Saya ingin melakukan analisis data ANCOVA mengenai kepadatan epifit tanaman. Pada awalnya, saya ingin tahu apakah ada perbedaan dalam kepadatan tanaman antara dua lereng, satu N dan satu S, tetapi saya memiliki data lain seperti ketinggian, keterbukaan kanopi dan ketinggian tanaman inang. Saya tahu bahwa kovariat saya haruslah dua lereng …

17 r  ancova