আমি মাইটগুলির কার্যকরী প্রতিক্রিয়ার ক্ষেত্রে গবেষণা করছি। আমি রজার্স টাইপ II ফাংশনের প্যারামিটারগুলি (আক্রমণ হার এবং হ্যান্ডলিংয়ের সময়) অনুমান করতে একটি রিগ্রেশন করতে চাই। আমার কাছে পরিমাপের একটি ডেটাসেট রয়েছে। আমি কীভাবে সর্বোপরি বিদেশী নির্ধারণ করতে পারি?
(Dateset একটি সহজ 2 কলামটি পাঠ্য নামক ফাইল: আমার রিগ্রেশন জন্য আমি আর নিম্নলিখিত স্ক্রিপ্ট (ক অ রৈখিক রিগ্রেশনের) ব্যবহার data.txtসাহায্যে ফাইল N0মান (প্রাথমিক শিকারী নম্বর) FRমান (24 ঘন্টা সময় খাওয়া শিকার সংখ্যা):
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ক্যালাসিক অবশিষ্টাংশের গ্রাফগুলি তৈরি করার জন্য আমি নিম্নলিখিত স্ক্রিপ্টটি ব্যবহার করি:
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")প্রশ্নাবলি
- কোন ডেটা পয়েন্টগুলি বিদেশী বলে আমি কীভাবে সর্বোত্তমভাবে নির্ধারণ করতে পারি?
- আমি কি আর তে পরীক্ষামূলক ব্যবহার করতে পারি যা উদ্দেশ্যমূলক এবং আমাকে দেখায় যে কোন ডেটা পয়েন্টগুলি বিদেশী?