«r» etiketlenmiş sorular

R'de, lm()işlevi şu şekilde çağırırsam : lm.1 = lm(response ~ var1 + var2 + var1 * var2) summary(lm.1) Bu bana olan tepkisi değişkenin doğrusal modelini verir var1, var2ve bunlar arasındaki etkileşimi. Ancak, etkileşim terimini tam olarak nasıl yorumlarız? Belgeler bunun var1ve arasındaki "çapraz" olduğunu söylüyor var2, ancak "çapraz" ın tam …

9 r  regression 

Bu yüzden R'deki lojistik regresyon modelleri ile çalışıyorum. İstatistiklere hala yeniyim, ancak şimdiye kadar regresyon modelleri için biraz anlayışım var gibi hissediyorum, ama yine de beni rahatsız eden bir şey var: Bağlantılı resme baktığımda, oluşturduğum örnek bir model için özet R baskılarını görüyorsunuz. Model veri kümesindeki email refound veya olmasın …

9 r  regression  logistic  p-value  interpretation 

ARIMA modelleri ile öngörüye adım attığım için, ARIMA'nın mevsimsellik ve sapmaya uygun bir tahminini nasıl geliştirebileceğimi anlamaya çalışıyorum. Verilerim aşağıdaki zaman serisidir (3 yıldan fazla, açık eğilim yukarı ve görünür mevsimsellik ile, 12, 24, 36 gecikmelerde otokorelasyon tarafından desteklenmiyor gibi görünüyor). > bal2sum3years.ts Jan Feb Mar Apr May Jun Jul …

9 r  time-series  forecasting  arima  seasonality 

benzeyen bir regresyon modelim var:Y=β0+β1X1+β2X2+β3X3+β12X1X2+β13X1X3+β123X1X2X3Y=β0+β1X1+β2X2+β3X3+β12X1X2+β13X1X3+β123X1X2X3Y = \beta_0+\beta_1X_1 + \beta_2X_2 + \beta_3X_3 +\beta_{12}X_1X_2+\beta_{13}X_1X_3+\beta_{123}X_1X_2X_3 ... veya R gösterimi ile: y ~ x1 + x2 + x3 + x1:x2 + x1:x3 + x1:x2:x3 Diyelim ki ve kategorik değişkenler ve sayısal. Karmaşık üç seviye ve standart kontrastlar yerine, test etmek gerekir:X1X1X_1X2X2X_2X3X3X_3X1X1X_1X1 a,X1 b,X1 cX1a,X1b,X1cX_{1a}, …

9 r  contrasts 

