2016-04-08 30 views
-1

Çok sayıda verinin üretildiği bir analitik laboratuarında çalışıyorum. Model uyumu yapıyoruz ve model katsayıları (intercept, a, b), r_squared ve rezidüel standart hata ile ilgileniyoruz. 70'e kadar bağımlı değişken (yanıt) sayısı oldukça büyüktür. Mylms'den satırların a ve f değerleri ile belirtilen model parametrelerinin olduğu bir tablo oluşturmak istiyorum. Birçok doğrusal modelden özelleştirilmiş sonuç tablosu oluşturma

benim davaya net bulunan çözümler adapte değildi:

çözümün bir parçası elbette özetidir (mylms) nerede parametreler saklanır. Bunu yapabilecek

#To generate a dataset 
x <- c(1:30) 
x2 <- x*x 
a <- x2 
b <- jitter(a, factor=10) 
c <- jitter(b, factor=30) 
d <- jitter(c, factor=40) 
e <- jitter(d, factor=50) 
f <- jitter(e, factor=60) 
z <- cbind(a,b,c,d,e,f) 

# The matrix y are the responses of x 
y <- as.matrix(z) 

# Define the weights for the regression 
w <- 1/x 

# create linear models 
mylms <- lm(y~x+x2, weights=w) 
+0

Böylelikle, alanların katsayılar olduğu en fazla 70 sütununuz var mı? Ve kaç tane satır var? Her zaman 6 (A-F)? – Simon

+1

Süpürge paketine göz atın: https://cran.r-project.org/web/packages/broom/index.html –

+0

Aslında, yanıt matrisinde y sütun sayısı değişecektir. 30 olabilir, ancak 70 veya 15 olabilir. Matris y'deki sütun sayısı, sonuç tablomdaki satırların sayısına eşittir. Ancak sonuç tablosunda her zaman 5 satır olacaktır: kesişme, a, b, r_squared ve artık standart hata. – ThomasG

cevap

0

:

mylms$residuals 

veya örneğin artıkların aralığını elde etmek için:

İşte

ben bir çözüm bulmak istiyoruz basitleştirilmiş bir örnektir

veya miktarları almak için örneğin:

apply(mylms$residuals,2,quantile,seq(0.1,0.9,by=0.1)) 
+0

Aslında rezidüel standart hatası arıyoruz, her bir artık kendisi değil! BOŞ R2, <- - – ThomasG

+0

I çözeltisi Rest.stabw ThomasG

İlgili konular