«r» etiketlenmiş sorular

Çok boyutlu bir veri setinin iki kısmı var, diyelim ki trainve test. Ve tren veri setine dayalı bir model oluşturmak ve test veri setinde doğrulamak istiyorum. Küme sayısı bilinmektedir. R'de k-ortalamaları kümelemeyi uygulamaya çalıştım ve kümelerin merkezlerini içeren bir nesne aldım: kClust <- kmeans(train, centers=N, nstart=M) R'de bulunan ve test …

14 r  clustering  k-means 

Kruskal-Wallis ve Mann-Whitney U testi için bir güç analizi yapmak mümkün müdür? Evetse, bunu yapan R paketleri / işlevleri var mı?

14 r  nonparametric  power-analysis  kruskal-wallis 

İki bağımsız örneğe uygulanan aynı testten elde edilen sonuçlara sahibim: x <- c(17, 12, 13, 16, 9, 19, 21, 12, 18, 17) y <- c(10, 6, 15, 9, 8, 11, 8, 16, 13, 7, 5, 14) Ve bir Wilcoxon rütbe toplamı testi yapmak istiyorum. istatistiğini elle hesapladığımda şunu elde ederim: …

14 r  wilcoxon-mann-whitney  wilcoxon-signed-rank 

Birisi bana iki oran arasındaki farkın güven aralığını doğru bir şekilde hesaplayıp hesaplamadığımı bilmeme izin veriyor mu? Örnek büyüklüğü 34'tür ve bunlardan 19'u kadın, 15'i erkektir. Bu nedenle, oranlardaki fark 0.1176471'dir. Farkın -0.1183872 ile 0.3536814 arasında olması için% 95 güven aralığını hesaplıyorum. Güven aralığı sıfırdan geçerken, fark istatistiksel olarak anlamlı …

14 r  confidence-interval 

İşte bu soruya ilişkin bağlamım: Anlayabildiğim kadarıyla, ağırlıklı veri ve surveypaketi kullanırken R'de sıradan bir en küçük kareler regresyonunu çalıştıramayız . Burada kullanmak zorundayız svyglm(), bunun yerine genelleştirilmiş doğrusal bir model çalıştırıyor olmalıyız (aynı şey olabilir mi? Burada neyin farklı olduğu konusunda bulanıkım). OLS ve lm()fonksiyon aracılığıyla, yorumunu anladığım bir …

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

Bu iki fonksiyon R'de bulunur, ancak farklılıklarını bilmiyorum. Çağrılırken onlar sadece aynı p-değerleri döndürür görünüyor wilcox.testile correct=FALSEve wilcox_testile (sikke paketinde) distribution="aymptotic". Diğer değerler için farklı p değerleri döndürürler. Ayrıca wilcox.test, parametre ayarlarından bağımsız olarak, veri kümem için her zaman W = 0 döndürüyor: x = c(1, 1, 1, 3, 3, …

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

R'de bir Poisson regresyon denkleminden nasıl veri üreteceğinizi merak ediyordum. Soruna nasıl yaklaşacağım konusunda kafam karıştı. Ben varsayarsak yüzden iki yordayıcılarını var ve X 2 dağıtılır N ( 0 , 1 ) . Ve kesme noktası 0'dır ve her iki katsayı 1'e eşittir. Sonra tahminim basitçe:X1X1X_1X2X2X_2N(0,1)N(0,1)N(0,1) log(Y)=0+1⋅X1+1⋅X2log⁡(Y)=0+1⋅X1+1⋅X2\log(Y) = 0+ 1\cdot …

14 r  regression  poisson-distribution  simulation 

