search for: dpender

Hi, I am having some problems using the if statement correctly. I have used it many times previously so I dona't know what is different with this case. Here is my problem: I have a 1X10 matrix of values as follows: > H.MC [,1] [1,] 4.257669 [2,] 7.023242 [3,] 4.949857 [4,] 5.107000 [5,] 4.257669 [6,] 4.257669 [7,] 4.257669 [8,] 4.257669 [9,] 4.257669 [10,] 4.257669

R community, I am trying to create an array of the time differences between datapoints for a very large set. For some reason for 4 of the values the difference has been calculated as NA. Looking at the individual points two of them are "1981-03-29 01:40:00" and "1981-03-29 02:00:00" This is the exact same format as the other points that return values of 20 mins as expected.

Hi, For this example: O <- c(0 0 0 2 0 0 2 0) I want to create an array every time O[i] > 0. The array should be in the form; R[j] <- array(-1, dim=c(2,O[i])) i.e. if O[i] > 0 4 times I want 4 R arrays. Does anyone have any suggestions? Thanks, Doug -- View this message in context: http://r.789695.n4.nabble.com/Create-Arrays-tp2996706p2996706.html Sent from the R help

Hi, I am using the uniroot function in order to carry out a bivariate Monte Carlo simulation using the logistics model. I have defined the function as: BV.FV <- function(x,y,a,A) (((x^(-a^-1)+y^(-a^-1))^(a-1))*(y^(a-1/a))*(exp(-((1^(-a^-1)+y^(-a^-1))^a)+y^-1)))-A and the procedure is as follows: Randomly generate values of A~(0,1), y0 = -(lnA)^-1 Where: A=Pr{X<xi|Y=yi-1} and a is the

R community, I am trying to write a code that fills in data gaps in a time series. I have no R or statistics background at all but the use of R is proving to be a large portion of my PhD research. So far my code identifies where and the number of new entries required but I do not know how to add additional rows or columns into an array. Any advice on how this can be done? Here is an example:

I am trying to convert an array from numeric values back to date and time format. The code I have used is as follows; for (i in 0:(length(DateTime3)-1)) { DateTime3[i] <- (strptime(start, "%m/%d/%Y %H:%M")+ i*interval) where start <- [1] "1/1/1981 00:00" However the created array (DateTime3) contains [1,] 347156400 347157600 347158800 347160000 347161200 347162400

R community, I am creating a bivariate return level plot by adding calculated return period values as lines onto an existing plot using the following code with the points representing the return periods. plot(H2,D2,pch="+",axes=TRUE) points(H.10,D.10, type="l",col="blue") points(H.20,D.20, type="l",col="green") points(H.50,D.50,

...tical kites in KiteChart (plotrix) (Jim Lemon) 10. Re: vertical kites in KiteChart (plotrix) (Michael Bedward) 11. Re: Lattice: arbitrary abline in multiple histograms (Dennis Murphy) 12. Re: compare histograms (Dennis Murphy) 13. Re: Data Gaps (Dennis Murphy) 14. Date Time Objects (dpender) 15. Re: compare histograms (Michael Bedward) 16. Re: Lattice: arbitrary abline in multiple histograms (Alejo C.S.) 17. arima (nuncio m) 18. Re: vertical kites in KiteChart (plotrix) (elpape) 19. Re: LME with 2 factors with 3 levels each (Dennis Murphy) 20. Re: Data Gaps (dpender) 21....

...ltiple .csv-files and assigning them to variable names (Gabor Grothendieck) 32. Contract programming position at Merck (NJ, USA) (Liaw, Andy) 33. doubt in climate variability analysis in R! (govindas@msu.edu) 34. I got lot in the permutation (cbbaniya@gmail.com) 35. Re: Clustering (dpender) 36. Re: Please help me about Monte Carlo Permutation (Chitra) 37. Re: Cox Proportional Models and Haplotypes (Brad McNeney) 38. Re: draw path diagram using dot (eshi) 39. Re: Clustering (David Winsemius) 40. Tukey post hoc comparison (glht?) after 3factorial mixed model (lmer) (...

...cs i/o when using Rscript versus GUI (MacQueen, Don) 48. Re: complex transformation of data (Henrique Dallazuanna) 49. Re: complex transformation of data (Henrique Dallazuanna) 50. Information (Akash) 51. Re: complex transformation of data (Den) 52. Storm Clustering using clusters in evd (dpender) 53. confidence interval (Francesco Petrogalli) 54. ordering a vector (Francesco Petrogalli) 55. How to find data that includes certain values (poppinkid) 56. Re: User input in R program (jverzani) 57. Looping with incremented object name and increment function (Michael Costello) 58. Re:...

Hi, I am using the clusters function in the evd package in order to determine storm events from a wave time series. So far I have the code working as I want it for wave height on its own but I would now like to include the period as well. The input data is in the form of: H t H t H t so every height measurement has a corresponding period. The storms are defined when the wave height exceeds a

Hi, I am using the clusters function in the evd package and was wondering if anyone knew how to set the row names as date/time objects. My data has 1 column of date and time and another of wave height (H). My two options are: 1 - When in import the data using read.table make sure that the row.names is set = 1 giving row names to the H data. 2 - Import H from the data file and set the row names

Hi, I am trying to determine the MLE of the following function: http://r.789695.n4.nabble.com/file/n3317341/untitled.bmp I have defined both parts of the equation as separate functions and looped over the t and G values to get summations of each part. The lamda function has 3 unknowns which I am trying to determine using MLE bub tin order to try and get the overall function working these

Hi, I am looking for some help regarding the use of the mle function. I am trying to get mle for 3 parameters (theta0, theta1 and theta2) that have been defined in the the log-likelihood equation as theta0=theta[1], theta1=theta[2] and theta2=theta[3]. My R code for mle is: mle(Poisson.lik, start=list(theta=c(20,1,1), method="Nelder-Mead", fixed=list(w=w, t1=t1, t2=t2)) But I keep