«r» etiketlenmiş sorular

Bir standart hata tahmini standart sapma σ ( θ ) kestiricinin bir parametre için .σ^( θ^)σ^(θ^)\hat \sigma(\hat\theta)θ^θ^\hat\thetaθθ\theta Kalıntıların tahmini standart sapmasına neden summary.lm"artık standart sapma" değil , "artık standart hata" (örn. R fonksiyonunun çıktısında) denir ? Burada hangi parametre tahminini standart bir hatayla donatıyoruz? Her bir kalıntıyı "kendi" hata terimi …

14 r  standard-error  residuals  terminology 

Birden fazla eksojen değişkenli çoklu zaman serisi ARIMA modeline dayalı bir tahmin yapmak istiyorum. Ne istatistiklere ne de UR'ye bağlı kalmak konusunda o kadar yetenekli olmadığım için mümkün olduğunca basit (3 aylık trend tahmini yeterlidir). Ben 1 bağımlı zaman serisi ve 3-5 tahmin zaman serisi, tüm aylık veri, boşluk, aynı …

14 r  time-series  arima 

Düzgün kullandığım bazı verilerim var loess. Düzeltilmiş çizginin bükülme noktalarını bulmak istiyorum. Mümkün mü? Eminim birisi bunu çözmek için süslü bir yöntem yapmış ... Yani ... sonuçta, R! Kullandığım yumuşatma işlevini değiştirme konusunda iyiyim. Sadece kullandım loessçünkü geçmişte kullandığım şey buydu. Ancak herhangi bir yumuşatma işlevi iyidir. Bükülme noktalarının kullandığım …

14 r  smoothing  loess 

Ben bir ROC eğrisi (hassasiyet ve özgüllük en üst düzeye çıkarıldığı değeri) için en uygun kesme noktası hesaplamak anlamaya çalışıyorum. aSAHPaketten veri kümesini kullanıyorum pROC. outcomeDeğişken iki bağımsız değişkenler tarafından açıklanabilir: s100bve ndka. EpiPaketin sözdizimini kullanarak iki model oluşturdum: library(pROC) library(Epi) ROC(form=outcome~s100b, data=aSAH) ROC(form=outcome~ndka, data=aSAH) Çıktı aşağıdaki iki grafikte gösterilmektedir: …

14 r  roc  sensitivity-analysis  sensitivity-specificity 

R dilinde yeniyim. Regresyonun dört varsayımını da yerine getiren çoklu doğrusal regresyon modelinden nasıl simüle edileceğini bilmek istiyorum. tamam teşekkürler. Diyelim ki bu veri kümesine dayalı olarak verileri simüle etmek istiyorum: y<-c(18.73,14.52,17.43,14.54,13.44,24.39,13.34,22.71,12.68,19.32,30.16,27.09,25.40,26.05,33.49,35.62,26.07,36.78,34.95,43.67) x1<-c(610,950,720,840,980,530,680,540,890,730,670,770,880,1000,760,590,910,650,810,500) x2<-c(1,1,3,2,1,1,3,3,2,2,1,3,3,2,2,2,3,3,1,2) fit<-lm(y~x1+x2) summary(fit) sonra çıktıyı alıyorum: Call: lm(formula = y ~ x1 + x2) Residuals: Min 1Q …

14 r  multiple-regression  simulation 

Şu anda Gamma dağıtımını kullanarak bir GLM takmak için R sözdizimini anlamakta bir sorunum var. Her satır 3 ortak değişken ( X1,X2,X3X1,X2,X3X_1, X_2, X_3 ), bir yanıt değişkeni ( YYY ) ve bir şekil parametresi ( KKK ) içeren bir veri kümesi var . Gamma dağılımının ölçeğini 3 değişkenin doğrusal …

14 r  generalized-linear-model  gamma-distribution  dglm 

Kareler (WSS) grup içi toplamını en aza indirmek için birimleri kümeleri halinde kümelemem gerekiyor , ancak kümelerin her birinin en az m birimi içermesini sağlamalıyım . R'nin kümeleme işlevlerinden herhangi biri , minimum küme boyutu kısıtlamasına tabi k kümeleri içinde kümelemeye izin veriyorsa herhangi bir fikir var mı? kmeans () …

