Displaying 1 result from an estimated 1 matches for "wrongmeanshift".
2005 Nov 09
2
help with legacy R code
...2]
good <- scan("rawdata.txt", list(x=0))
xmin = 0
newx = good$x
xmean = mean(newx)
xmax = max(newx)+0.15
goodhist <- hist(newx, br=seq(from=0,to=xmax,by=0.15), probability=T,
col="lightyellow")
initmean <- (min(newx)+max(newx))/2
totalx <- length(newx)
wrongmeanshift <- 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 <- wrongs...