«r» etiketlenmiş sorular

Verilerim, istihdam edilen nüfusun zaman dizisi, L ve zaman aralığı, yıl. n.auto=auto.arima(log(L),xreg=year) summary(n.auto) Series: log(L) ARIMA(2,0,2) with non-zero mean Coefficients: ar1 ar2 ma1 ma2 intercept year 1.9122 -0.9567 -0.3082 0.0254 -3.5904 0.0074 s.e. NaN NaN NaN NaN 1.6058 0.0008 sigma^2 estimated as 1.503e-06: log likelihood=107.55 AIC=-201.1 AICc=-192.49 BIC=-193.79 In-sample error …

9 r  regression  arima 

@whuber, bir zaman noktası için çok değişkenli sonuçların ( , ve ) nasıl simüle edileceğini göstermiştir .y1y1y_1y2y2y_2y3y3y_3 Bildiğimiz gibi, boyuna veriler genellikle tıbbi çalışmalarda ortaya çıkar. Sorum, R'de çok değişkenli sonuçların tekrarlanan ölçümlerin nasıl simüle edileceğidir? Örneğin, iki farklı tedavi grubu için 5 farklı zaman noktasında , ve tekrar tekrar …

9 r  repeated-measures  simulation 

Bunun muhtemelen temel bir soru olduğunu biliyorum ... Ama cevabı bulamıyorum. Ben bir Poisson ailesi ile bir GLM takıyorum ve daha sonra tahminlere bir göz atmaya çalıştım, ancak ofset dikkate alınmış gibi görünüyor: model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003), offset=(log(population)), data=data, subset=28:36, family=poisson()) predict (model_glm, type="response") Oranları değil davaları alıyorum ... Ben de denedim model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003)+ …

9 r  generalized-linear-model  prediction  offset 

Doğrusal bir SVM modeli için denklemi verebilecek bir pakete ihtiyacım var. Şu anda e1071'i şu şekilde kullanıyorum : library(e1071) m = svm(data, labels, type='C', kernel='linear', cost=cost, probability=FALSE, scale=scale) w = t(m$coefs) %*% data[m$index,] #Weight vector b = -model$rho #Offset Ancak, nasıl e1071::svm()pozitif ve negatif sınıfları seçtiğinden emin değilim , bu …

9 r  svm  e1071 

cor()Ve cor.test()işlevlerini kullanmada sorun yaşıyorum . Ben sadece iki matris (sadece sayısal değerler ve aynı satır ve sütun sayısı) var ve ben korelasyon numarası ve karşılık gelen p-değeri istiyorum. Kullandığımda cor(matrix1, matrix2)tüm hücreler için korelasyon katsayılarını elde ediyorum. Kor sonucu sadece tek bir sayı istiyorum. Buna ek olarak cor.test(matrix1, matrix2), …

9 r  correlation 

Bazı verileri sınıflandırmak için bir SVM oluşturmak için R kernlab paketini kullanıyorum . SVM, iyi bir doğruluk 'tahminleri' sağlaması açısından iyi çalışıyor, ancak giriş değişkenleri listem istediğimden daha büyük ve farklı değişkenlerin göreceli önemi konusunda emin değilim. En iyi eğitimli / en uygun SVM üreten girdi değişkenlerinin alt kümesini seçmek …

9 r  machine-learning  svm  genetic-algorithms 

Bunun gibi hayatta kalma verileriniz olduğunu varsayalım: obs <- data.frame( time = c(floor(runif(100) * 30), floor((runif(100)^2) * 30)), status = c(rbinom(100, 1, 0.2), rbinom(100, 1, 0.7)), group = gl(2,100) ) Standart bir log rank testi yapmak için survdiff(Surv(time, status) ~ group, data = obs, rho = 0) sağ? Peki ya …

9 r  survival  wilcoxon-signed-rank 

