Akarların fonksiyonel tepkisi alanında araştırma yapıyorum. Rogers tip II fonksiyonunun parametrelerini (saldırı hızı ve işleme süresi) tahmin etmek için bir gerileme yapmak istiyorum. Bir ölçüm veri setim var. Aykırı değerleri en iyi nasıl belirleyebilirim?
Regresyonum için R'de (doğrusal olmayan bir regresyon) aşağıdaki komut dosyasını kullanıyorum: (tarih seti, değerler (ilk av sayısı) ve değerler (24 saat boyunca yenen av sayısı ) data.txtile dosya adı verilen basit bir 2 sütunlu metin dosyasıdır :N0FR
library("nlstools")
dat <- read.delim("C:/data.txt")
#Rogers type II model
a <- c(0,50)
b <- c(0,40)
plot(FR~N0,main="Rogers II normaal",xlim=a,ylim=b,xlab="N0",ylab="FR")
rogers.predII <- function(N0,a,h,T) {N0 - lambertW(a*h*N0*exp(-a*(T-h*N0)))/(a*h)}
params1 <- list(attackR3_N=0.04,Th3_N=1.46)
RogersII_N <- nls(FR~rogers.predII(N0,attackR3_N,Th3_N,T=24),start=params1,data=dat,control=list(maxiter= 10000))
hatRIIN <- predict(RogersII_N)
lines(spline(N0,hatRIIN))
summary(RogersII_N)$parameters
Calssic Artık grafiklerini çizmek için aşağıdaki komut dosyasını kullanıyorum:
res <- nlsResiduals (RogersII_N)
plot (res, type = 0)
hist (res$resi1,main="histogram residuals")
qqnorm (res$resi1,main="QQ residuals")
hist (res$resi2,main="histogram normalised residuals")
qqnorm (res$resi2,main="QQ normalised residuals")
par(mfrow=c(1,1))
boxplot (res$resi1,main="boxplot residuals")
boxplot (res$resi2,main="boxplot normalised residuals")
Sorular
- Hangi veri noktalarının aykırı olduğunu en iyi nasıl belirleyebilirim?
- R'de kullanabileceğim objektif olan ve hangi veri noktalarının aykırı olduğunu gösteren testler var mı?