similar to: ML fit of gamma distribution to grouped data

Good day everyone. I have a large dataset of 1 min wind speeds covering 5 years. How can I make an array of maximum daily values? The vectors I have are: 'VDATE' with dates in format '%Y-%m-%d' (like '1992-10-28') and 'WS' with wind speed data (same number of elements as VDATE). I want an array with 2 columns: Max daily wind speed and corresponding day. Has

Good day everyone, I want to use zoo(aggregate) to calculate 15-min mean values from a wind dataset which has 1-min values. The data I have looks like this: vector VDATE vector WS 1 1998-10-22:02:11 12.5 2 1998-10-22:02:12 10.1 3 1998-10-22:02:13 11.2 4 1998-10-22:02:14 10.5 5 1998-10-22:02:15 11.5 . . . n

Good day everyone, Can anyone suggest an effective method to solve the following problem: I have 2 dataframes D1 and D2 as follows: D1: dates ws wc pwc 2005-10-19:12:00 10.8 80 81 2005-10-20:12:00 12.3 5 15 2005-10-21:15:00 12.3 3 15 2005-10-22:15:00 11.3 13 95 2005-10-23:12:00 12.3 13 2 2005-10-24:15:00 10.3 2 95 2005-10-25:15:00 10.3 2 2 D2:

Good day, Is there any way to save a plot produced by R in a LINUX (Debian) machine? The window opened by R to put the plot in, does not give any option to save it (there are options to move, close, minimise it, etc. but not to save it). How do you do that? Thanks, Augusto -------------------------------------------- Augusto Sanabria. MSc, PhD. Mathematical Modeller Risk Research Group

Good day everyone, I want to assess the error when fitting a Gram-Charlier CDF to some data 'ws', that is, I want to calculate: Err = |ecdf(ws) - GCh_ser(ws)| The problem is, I cannot get the F(x) values from the ecdf. 'Summary(ecdf())' returns some of the x-axis values, but how do you get the F(x) values? Thank you for any help you can provide. Regards, Augusto

Good day everyone, I am new in R programming (my question may sound trivial to you): is there any way to ask the user to enter a string within an R process, say a filename, make R to recognise it and open the given file? It is a simple exercise in other languages. I am using R2.1.1 in a LINUX machine. Thanks for your help. Augusto -------------------------------------------- Augusto

Good day everyone, I want to use the Gram-Charlier series expansion to model some data. To do that, I need functions to: 1) Calculate 'n' moments from given data 2) Transform 'n' moments to 'n' central moments, or 3) Transform 'n' moments to 'n' cumulants 4) Calculate a number of Hermite polynomials Are there R-functions to do any of the above?

Does anyone know of a resource for learning the basics of how to manage and manipulate dates and times in R? I have been reading "Introductory Statistics with R" by Peter Dalgaard which is fantastic. But alas, I could find no reference to date and time. I have looked at the reference manual but it is particularly unapproachable. So rather than dense technical talk I would rather see a

Dear list members, I have a set of claims data, which are in ranges and the shape of the distribution is relatively different. I have looked through R help threads and found out that an ideal way is suggested for the gamma distribution ML fitting for grouped data. I just wonder if there is any method that works for lognormal or pareto distribution? An example would be:

Hi all, I would like to fit a gamma pdf to my data using the method of RSS (Residual Sum-of-Squares). Here are the data: x <- c(86, 90, 94, 98, 102, 106, 110, 114, 118, 122, 126, 130, 134, 138, 142, 146, 150, 154, 158, 162, 166, 170, 174) y <- c(2, 5, 10, 17, 26, 60, 94, 128, 137, 128, 77, 68, 65, 60, 51, 26, 17, 9, 5, 2, 3, 7, 3) I have typed the following code, using nls method:

Dear R-lister Could anyone know how to make a multi-line title for a plot with expression()? In my plot, the title should be writeen in two lines (because it is two long for one line) and it should use a mathematical expression. I tried to use "\n", but "\n" is ignored in expression() call: hist(diffChangeRequestHintsBeforeAnswering[,4], br = 50, xlab = "Skill Change in

Hello everyone, With th bar histogram (number of occurrences) hist<-c(24,7,4,1,2,1,1) of seven equally spaces classes ]1-4], ]5-8], ]9-12], ]13-16], ]17-20], ]21-24], ]25-28], I obtained shape=0.8276 and rate=0.1448. I would like to know how to build the continuous pdf of a this gamma distribution knowing these two estimated parameters such that I will be able to predict the pdf of any

Hi, I am new to R and I'm looking to perform some method comparison analysis and would be grateful for some help regarding package installation. There is a package 'MethComp' that I want to use (http://staff.pubhealth.ku.dk/~bxc/MethComp/Courses/Sthm.2007/.) Unfortunately I am unable to install it. I have tried numerous approaches but to no avail. this is the first time I have tried to

Hi, all How can generate a sample from truncated inverse gamma distribution in R? thanks

Colleagues, I am using R 2.4.0 on both a Mac (10.4.8) and Linux (RedHat 9). To read data from an Excel spreadsheet, I do "save as" in Excel, then select the "Text (tab-delimited)" format. The resulting file uses a tab separator and I can usually read the file using read.delim. Sometimes, the header row contains embedded carriage returns. When I view the file,

Hi all, I have a very simple question that I can't seem to find the answer to. How do I extract rows that meet a certain criteria from a data frame and group them into a new data frame? For example, if I want to make a new data frame that only includes rows of data for which the p values (given by one of the columns in the data frame) are less than a certain value, how do I do this? It seems

Dear All: I create a file named :"P_Value" with only one simple function: P_Value <- function( Table ) { S = fisher.test(Table, alternative = "two.sided"); return(S$p.value); } However, it seems that it's impossible to use this function directly, because no matter where I save this file, R always reports

Dear Sir/Madam, I apologize for any cross-posting. I got a simple question, which I thought the R list may help me to find an answer. Suppose we have Y_1, Y_2, ., Y_n ~ Poisson (Lambda_i) and Lambda_i ~Gamma(alpha_i, beta_i). Empirical Bayes Estimator for hyper-parameters of the gamma distr, i.e. (alpha_t, beta_t) are needed. y=c(12,5,17,14) n=4 What about a Hierarchal B ayes

Hi R users, Im trying to install XML source package (in window platform) using the following command. *install.packages(pkgs = "C:\\Documents and Settings\\skang\\Desktop\\Softwares\\R\\Library\\XML_2.6-0.tar.gz", repos = NULL, type = "source")* and get the following result. gzip: stdin: unexpected end of file /cygdrive/c/Rtools/bin/tar: Unexpected EOF in archive

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20071006/065906cc/attachment.pl