Displaying 1 result from an estimated 1 matches for "gplistrang".
Did you mean:
gplistrange
2011 Jul 16
2
ecdf() to nls() - how to transform data?
...ipedia.org/wiki/Gumbel_distribution
# ( - ( x - mue ) / beta )
# ( -e )^
# F(x) = e^
# formula for fitting-function
# formula: y ~ I( exp(1) ^ ( - exp(1) ) ^ ( - ( x - mue ) / beta ) )
# data: ecdf( hgd$V1 )
# start: list( mue=0.1, beta=0.025 )
gpfunction <- ecdf( hgd$V1 )
gplistrange <- seq( 0.0, 1, by=0.001 )
gplist <- gpfunction( gplistrange )
print( gplist )
print( class( gplist ) )
print("---")
#res <- nls( y ~ I( exp(1) ^ ( - exp(1) ) ^ ( - ( x - mue ) / beta ) ), df,
list( mue=0.1, beta=0.025 ), trace=TRUE )
#print( summary( res ) )
--
View this mess...