«r» etiketlenmiş sorular

Bir lmer () modelinden bir tahmin çevresinde bir tahmin aralığı almak istiyorum. Bununla ilgili bazı tartışmalar buldum: http://rstudio-pubs-static.s3.amazonaws.com/24365_2803ab8299934e888a60e7b16113f619.html http://glmm.wikidot.com/faq fakat rastgele etkilerin belirsizliğini hesaba katmamış gibi görünüyorlar. İşte özel bir örnek. Altın balıkla yarışıyorum. Son 100 yarış hakkında verilerim var. RE tahminlerimin ve FE tahminlerinin belirsizliğini dikkate alarak 101'ini tahmin …

37 r  mixed-model  prediction  prediction-interval  lme4-nlme 

Baltık Denizi'nin deniz yüzeyinin sıcaklığının her yıl aynı olduğunu varsaymak istiyorum, sonra bunu bir fonksiyon / doğrusal modelle tanımlamak istiyorum. Benim düşüncem, yılı bir ondalık sayı olarak girmek (ya da num_months / 12) ve sıcaklığın o zaman ne olacağı hakkında bilgi almaktı. R'deki lm () fonksiyonuna atıldığında, sinüzoidal verileri tanımıyor, …

37 r  regression  time-series  lm 

Tedavi değişkenini ( Curevs. No Cure) aldıktan sonra sonuç değişkeninin iyileştirildiği lojistik bir regresyon oluşturdum . Bu çalışmada tüm hastalar tedavi gördü. Diyabetin bu sonuçla ilişkili olup olmadığını görmekle ilgileniyorum. R'de lojistik regresyon çıktım şöyle gözüküyor: Call: glm(formula = Cure ~ Diabetes, family = binomial(link = "logit"), data = All_patients) …

37 r  hypothesis-testing  logistic  generalized-linear-model  odds-ratio 

İki farklı fonksiyon ( prcompve princomp) kullanarak R ile temel bileşen analizi (PCA) yaptım ve PCA puanlarının işaret olarak farklı olduğunu gözlemledim. Nasıl olabilir? Bunu düşün: set.seed(999) prcomp(data.frame(1:10,rnorm(10)))$x PC1 PC2 [1,] -4.508620 -0.2567655 [2,] -3.373772 -1.1369417 [3,] -2.679669 1.0903445 [4,] -1.615837 0.7108631 [5,] -0.548879 0.3093389 [6,] 0.481756 0.1639112 [7,] 1.656178 …

37 r  pca  factor-analysis 

Bu sorunun zaten burada bir cevabı var: Rasgele sayı üretecindeki bir tohum tam olarak nedir? 3 cevaplar Birinin set.seed()R'de sahte rastgele sayı üretimi için kullandığını fark ediyorum . Ayrıca aynı sayıyı kullanarak, set.seed(123)sonuçları yeniden üretebileceğiniz sigortaları da anlıyorum . Ama anlamadığım, değerlerin kendileri ne anlama geldiğidir. Birkaç fonksiyonla oynuyorum ve …

37 r  random-generation 

Verilerimi çeşitli modellere sığdırmaya çalışıyordum ve fitdistrkütüphanedeki işlevin bana en iyi MASSsonucu Rverdiğini anladım Negative Binomial. Şimdi wiki sayfasından, tanım şöyle verilir: NegBin (r, p) dağılımı, son denemede başarı ile sonuçlanan k + r Bernoulli (p) denemelerinde k başarısızlık ve r başarı olasılığını tanımlar. RModel uydurma gerçekleştirmek için kullanmak bana …

37 r  distributions  modeling  negative-binomial 

Sınıflandırma ve regresyon için tahmine dayalı modeller oluşturmak için R'deki şapka paketini kullanıyorum . Caret, çapraz doğrulama veya önyükleme kayışı yoluyla model hiper parametrelerini ayarlamak için birleştirilmiş bir arayüz sağlar. Örneğin, sınıflandırma için basit bir 'en yakın komşular' modeli oluşturuyorsanız, kaç tane komşu kullanmalısınız? 2? 10? 100? Caret, verilerinizi yeniden …

37 r  time-series  cross-validation 