Aşağıdaki PCA biplotunu düşünün: library(mvtnorm) set.seed(1) x <- rmvnorm(2000, rep(0, 6), diag(c(5, rep(1,5)))) x <- scale(x, center=T, scale=F) pc <- princomp(x) biplot(pc) Bir sürü kırmızı ok çizilir, ne anlama geliyorlar? "Var1" ile etiketlenmiş ilk okun veri kümesinin en değişken yönünü göstermesi gerektiğini biliyordum (eğer bunları 2000 veri noktası olarak düşünürsek, …

14 r  pca  linear-algebra  biplot 

Andrew Ng'in Stanford'daki ücretsiz çevrimiçi makine öğrenimi kursundaki videolardan geçiyorum . Degrade İnişi, bunu yapmak için Octave'de lineer regresyon ve yazma fonksiyonlarını çözmek için bir algoritma olarak tartışır. Muhtemelen R'de bu işlevleri yeniden yazabilirim, ama sorum şu ki, lm () işlevi zaten doğrusal regresyonun çıktısını vermiyor mu? Neden kendi degrade …

14 r  regression  machine-learning  gradient-descent 

Birisi bana Cox modelimi düz İngilizce olarak açıklayabilir mi? Bu işlevi kullanarak tüm verilerime aşağıdaki Cox regresyon modelini yerleştirdim cph. Verilerim adlı bir nesneye kaydedilir Data. Değişkenler w, xve ysüreklidir; ziki seviyeli bir faktördür. Zaman ay olarak ölçülür. Bazı hastalarım değişkenleri için veri eksik z( Not : Dr.Harrell'in aşağıda, modelime …

14 r  cross-validation  survival  cox-model 

Birisi bana diğer rastgele etkileri de tahmin etmem gerektiğinde R'nin parçalı doğrusal bir modelde (sabit veya rastgele bir parametre olarak) kırılma noktasını nasıl tahmin edeceğini söyleyebilir mi? Aşağıda bir kırılma noktası için rastgele bir eğim varyansı ve rastgele bir y kesme noktası varyansı ile bir hokey sopası / kırık çubuk …

14 r  mixed-model  lme4-nlme  change-point  piecewise-linear 

Başlık her şeyi söylüyor ve kafam karıştı. Aşağıdaki R tekrarlanan önlemler aov () çalışır ve eşdeğer bir lm () çağrısı olduğunu düşündüğüm çalışır, ancak farklı hata kalıntıları (kareler toplamları aynı olmasına rağmen) döndürür. Açıkça aov () 'dan kalıntılar ve takılmış değerler modelde kullanılan değerlerdir, çünkü karelerin toplamları modelin her birine …

14 r  anova  repeated-measures  linear-model 

R öğrenmeye dalmak üzereyim ve öğrenme projem, tahmini bir denklem geliştirmek için bir veri kümesine karışık veya rastgele etkiler regresyonunun uygulanmasını gerektirecektir. Bu yazıda yazarın endişesini paylaşıyorum Karışık efekt modelleri için nlme veya lme4 R kütüphanesi nasıl seçilir? NLME veya LME4'ün kendimi tanımak için daha iyi bir paket olup olmadığını …

14 r  regression  random-effects-model 

Birisi bana boyuna verileri kullanma aov()ve lme()analiz etme arasındaki farkı söyleyebilir ve bu iki yöntemden elde edilen sonuçları nasıl yorumlayabilir? Aşağıda, aynı kullanarak veri kümesi analiz aov()ve lme()ve 2 farklı sonuçlar aldık. İle aov()ben Tedavi etkileşimi süresinde önemli bir sonuç var, ama doğrusal karma modelin yerleştirilmesi, Tedavi etkileşimi zaman önemsizdir. …

14 r  mixed-model  repeated-measures  panel-data 

Diyelim ki, konum, hız ve ivme için bazı ölçümler içeren bir veri kümem var. Hepsi aynı "koşu" dan geliyor. Doğrusal bir sistem kurabilirim ve tüm bu ölçümlere bir polinom sığdırabilirim. Ama aynısını spline ile yapabilir miyim? Bunu yapmanın 'R' yolu nedir? İşte sığdırmak istediğim bazı simüle veriler: f <- function(x) …

14 r  model  fitting  splines  derivative