[R] sum to infinity

2009-03-25 Thread Roslina Zakaria
Hi r-users, How do we evaluate the summation of (1/m!) from 0 to infinity (for example). Any help is very much appreciated. Thank you. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] newton method

2009-03-22 Thread Roslina Zakaria
Hi R-users, Does R has a topic on newton's method? Thank you for the info. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

[R] extract positive pairs

2009-01-27 Thread Roslina Zakaria
Hi, I have a data below and would like to search for positive pairs only and form a new data set.    X1             X2 31.0 9.0 11.0 1.0 1.0 0.0 0.0 0.0 8.0 0.0 0.0 0.0 2.0 2.0 18.0 3.0 0.0 0.0 0.0 0.0 0.0 0.0 10.0 0.0 6.0 0.0 ... The new data will be X1'   X2' 31.0 9.0

[R] recursive term

2008-12-09 Thread Roslina Zakaria
Hi,   I would like to write a function for this formula: F(a,b,c,z)=1+(ab)/(1!c)+ +(a(a+1)b(b+1))/(2!c(c+1))+ +(a(a+1)(a+2)b(b+1)(b+2))/(3!c(c+1)(c+2))+…   I wrote this function but not sure what is not right:   hypergeo_sum - function (a,b,c,z,n) { for (i in 1:n)   { aa - 1+(a*b*z)/c   aa[i] -

[R] hypergeometric function

2008-12-03 Thread Roslina Zakaria
Hi,   I hope somebody can help me on how to use the hypergeometric function.  I did read through the R documentation on hypergeometric but not really sure what it means.   I would like to evaluate the hypergeometric function as follows: F((2*alpha+1)/2, (2*alpha+2)/2 , alpha+1/2, betasq/etasq).  

[R] histogram

2008-05-16 Thread Roslina Zakaria
Hi r-expert, I would like to plot histogram using frequency not density. But I got the following warning.  obs.hist - hist(jan_data2[,4],right=FALSE,breaks=c(0,5,10,15,20,100),freq=TRUE, + xlab=Rain amt (mm),ylim=c(0,3000), + main=Frequency of observed, Jan (1901-1990),

[R] How to swap and rearrange rows?

2008-05-16 Thread Roslina Zakaria
Dear R-expert, How to swap and rearrange the row so that I will have Jan-Dec in order? est31     p0 est.alpha est.beta  est.rate Jan  0.8802867 0.7321440 7.241757 0.1380880 Mar  0.8598566 0.7096567 7.376367 0.1355681 May  0.6204301 0.8657272 6.036106 0.1656697 July 0.5032258 0.9928488

[R] How to paste graph from R in Latex?

2008-05-16 Thread Roslina Zakaria
Dear R-expert, Is it possible to save graph from R into Latex document? I can see save as metafile , PNG, pdf etc, but I'm not sure which one to use. Thank  you so much for your help. __ R-help@r-project.org mailing list

[R] plot 2 bars

2008-05-13 Thread Roslina Zakaria
 Hi r-expert, How do I plot this 2 histogram side by side so that it is easy to compare for bin range 0,5,10,15, 20 and more     obs_data  pre_gam [1,] 2695 2677.284 [2,]   43   61.101 [3,]   20   29.016 [4,]   14   11.160 [5,]   18   11.439 Thanks in advance.

[R] How to eliminate perticular date

2008-05-11 Thread Roslina Zakaria
Hi R-expert, I try to eliminate 29 Feb but I got an error message below: feb_data1 - Pooraka_data[Pooraka_data$Month==2,] feb_28days - feb_data1 [feb_data1$Day==28,] feb_29days - feb_data1 [feb_data1$Day==29,] ## delete 29 Feb feb_no_29 - feb_data1 [-(feb_29days),] feb_no_29 - feb_data1

[R] histogram

2008-05-08 Thread Roslina Zakaria
Dear R-expert, For histogram function, can we get the table of bin and frequency like in excel, together with the histogram? Therefore, we can check the number of data included. Thank you so much for your attention and help.

[R] convert S-Plus function to R function

2008-05-06 Thread Roslina Zakaria
Hi R-expert, If I have this daily rainfall data, how do call a particular set of data? Year,Month,Day,Amount 1981,1,1,0 1981,1,2,0 1981,1,3,78 1981,1,4,22 1981,1,5,2 1981,1,6,0 1981,1,7,0 1981,1,8,0 1981,1,9,0 1981,1,10,10 1981,1,11,0 1981,1,12,108 1981,1,13,328 1981,1,14,10 1981,1,15,0

[R] to extract particular date/data

2008-05-01 Thread Roslina Zakaria
Hi R-expert, If I have this daily rainfall data, how do call a particular day? Year,Month,Day,Amount 1900,12,22,1.3 1900,12,23,0 1900,12,24,0 1900,12,25,0 1900,12,26,0 1900,12,27,0 1900,12,28,0 1900,12,29,4.8 1900,12,30,0.3 1900,12,31,0.5 1901,1,1,0 1901,1,2,3 1901,1,3,0 1901,1,4,0.5 1901,1,5,0

[R] looping problem

2008-04-24 Thread Roslina Zakaria
Dear r-expert, I would like to generate 30 sets of random numbers from uniform distribution (0,1). Each set of random numbers should have 90 data.  Here is my code: rand.no - function(n,itr) { for (i in 1:itr)   {rand.1 - runif(n,0,1)   if (i ==1) rand.2 - rand.1   else rand.2 -

[R] looping problem

2008-04-13 Thread Roslina Zakaria
Hi R-users, I would like to do looping for this process below to estimate alpha beta from gamma distribution: Here are my data: day_data1 - 123456 789 10 11 12 13 14 15 16 17 18 19 20 21 22 23 1943 48.3 18.5 0.0 0.0 18.3 0.0

[R] chi-square test

2008-04-08 Thread Roslina Zakaria
Hi R-users, I would like to find the goodness of fit using Chi-suare test for my data below: xobs=observed data, xtwe=predicted data using tweedie, xgam=predicted data using gamma xobs - c(223,46,12,5,7,17) xtwe - c(217.33,39,14,18.33,6.67,14.67) xgam - c(224.67,37.33,12.33,15.33,5.33,15)

[R] tweedie model

2008-04-07 Thread Roslina Zakaria
Hi R-users, I have this code below for tweedie inverse, but I don’t understand why it doesn’t give zero for some of r1 values?   r1 - runif(31,0,1) p0 - 0.72 ; p - 1.63; mu - 1.48; phi - 9.61   tweedie.inv - function(rand,p0,p,mu,phi) {  y - rand    ind - rand p0    y[ind] - 0    y[!ind] -

[R] solver in R

2008-03-31 Thread Roslina Zakaria
Dear R-users, Is there any function in R that works similarly like solver in Excel. I have a set of daily rainfall data and I would like to estimate alpha and beta for the gamma function. Here is my daily rainfall data: [1] 0.2 1.2 0.0 0.0 0.0 0.0 0.0 0.6 0.0 0.0 0.0 0.0 0.0 0.0

[R] mean for each row

2008-03-16 Thread Roslina Zakaria
Hi r-users, How do find the mean for each row? Thank you in advance for your help. 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 Day Totals 10.0 0.0 0.0 0.0 0.6 0.0 8.4 0.0 29.4 0.0 38.4 20.0 0.0 1.8 0.0 22.4 0.0 0.2 0.4 0.8 0.0

[R] tweedie model

2008-03-05 Thread Roslina Zakaria
Hi R users, Anybody using tweedie model to analyze rainfall data? Thanks in advance for your attention. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] quantile-residual plot

2008-02-28 Thread Roslina Zakaria
Hi, What is quantile residuals? Thank you so much in advance. Be a better friend, newshound, and __ R-help@r-project.org mailing list

<    1   2