Bir yan hobi olarak, zaman serilerini tahmin ediyorum (özellikle R kullanarak). Verilerim için, günde yaklaşık 4 yıl öncesine kadar her gün ziyaretim var. Bu verilerde bazı farklı modeller vardır: Pazartesi-Cum çok ziyaret (Pazartesi / Sal en yüksek), ancak Sat-Sun önemli ölçüde daha az ziyaret var. Yılın belirli zamanlarında düşüş (yani …

9 r  time-series  seasonality  forecasting 

Kapalı. Bu soru konu dışı . Şu anda cevapları kabul etmiyor. Bu soruyu geliştirmek ister misiniz? Sorunuzu güncelleyin o yüzden -konu üzerinde Çapraz doğrulanmış için. 3 yıl önce kapalı . R kodumdan ve ggplot2 görselleştirmeden tekrarlayan raporlar oluşturmayı öğrenmekle çok ilgileniyorum. LaTeX'in olası bir cevap gibi göründüğünü ve en çok …

9 r 

R'de bir Parçacık Filtresi çalıştırmak için basit bir kod örneği arıyorum. Pompa paketi, durum uzay matematik bitini destekliyor gibi görünüyor, ancak örnekler, özellikle benim gibi basit bir OO geliştiricisi için programlı olarak takip etmek biraz zor gözlemlenen verilerin bir pompa nesnesine nasıl yükleneceği. Buradaki örnekler: http://cran.r-project.org/web/packages/pomp/vignettes/intro_to_pomp.pdf Diyelim ki girdi olarak …

9 r 

Kapalı. Bu soru konu dışı . Şu anda cevapları kabul etmiyor. Bu soruyu geliştirmek ister misiniz? Sorunuzu güncelleyin o yüzden -konu üzerinde Çapraz doğrulanmış için. 9 ay önce kapalı . Bu soruda belirtilen sorun R paket glmnet 1.7.3 sürümünde giderilmiştir. Ben aile = multinomial ile glmnet çalışan bazı sorunlar yaşıyorum, …

9 r  multinomial  glmnet 

Bu örnekte kullandığım termplotun kendi biraz geliştirilmiş versiyonunu oluşturdum, burada bulabilirsiniz . Daha önce SO üzerine yayınladım, ancak daha çok düşündüğümde, bunun muhtemelen gerçek kodlamadan daha çok Cox Orantılı tehlikeler modelinin yorumlanması ile ilgili olduğuna inanıyorum. Sorun Bir Tehlike Oranı grafiğine baktığımda, güven aralığının doğal olarak 0 olduğu bir referans …

9 r  survival  cox-model 

Ölçülen değişkenin ayrık pozitif tamsayıları (sayıları) olduğu bazı zaman serisi verileri var. Zaman içinde yukarı yönlü bir eğilim olup olmadığını test etmek istiyorum (ya da değil). Bağımsız değişken (x) 0-500 aralığında ve bağımlı değişken (y) 0-8 aralığındadır. Buna y = floor(a*x + b)sıradan en küçük kareler (OLS) kullanarak formun gerilemesini …

9 r  regression  python 

4 seviyeli bir değişkenim varsa, teoride 3 kukla değişken kullanmam gerekir. Uygulamada, bu aslında nasıl yapılır? 0-3 mü kullanıyorum, 1-3 mü kullanıyorum ve 4'ü boş mu bırakıyorum? Herhangi bir öneri? NOT: R'de çalışacağım. GÜNCELLEME: Yalnızca AD'ye karşılık gelen 1-4 kullanan bir sütun kullanırsam ne olur? Bu işe yarar mı yoksa …

9 r  regression  categorical-data  categorical-encoding 

4 kategorisi olan bazı tüketici verilerini görselleştirmeye çalışıyorum. Kullanıcılar farklı kategoriler arasında geçiş yapmakta serbesttir. Her bir birey için son üç veya dört anahtarı görselleştirmek istiyorum. Bu nedenle, 4 yığın oranına sahip bir sütun içeren bir arsa ile başlayacağız. Bundan sonra, her kategori insanların önceki vesileyle yaptıkları şeylere ayrılırken 16'mız …

9 r  data-visualization  sankey-diagram