«r» etiketlenmiş sorular

İngiltere'de kafa travması veri tabanındaki retrospektif verilere dayanan bir lojistik regresyon modeli geliştiriyorum. Kilit sonuç 30 günlük mortalite ("Hayatta Kalma" ölçüsü olarak ifade edilir). Önceki çalışmalarda sonuç üzerinde önemli etkisi olduğuna dair kanıtları yayınlanmış diğer önlemler şunlardır: Year - Year of procedure = 1994-2013 Age - Age of patient = …

12 r  logistic  generalized-linear-model  goodness-of-fit  regression-strategies 

R version 3.1.1 (2014-07-10) -- "Sock it to Me" > bl <- c(140, 138, 150, 148, 135) > fu <- c(138, 136, 148, 146, 133) > t.test(fu, bl, alternative = "two.sided", paired = TRUE) Error in t.test.default(fu, bl, alternative = "two.sided", paired = TRUE) : data are essentially constant Sonra …

12 r  t-test 

Ben, R'ye daha aşina bir kullanıcıyım ve dört habitat değişkeni için 5 yaşın üzerindeki yaklaşık 35 kişi için rastgele eğimler (seçim katsayıları) tahmin etmeye çalışıyorum. Yanıt değişkeni, bir konumun "kullanıldı" (1) veya "kullanılabilir" (0) habitat (aşağıda "kullanım") olup olmadığıdır. Windows 64 bit bilgisayar kullanıyorum. R sürüm 3.1.0, aşağıdaki veri ve …

12 r  binomial  sas  random-effects-model  lme4-nlme 

Özet: Standart normal dağılım yerine lojistik regresyon katsayılarının testleri için dağılımının (rezidüel sapmaya dayalı serbestlik dereceleriyle) kullanımını destekleyen herhangi bir istatistiksel teori var mı ?ttt Bir süre önce SAS PROC GLIMMIX'a bir lojistik regresyon modeli takarken, varsayılan ayarların altında lojistik regresyon katsayılarının standart normal dağılım yerine dağılımı kullanılarak test edildiğini …

12 r  logistic  mathematical-statistics  sas  degrees-of-freedom 

Paketin skeçini kullanarak kafamı R'deki arabuluculuk paketinin etrafına getirmeye çalışıyorum. mediate()Fonksiyonun çıktısını anlamak için uğraşıyorum . require("mediation") require("sandwich") data("framing") med.fit <- lm(emo ~ treat + age + educ + gender + income, data = framing) out.fit <- glm(cong_mesg ~ emo + treat + age + educ + gender + income, …

12 r  mediation 

Ben bir laboratuvar için araştırma görevlisiyim (gönüllü). Ben ve küçük bir grup, büyük bir çalışmadan alınan bir veri kümesi için veri analizi ile görevlendirildik. Ne yazık ki, veriler bir tür çevrimiçi uygulama ile toplanmış ve verilerin en kullanışlı biçimde çıkması için programlanmamıştır. Aşağıdaki resimler temel sorunu göstermektedir. Bunun "Yeniden Şekillendirme" …

12 r  excel  data-cleaning 

Verilerim burada açıklanmıştır Tekrarlanan önlemler ANOVA takarken aov "Hata () modeli tekil hata" neden olabilir? lmerBenim temel durum kullanarak kullanarak bir etkileşimin etkisini görmeye çalışıyorum : my_null.model <- lmer(value ~ Condition+Scenario+ (1|Player)+(1|Trial), data = my, REML=FALSE) my.model <- lmer(value ~ Condition*Scenario+ (1|Player)+(1|Trial), data = my, REML=FALSE) Çalıştırmak anovabana önemli sonuçlar …

12 r  lme4-nlme 

Çok büyük bir veri setim var ve yaklaşık% 5 rasgele değerler eksik. Bu değişkenler birbiriyle ilişkilidir. Aşağıdaki örnek R veri kümesi sadece yapay korelasyonlu verilere sahip bir oyuncak örneğidir. set.seed(123) # matrix of X variable xmat <- matrix(sample(-1:1, 2000000, replace = TRUE), ncol = 10000) colnames(xmat) <- paste ("M", 1:10000, …