Kimse kmeans kullanarak kümelenmenin uygun olup olmadığını belirlemek için iyi bir yöntem biliyor mu? Yani, numuneniz aslında homojen ise? Bir karışım modeli gibi bir şey biliyorum (R'de mclust aracılığıyla) 1: k küme durumu için uygun istatistikler sağlayacaktır, ancak kmeans değerlendirmek için tüm teknikler en az 2 küme gerektirir gibi görünüyor. …

9 r  clustering  k-means 

Lojistik regresyon için bazı test verileri oluşturmaya çalışıyordum ve bu yazıyı buldum Lojistik regresyon için yapay veriler nasıl simüle edilir? Güzel bir cevap ama sadece sürekli değişkenler yaratıyor. Bağlantıdaki ile aynı örnek için y ile ilişkili 5 seviyeli (ABCDE) kategorik bir x3 değişkenine ne dersiniz?

9 r  logistic  simulation 

Matris çarpımı kullanarak Jaccard katsayısını hesaplamanın olası bir yolu olup olmadığını bilmek istiyorum. Bu kodu kullandım jaccard_sim <- function(x) { # initialize similarity matrix m <- matrix(NA, nrow=ncol(x),ncol=ncol(x),dimnames=list(colnames(x),colnames(x))) jaccard <- as.data.frame(m) for(i in 1:ncol(x)) { for(j in i:ncol(x)) { jaccard[i,j]= length(which(x[,i] & x[,j])) / length(which(x[,i] | x[,j])) jaccard[j,i]=jaccard[i,j] } } …

9 r  matrix  binary-data  association-measure  similarities 

Son zamanlarda bir R betiğinde hata ayıkladım ve çok garip bir şey buldum, yazar kendi p-değeri işlevini tanımladı pval <- function(x, y){ if (x+y<20) { # x + y is small, requires R.basic p1<- nChooseK(x+y,x) * 2^-(x+y+1); p2<- nChooseK(x+y,y) * 2^-(x+y+1); pvalue = max(p1, p2) } else { # if …

9 r  hypothesis-testing  p-value 

Aralık sansürlenmiş hayatta kalma süreleri göz önüne alındığında, aralık sansürlü Cox PH modelini nasıl yapabilirim R? Bir rseek araması intcoxartık Rdepoda bulunmayan paketi açar . Neredeyse paketteki coxphfonksiyonun survivalsansürlenmiş sağkalım verilerini aralıksız işleyemediğine inanıyorum. Ayrıca, verileri çarpıtmak ve sonra coxphişlevi kullanmak istemiyorum . Bu yöntem, aralık sansürünün belirsizliğini göz ardı …

9 r  survival  cox-model  interval-censoring 

Ben acemiyim R, ses tahmini R tahmin paketinde nasıl kullanılacağını açıklar mısınız? İlk dönem sayısını ve düzeltme sabitini seçmek istiyorum. d <- c(3,4,41,10,9,86,56,20,18,36,24,59,82,51,31,29,13,7,26,19,20,103,141,145,24,99,40,51,72,58,94,78,11,15,17,53,44,34,12,15,32,14,15,26,75,110,56,43,19,17,33,26,40,42,18,24,69,18,18,25,86,106,104,35,43,12,4,20,16,8) 70 dönemim var, başlangıç ​​için 40 Dönem ve örneklem dışı için 30 Dönem kullanmak istiyorum. ses(d, h=30, level=c(80,95), fan=FALSE,initial=c("simple"), alpha=.1) Doğru mu?

9 r  time-series  forecasting 

Sadece bir öngörücü (örneğin (x, y)) ile bazı verilere doğrusal bir model sığdırmaya çalışıyorum. Veriler, küçük x değerleri için, y değerlerinin düz bir çizgiye sıkıca oturduğu, ancak x değerleri arttıkça y değerlerinin daha uçucu hale geleceği şekildedir. İşte bu tür verilere bir örnek (R kodu) y = c(3.2,3.4,3.5,3.8,4.2,5.5,4.5,6.8,7.4,5.9) x = …

9 r  regression  variance  data-transformation 

Aşağıdaki veri kümesi var: https://dl.dropbox.com/u/22681355/ORACLE.csv ve 'Tarih' ile 'Aç' günlük değişiklikleri çizmek istiyorum, bu yüzden aşağıdakileri yaptım: oracle <- read.csv(file="http://dl.dropbox.com/u/22681355/ORACLE.csv", header=TRUE) plot(oracle$Date, oracle$Open, type="l") ve aşağıdakileri alıyorum: Şimdi bu açıkça en güzel arsa değil, bu yüzden böyle ayrıntılı verileri çizerken kullanmak için doğru yöntemin ne olduğunu merak ediyorum?

9 r  data-visualization 

Tamam, ben lineer regresyonu anlamaya çalışıyorum. Bir veri setim var ve her şey yolunda görünüyor, ama kafam karıştı. Bu benim doğrusal model özeti: Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.2068621 0.0247002 8.375 4.13e-09 *** temp 0.0031074 0.0004779 6.502 4.79e-07 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 …

9 r  regression 

İki grup arasındaki insidans oranlarını karşılaştırmak istiyorum (biri hastalığı olmayan ve diğeri ile). İnsidans oranı oranını (IRR), yani insidans oranı grup B / insidans oranı grup A'yı hesaplamayı ve sonra bu oranın 1'e eşit olup olmadığını test etmeyi ve son olarak IRR için% 95 CI aralıklarını hesaplamayı planlıyordum. Bir kitapta% …

9 r  poisson-distribution  epidemiology  incidence-rate-ratio 

İstatistikler hakkında fazla bir şey bilmiyorum, bu yüzden bana katlan. Diyelim ki 1000 kişilik bir işçim var. En zor çalışanın kim olduğunu bulmak istiyorum ama sadece bir saat değerinde 1-100 kişilik gruplar halinde yapılan iş miktarını ölçebilirim. Her bir işçinin her zaman aynı miktarda çalışma yaptığını varsayarsak, çok sayıda deneme …

9 r  regression  data-mining