[R] Error message in `fitdistr` function in the MASS package

2014-02-26 Thread kmmoon100
I have one issue at the moment while I am using a function, fitdistr in package MASS. I have 12 years of wind data with 30 min average interval. My data look like below: Timewindspeed 1/01/2005 0:00 15 1/01/2005 0:30 11 1/01/2005 1:00 11 1/01/2005 1:30 13 1/01/2005 2:00 15

[R] How can I apply “Sapply” in R with multiple codes in one function?

2013-12-09 Thread kmmoon100
Hello everybody, I am a new R user..(still;;;). I have a simple "sapply" function example for calculating mean and sd of splited data scale. My data contains half hourly wind speed with direction. I want to know daily weibull distribution for my study for 13 years. That is why my dataset is split

[R] How can I get seasonal variation table from generalized additive mixed models in R?

2013-11-18 Thread kmmoon100
Hello everyone, I used a function called gamm in mgcv in R program to investigate seasonal variation of wind speed based on 13 years of measurement dataset. I would like to use this variation in my formula for my study but as I am new in R, it's extremely hard to find values on y-axis.. too many s

[R] Plotting multiple weibull distributions in one graph

2013-11-12 Thread kmmoon100
Hello, I am trying to plot multiple weibull distributions in one graph. Does anyone know how to produce? I already have different parameters for weibull distributions (shape and scale) Here are examples. Shape Scale 20112.455 21.657 20102.328 21.486 20092.336 22.642 2

Re: [R] Comparison of two weibull distributions

2013-11-04 Thread kmmoon100
Thank you so much for your explanation. I might be in a trouble again with processing log-likelihood analysis. If it happens, may I ask your instuctions next time? -- View this message in context: http://r.789695.n4.nabble.com/Comparison-of-two-weibull-distributions-tp4679632p4679680.html Sen

[R] Formatting time stamps for weibull distribution

2013-11-04 Thread kmmoon100
Hello everyone, I have been using bReeze package for weibull distribution analysis. In this package, there is a function called 'formatTS' for converting time stamps from string to POSIXlt. (http://rgm3.lab.nig.ac.jp/RGM/R_rdfile?f=bReeze/man/formatTS.Rd&d=R_CC) This is a part of process to plot a

[R] Comparison of two weibull distributions

2013-11-03 Thread kmmoon100
Hello, How can I do a test of two weibull distributions? I have two weibull distribution sets from two wind datasets in order to check whether they are same. I thought 2 sample t-test would be applicable but I couldn't find any ways to do that on the Internet. Does anyone know what type of test is

Re: [R] Optimization failed in fitdistr (Weibull distribution)

2013-10-28 Thread kmmoon100
Hi Berend, Thank you for your reply. How can I use dput function for this type of data? I looked up the description of the function but I still can't understand how to use it for solving my error. Regards, Kangmin. -- View this message in context: http://r.789695.n4.nabble.com/Optimization-f

[R] "Optimization fail" error from fitdistr (Weibull distribution)

2013-10-28 Thread kmmoon100
Hello everyone, This is Kangmin. I am trying to produce shape and scale of my wind data. My data is based on wind speed frequency with 1km/hr increment. data is described below. Windspeed (km/h)Frequency 1 351 2 147 3 317 4 378 5 527 6 667 7 865 8 970 9 987 10 907 11 905 12 64

[R] Optimization failed in fitdistr (Weibull distribution)

2013-10-28 Thread kmmoon100
Hello everyone, This is Kangmin. I am trying to produce shape and scale of my wind data. My data is based on wind speed frequency with 1km/hr increment. data is described below. Windspeed (km/h)Frequency 1 351 2 147 3 317 4 378 5 527 6 667 7 865 8