12 r  random-forest  missing-data  data-imputation  multiple-imputation  large-data  definition  moving-window  self-study  categorical-data  econometrics  standard-error  regression-coefficients  normal-distribution  pdf  lognormal  regression  python  scikit-learn  interpolation  r  self-study  poisson-distribution  chi-squared  matlab  matrix  r  modeling  multinomial  mlogit  choice  monte-carlo  indicator-function  r  aic  garch  likelihood  r  regression  repeated-measures  simulation  multilevel-analysis  chi-squared  expected-value  multinomial  yates-correction  classification  regression  self-study  repeated-measures  references  residuals  confidence-interval  bootstrap  normality-assumption  resampling  entropy  cauchy  clustering  k-means  r  clustering  categorical-data  continuous-data  r  hypothesis-testing  nonparametric  probability  bayesian  pdf  distributions  exponential  repeated-measures  random-effects-model  non-independent  regression  error  regression-to-the-mean  correlation  group-differences  post-hoc  neural-networks  r  time-series  t-test  p-value  normalization  probability  moments  mgf  time-series  model  seasonality  r  anova  generalized-linear-model  proportion  percentage  nonparametric  ranks  weighted-regression  variogram  classification  neural-networks  fuzzy  variance  dimensionality-reduction  confidence-interval  proportion  z-test  r  self-study  pdf 

step()R'deki fonksiyon tarafından kademeli regresyon için kullanılan kritik p-değeri nedir? Sanırım 0.15, ama varsayım doğru mu? Kritik p değerini nasıl değiştirebilirim?

12 r  regression  p-value  stepwise-regression 

Doğrusal bir modelimiz olduğunu Model1ve vcov(Model1)aşağıdaki matrisi verdiğimizi varsayalım : (Intercept) latitude sea.distance altitude (Intercept) 28.898100 -23.6439000 -34.1523000 0.50790600 latitude -23.643900 19.7032500 28.4602500 -0.42471450 sea.distance -34.152300 28.4602500 42.4714500 -0.62612550 altitude 0.507906 -0.4247145 -0.6261255 0.00928242 Bu örnek için, bu matris gerçekte ne gösteriyor? Modelimiz ve bağımsız değişkenleri için güvenle hangi varsayımlarda …

12 r  interpretation  multicollinearity  assumptions  covariance-matrix 

Ben tarafından optimize edilmiş işlevi için genelleştirilmiş doğrusal olmayan en küçük kareler regresyonu için günlük olasılığını hesaplamaya çalışıyorum R paketin fonksiyon , (Brown hareketi varsayarak aa filogenetik ağaç mesafe tarafından meydana varyans kovaryans matrisi kullanılarak gelen paket). Aşağıdaki yeniden üretilebilir R kodu, gnls modeline x, y verileri ve 9 taksonlu …

12 r  maximum-likelihood  least-squares  nonlinear-regression  mixed-model 

Soru Üç grubun test puanları R'de ayrı vektörler olarak kaydedilir. set.seed(1) group1 <- rnorm(100, mean = 75, sd = 10) group2 <- rnorm(100, mean = 85, sd = 10) group3 <- rnorm(100, mean = 95, sd = 10) Bu gruplar arasındaki medyanlarda önemli bir fark olup olmadığını bilmek istiyorum . …

12 r  hypothesis-testing  multiple-comparisons  mean  median 

1000 belgeyi ve içinde görünen tüm kelimeleri temsil eden bir veri setim var. Böylece satırlar belgeleri ve sütunlar kelimeleri temsil eder. Yani, örneğin, hücrede değer kez kelimenin açılımı belgede oluşur . Şimdi, tf / idf yöntemini kullanarak kelimelerin 'ağırlıklarını' bulmak zorundayım, ama aslında bunu nasıl yapacağımı bilmiyorum. Birisi bana yardım …

12 r  data-mining  feature-selection 

Ben işlevi bir izlenimi lmer()de lme4paketin p-değerleri (bkz vermedi lmer, p-değerleri ve tüm bu ). MCMC oluşturulan p değerleri yerine şu soruya göre kullanıyorum: Karışık model ve bu soruda önemli etkisilme4 : içinde paketinde gelen p-değerleri bulunamıyorlmer()lm4R . Son zamanlarda modelimingetSummary.mer() sabit efektlerini bir csv dosyasına almak için memisc ve …

12 r  mixed-model  p-value  mcmc  lme4-nlme 

Bazı kod ve çıktı var ve bir model oluşturmak istiyorum. Bu çıktıyı kullanarak bir model oluşturmak nasıl bilmiyorum: require("splines") x <- c(0.2, 0.23, 0.26, 0.29, 0.33, 0.46, 0.53 ) y <- c(0.211, 0.2026, 0.2034, 0.2167, 0.2177, 0.19225, 0.182) fit <- lm(y ~ ns(x,3)) summary(fit) Not ns()doğal kübik eğri B-parça baz …

12 r  splines