14 r  clustering 

Ben tam olarak coxph "coef" ve "(exp) coef" çıktı anlamını ayırt etmeye çalışıyorum. Görünüşe göre "(exp) coef", modeldeki ilk değişkenin komutta atanan gruba göre karşılaştırmalarıdır. Coxph işlevi "coef" ve "(exp) coef" değerlerine nasıl ulaşır? Ek olarak, sansür söz konusu olduğunda coxph bu değerleri nasıl belirler?

14 r  survival  interpretation 

Ben dput kullanarak burada yüklediğim sayıların bir vektör (... / code / MyData.Rdata) var. Ben bca ci almak istiyorum bu yüzden bu kodu yazdım: my.mean <- function(dat, idx){ return (mean(dat[idx], na.rm = TRUE)) } boot.out<-boot(data=my.data, statistic = my.mean, R=1000) Ama aşağıdakileri çalıştırdığımda bunu alıyorum: > boot.ci(boot.out) Error in bca.ci(boot.out, conf, …

14 r  bootstrap 

hayatta kalma süresinin bir örneğine dayanarak , Kaplan-Meier tahmincisini kullanarak, belirli bir için , hayatta kalma süresinin olasılığını tahmin etmek istiyorum . Bunu yapmak mümkün mü ? Lütfen mutlaka bir etkinlik zamanı olmadığını unutmayın .nnnttttttRttt

14 r  kaplan-meier 

Ben sadece negatif bir binom GLM koştu ve bu çıktı: Call: glm.nb(formula = small ~ method + site + depth, data = size.dat, init.theta = 1.080668549, link = log) Deviance Residuals: Min 1Q Median 3Q Max -2.2452 -0.9973 -0.3028 0.3864 1.8727 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 1.6954 …

14 r  categorical-data  generalized-linear-model  interpretation  negative-binomial 

Aşağıdaki örnekte > m = matrix(c(3, 6, 5, 6), nrow=2) > m [,1] [,2] [1,] 3 5 [2,] 6 6 > (OR = (3/6)/(5/6)) #1 [1] 0.6 > fisher.test(m) #2 Fisher's Exact Test for Count Data data: m p-value = 0.6699 alternative hypothesis: true odds ratio is not equal to …

14 r  odds-ratio  fishers-exact 

Bir senaryo kullanıyorum. Çekirdek kayıtlar içindir. Belirli bir derinlikte (ilk sütunda) sütunlarda farklı elementel kompozisyonlar gösteren bir veri çerçevesi var. Onunla bir PCA yapmak istiyorum ve seçmem gereken standardizasyon yöntemi hakkında kafam karıştı. Sizden kimse kullandı clr()için verilerinizi hazırlamak için prcomp()? Yoksa çözümlerimi karıştırıyor mu? Öznitelik ölçeğini kullanmaya ek olarak …

13 r  pca  normalization  compositional-data 

Douglas Bates, aşağıdaki modellerin "vektör değerli rasgele efektler için varyans-kovaryans matrisinin bileşik simetri adı verilen özel bir formu varsa" eşdeğer olduğunu belirtir ( bu sunumdaki slayt 91 ): m1 <- lmer(y ~ factor + (0 + factor|group), data) m2 <- lmer(y ~ factor + (1|group) + (1|group:factor), data) Özellikle Bates …

13 r  anova  mixed-model  repeated-measures  lme4-nlme 

Regresyon analizinin yorumlanması üzerinde çok az arka plana sahibim ama r, r kare ve artık standart sapmanın anlamı hakkında gerçekten kafam karıştı. Tanımları biliyorum: Karakterizasyonları r dağılım grafiğindeki iki değişken arasındaki doğrusal ilişkinin gücünü ve yönünü ölçer R-kare, verilerin yerleştirilmiş regresyon hattına ne kadar yakın olduğunun istatistiksel bir ölçüsüdür. Kalıntı …

13 r  regression  regression-coefficients  linear  pearson-r