Şu anda lme4 R paketini kullanıyorum . Rasgele efektli doğrusal karışık efekt modellerini kullanıyorum: library(lme4) mod1 <- lmer(r1 ~ (1 | site), data = sample_set) #Only random effects mod2 <- lmer(r1 ~ p1 + (1 | site), data = sample_set) #One fixed effect + # random effects mod3 <- lmer(r1 …

37 r  regression  mixed-model  random-effects-model 

EffectsPaket için çok hızlı ve kolay bir yol sağlar karışık etki modeli sonuçlar doğrusal çizilmesi ile elde edilen lme4bir paket . effectFonksiyon hesaplar güven aralığı (GA) çok hızlı bir şekilde, ama nasıl güvenilir bunlar güven aralıkları nelerdir? Örneğin: library(lme4) library(effects) library(ggplot) data(Pastes) fm1 <- lmer(strength ~ batch + (1 | …

36 r  mixed-model  confidence-interval  effects  lme4-nlme 

Durumum için hangi çapraz doğrulama yönteminin en iyisi olduğunu bulmaya çalışıyorum. Aşağıdaki veriler sadece (R cinsinden ) konu üzerinde çalışmak için bir örnektir, fakat benim gerçek Xverilerim ( xmat) birbirleriyle ve ydeğişken ( ymat) değişkeniyle farklı derecelerle ilişkilendirilmektedir . R kodu verdim, ancak sorum R ile ilgili değil, yöntemlerle ilgili. …

36 r  regression  cross-validation  linear-model 

Sahip olduğum bazı verilere ikinci dereceden bir polinom uyumu yaratmaya çalışıyorum. Diyelim ki bu uyumu şöyle çizdim ggplot(): ggplot(data, aes(foo, bar)) + geom_point() + geom_smooth(method="lm", formula=y~poly(x, 2)) Alırım: Bu yüzden, ikinci dereceden bir form uyumu oldukça iyi çalışıyor. R ile hesaplarım: summary(lm(data$bar ~ poly(data$foo, 2))) Ve anladım: lm(formula = data$bar …

36 r  regression  interpretation  regression-coefficients 

Eşit olmayan örneklem büyüklükleri ( , , ) içeren 3 grup yosun biyokütlesinden ( , , ) verilerim var ve bu grupların aynı popülasyondan olup olmadıklarını karşılaştırmak istiyorum.AAAC n A = 15 n B = 13 n C = 12BBBCCCnA=15nA=15n_A=15nB=13nB=13n_B=13nC=12nC=12n_C=12 Tek yönlü ANOVA kesinlikle gidilecek yol olacaktır, ancak verilerim üzerinde …

36 r  anova  data-transformation  heteroscedasticity 

plot(density(rexp(100)) Açıkçası sıfırın solundaki tüm yoğunluk yanlılığı temsil ediyor. İstatist olmayanlar için bazı verileri özetlemeye çalışıyorum ve negatif olmayan verilerin neden sıfıra sol tarafa yoğunluğa sahip olduğuna dair sorulardan kaçınmak istiyorum. Grafikler randomizasyon kontrolü içindir; Değişkenlerin tedavi ve kontrol gruplarına göre dağılımını göstermek istiyorum. Dağılımlar genellikle üsteldir. Histogramlar çeşitli nedenlerden …

36 r  pdf  gamma-distribution  kernel-smoothing 

50.000'den fazla değişkenli bir modelde Kement veya ridge regresyon kullanmak istiyorum. Bunu R'deki yazılım paketini kullanarak yapmak istiyorum. Büzülme parametresini ( ) nasıl tahmin edebilirim ?λλ\lambda Düzenlemeler: İşte anladığım nokta: set.seed (123) Y <- runif (1000) Xv <- sample(c(1,0), size= 1000*1000, replace = T) X <- matrix(Xv, nrow = 1000, …

36 r  lasso  ridge-regression  high-dimensional 

Amazon EC2 üzerine kurulu kümelerde R, python veya oktav betiklerini çalıştırmak için bir platform sağlayan küçük bir şirket listesine sahibim. Eklemem gereken başka isimler var mı? Cloudnumbers Opani crdata

36 r  machine-learning