Displaying 1 result from an estimated 1 matches for "fit_wind_speed".
2011 Apr 23
0
MASS fitdistr call in plyr help!
...locations. The data is
a data frame with two columns: site and wind speed. I want to split
the data on site and call a function to find the shape and scale
parameters of a weibull distribution fit.
The end result is a plot with x-axis = shape and y-axis = scale.
Currently my code looks like:
fit_wind_speed<-function(x){
x<-replace(x,x<=0,0.0001)
temp<-fitdistr(na.exclude(x[,1]),"weibull")
l<-length(names(x))
for(i in 1:l){
temp[i]<-(fitdistr(na.exclude(x[,i]),"weibull"))
}
temp
}
wind_speed_wide_dataframe<-function(x){
mini<-m...