Displaying 1 result from an estimated 1 matches for "fitrat".
Did you mean:
firat
2005 Nov 09
2
help with legacy R code
...wrongmean + 0.2
wrongper <- pgamma(wrongmeanshift, wrongshape, wrongrate)
nfalseundermean <-
which(abs(newx-wrongmeanshift)==min(abs(newx-wrongmeanshift)))
initnfalse <- nfalseundermean / wrongper
fitmean <- -1
fitsd <- 0
fitnfalse <- initnfalse
fitshape <- wrongshape
fitrate <- wrongrate
curve((fitnfalse*dgamma(x,fitshape,fitrate))/totalx, add=T, col="red",
lwd=2)
breaksllength <- length(goodhist$breaks)
endi = breaksllength - 1
binprob = c(1)
for (i in 1:endi) {
expnegative <- fitnfalse * (pgamma(goodhist$breaks[i+1],wrongshape,
wrongrate)-...