[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 http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[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 commented, minimal, self-contained, reproducible code.


[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 
11.0 1.0 
2.0 2.0 
18.0 3.0 

I tried to write the function as:

y1y2 -  read.csv(genX1X2.csv, header=FALSE)
(y1y2[,1]  0 )    (y1y2[,2]0)
cbind(y1y2[,1],y1y2[,2])
 
Thank you for your help.




__
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 commented, minimal, self-contained, reproducible code.


[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] - (aa-1)*(a+i)*(b+i)*z/((i+1)*(c+i))
  comb_sum - aa + aa[i]+ aa[i+2]
  }
  comb_sum
}
 
hypergeo_sum (1.25,1.75,1.25,0.5,3)
 
output:
 hypergeo_sum (1.25,1.75,1.25,0.5,3)
[1] 2.394531   NA 1.039062

 
The answer should be 2.852539.
 
Or can you suggest any R book that I can refer to. Thank you so much for your 
help.




__
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 commented, minimal, self-contained, reproducible code.


[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).
 
where 
alpha - .75; beta1 - 7 ; beta2 - 5.5;
etasq - ((beta1+beta2)/(2*beta1*beta2*(1-rho))) ^2
betasq - ((beta1-beta2)^2+4*beta1*beta2*rho)/(4*beta1^2*beta2^2*(1-rho)^2) 
 
I’m not sure which function should be used- either phyper or  qhyper or dhyper
 
Thank you so much for your help.


   
__
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 commented, minimal, self-contained, reproducible code.


[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), Pooraka)
Warning message:
In plot.histogram(r, freq = freq1, col = col, border = border, angle = angle,  :
  the AREAS in the plot are wrong -- rather use freq=FALSE

I don't really understand the warning message and hope anybody can help me.
Thanks in advance for your help.  





__
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 commented, minimal, self-contained, reproducible code.


[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 4.027408 0.2482986
Aug  0.5322581 0.9625738 4.103121 0.2437169
Oct  0.6792115 0.8526226 5.105218 0.1958780
Dec  0.8397849 0.7490287 7.070349 0.1414357
 est30
    p0 est.alpha est.beta  est.rate
Apr  0.7296296 0.7929348 6.303877 0.1586325
Jun  0.5574074 0.8588608 5.695905 0.1755647
Sept 0.607 0.9031150 4.594891 0.2176330
Nov  0.7725926 0.7600906 5.636366 0.1774193
 est28
  p0 est.alpha est.beta  est.rate
Feb 0.877262 0.6567584 8.708051 0.1148363
Thank you so much.




__
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 commented, minimal, self-contained, reproducible code.


[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
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[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-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 commented, minimal, self-contained, reproducible code.


[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 [-(feb_29days),]
Error in xj[i] : invalid subscript type 'list'
Thank you so much for your attention.


  


[[elided Yahoo spam]]

__
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 commented, minimal, self-contained, reproducible code.


[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.


  


[[elided Yahoo spam]]

__
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 commented, minimal, self-contained, reproducible code.


[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
1981,1,16,6
1981,1,17,0
1981,1,18,28
1981,1,19,0
1981,1,20,0
1981,1,21,0
1981,1,22,0
1981,1,23,0
1981,1,24,0
1981,1,25,32
1981,1,26,0
1981,1,27,2
1981,1,28,0
1981,1,29,2
1981,1,30,0
1981,1,31,0
1981,2,1,0
1981,2,2,118
1981,2,3,0
1981,2,4,0
1981,2,5,0
1981,2,6,12
1981,2,7,0
1981,2,8,0
1981,2,9,46
1981,2,10,4
1981,2,11,0
1981,2,12,0
1981,2,13,0
1981,2,14,0
1981,2,15,0
1981,2,16,0
1981,2,17,0
1981,2,18,0
1981,2,19,0
1981,2,20,0
1981,2,21,4
1981,2,22,16
1981,2,23,0
1981,2,24,0
1981,2,25,2
1981,2,26,0
1981,2,27,0
1981,2,28,0
1981,3,1,0
1981,3,2,0
1981,3,3,0
1981,3,4,0
1981,3,5,0
1981,3,6,0
1981,3,7,0
1981,3,8,0
1981,3,9,106
1981,3,10,10
1981,3,11,0
1981,3,12,0
1981,3,13,0
1981,3,14,0
1981,3,15,32
...
For example, if I would like to extract january and february data, and use 
1/1/1901 as the origin. How do I use the dates function in R ? I did try to 
lookup in the manual, but not sure how to apply it?
This is my code in S-plus.
define.date1-function(dt1,mt1,mt2,nn,da)
{  mt2-mt2+1
   start-julian(mt1, 1, nn, origin=c(month=1, day=1, year=1971))+1
   end-julian(mt2, 1, nn, origin=c(month=1, day=1, year=1971))+da
   a-dt1[start:end,]
   am-as.matrix(a[,5])
}
Thank you so much for your kind attention and help.


  

Be a better friend, newshound, and 
know-it-all with Yahoo! Mobile.  Try it now.  
http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ

__
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 commented, minimal, self-contained, reproducible code.


[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
1901,1,6,0
...
I used to use julian.date in S-Plus.
Thank you so much for your kind attention and help.


  


[[elided Yahoo spam]]

__
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 commented, minimal, self-contained, reproducible code.


[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 - cbind(rand.2,rand.1)
  }
  rand.2
}
Question: The code is okay is just that it give me 31 sets instead of 30 sets. 
Can anybody any adjustment to my code.
Your help is really appreciated.
Thank you.


  

Be a better friend, newshound, and 
know-it-all with Yahoo! Mobile.  Try it now.  
http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ

__
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 commented, minimal, self-contained, reproducible code.


[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   0.0  0.0  0.0  0.0  0.0  0.0  2.0  0.0 0.0 
 0.0  0.0  0.0  0.0  0.0  0.0  0.8  2.8
1944  0.0  0.0  0.0  0.0  0.0  0.0   0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 0.0 
 0.0  6.6  0.0  0.0  0.0  0.0  0.0  0.0
1945  5.3  0.0  0.0  0.0  0.0  2.5   0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 0.0 
 0.0  5.3  0.0  0.0  0.0  0.0  0.0  0.0
1946  0.0  0.0  0.0  0.0  0.0  2.3   0.0  0.0  0.0  0.0  4.8  0.3  1.5  0.0 0.8 
 0.0  0.0  5.8 70.6 12.4  0.5 23.6  0.0
1947  0.0  0.0  0.0  0.0  0.0  0.0   0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0 0.0 
 2.0  0.0  0.0  0.0  2.8  0.0  0.0  0.0
1948  0.3 20.1  0.0  0.0  0.0  0.0   0.0  0.0  0.0  0.0  1.5  0.5  0.0  0.0 0.0 
 0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0

## To extract all the positive values
x1 - day_data1[,1]  
x2 - day_data1[,2]
x3 - day_data1[,3]
x4 - day_data1[,4]
x5 - day_data1[,5]
tol - 1E-6  
a1 - x1[x1tol]
a2 - x2[x2tol]
a3 - x3[x3tol]
a4 - x4[x4tol]
a5 - x5[x5tol]

library(MASS)
## Example January Charleville 1943-2007
fitdistr(a1,dgamma, list(shape = 1, rate = 0.1), lower = 0.01)
fitdistr(a2,dgamma, list(shape = 1, rate = 0.1), lower = 0.01)
fitdistr(a3,dgamma, list(shape = 1, rate = 0.1), lower = 0.01)
fitdistr(a4,dgamma, list(shape = 1, rate = 0.1), lower = 0.01)
fitdistr(a5,dgamma, list(shape = 1, rate = 0.1), lower = 0.01)

Here is my code: 

alpha.beta - function(data,n)
{  tol - 1E-6
  { for (i in 1:n)
xi - data[,i]
ai - xi [xi  tol]
fit - fitdistr(ai,dgamma, list(shape = 1, rate = 0.1), lower = 0.01)
  }
  fit
}

I’m not sure what went wrong since it gives only one output by right 31 outputs.
Thank you for your attention. 

__
t spam protection around 
http://mail.yahoo.com

__
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 commented, minimal, self-contained, reproducible code.


[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)

 chisq.test(xobs, xtwe = xtwe, rescale.p = TRUE)
Error in chisq.test(xobs, xtwe = xtwe, rescale.p = TRUE) : 
  unused argument(s) (xtwe = c(217.33, 39, 14, 18.33, 6.67, 14.67))
chisq.test(x, p = p, rescale.p = TRUE)
I'm not sure what's wrong with it. 
Thank you so much for your help.


  

You rock. That's why Blockbuster's offering you one month of Blockbuster Total 
Access, No Cost.  
http://tc.deals.yahoo.com/tc/blockbuster/text5.com

__
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 commented, minimal, self-contained, reproducible code.


[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] - qtweedie((rand[!ind]-p0)/(1-p0),p,mu,phi)
   cbind (prob=rand,rain.amt=y)
}
 
tweedie.inv(r1,p0,p,mu,phi)
 
Thank you in advance for your attention.


  

You rock. That's why Blockbuster's offering you one month of Blockbuster Total 
Access, No Cost.  
http://tc.deals.yahoo.com/tc/blockbuster/text5.com

__
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 commented, minimal, self-contained, reproducible code.


[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  0.0  
0.0  4.8  1.6  0.0  1.8  1.8  0.0  2.6
 [24] 33.0 19.0  0.0  0.0  0.6  0.0  0.0  6.8  0.2  0.0  0.0  0.0  0.0  0.2  
0.0  0.0  0.0  0.0  0.0  0.0  1.6  5.4  0.0
 [47]  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.2  0.0  0.0  6.0  6.6  
0.0  0.0  0.0  0.0  0.2  0.0  0.0  0.0  0.0
 [70]  0.0  0.0  0.0  0.0  0.0  0.0  0.0 12.8  0.0  0.0  0.0  0.0  0.0  0.0  
0.2  0.0  3.8  2.2  1.2 12.2  3.0  0.0  0.0
 [93]  0.0  0.0  5.0  3.6  0.0  0.0  0.0  1.8  1.8  1.2  1.2  1.6  3.6  0.2  
0.8  0.0  0.0  0.0  0.0  0.0  1.4 16.6  0.0
[116]  1.6  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  
0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.8  5.4
[139]  1.6  5.8  0.4 34.2  0.0  0.4  0.4  1.4  1.8  0.0  0.0  0.2  0.0  2.4 
24.2  0.4  0.0  0.0  0.0  1.6  0.0  0.0  0.0
[162]  0.0  0.2  0.0  0.0  0.0  4.4  0.0  4.4 21.2  3.0  0.4  0.0  0.0  0.0  
0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
[185]  0.0  0.8  0.8  3.6  0.0  1.2  0.0  0.0  0.0  0.0  0.0  0.0 13.8  1.6  
6.6  5.0  0.0  1.2  0.0  1.0  1.2  2.4  6.0
[208]  2.8  7.0  0.2  0.6  0.0  0.0  0.0  0.0  0.0  3.2  0.2  2.4  0.0  0.0  
0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  2.0
[231]  0.0  0.0  0.0  0.0  0.0  0.0  0.0  1.2  6.6  2.4  1.0  0.0  0.0  0.0  
0.0  0.0  0.0  0.0  0.0  1.6  0.0  0.0  0.0
[254]  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.4  2.6  0.0  0.6  
4.8  0.0  0.2  3.8 17.0  0.2  0.0  0.0  0.0
[277] 43.6  2.6  0.0  0.0  0.0  0.4  8.2 14.4  0.2  0.0  0.0  0.0  1.8  0.0  
0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  1.2
[300]  0.0

Thank you in advance for your help.


  

Looking for last minute shopping deals?

__
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 commented, minimal, self-contained, reproducible code.


[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   25.6
37.8  0.0  0.0 17.6  1.4  0.0  0.0  0.0  0.0  0.0   26.8
42.2  0.8  0.4  0.0  0.2 11.2  1.4 33.2  0.0  0.0   49.4
50.2  1.8  0.0  1.0  0.0  0.2  0.0 12.2  0.0 19.2   34.6
60.0  0.0  0.0  1.0  0.0  0.0  0.0  2.2  0.0 14.6   17.8
70.0  0.0  0.0  0.0  3.6  0.2  0.0  2.0  0.0  0.26.0
80.0  0.0 10.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0   10.0
90.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.00.0
10   1.0  1.4  0.0  0.0  0.0  1.0  0.4  0.0  0.0  0.03.8
11   0.0  8.2  0.0  0.0  0.0  0.0  8.4  0.0  0.0  0.4   17.0
12  10.8  0.8  0.0  0.0  0.0  0.0  1.0  0.0  0.0  4.2   16.8
13  32.8  0.0  0.0  0.8  0.0  0.0  0.0  0.0  0.2  0.2   34.0
14   1.0  0.0  1.6  0.2  0.0  0.0  0.0  0.0  0.0  0.02.8
15   0.0  0.0  0.0  2.2  1.4  0.0  0.0  0.0  0.0  0.03.6
16   0.6  0.0  0.0  0.6 22.0  0.0  0.0  0.0  0.0  0.0   23.2
17   0.0  0.0  0.0  0.0  0.0  0.0  2.8  0.0  0.0  0.02.8
18   2.8  0.0  0.0  0.0  0.0  0.0  8.2  0.0  0.0  8.2   19.2
19   0.0  0.0  0.0  0.0  0.0  0.2  0.0  0.0  0.0  0.20.4
20   0.0  8.2  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.08.2
21   0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.00.0
22   0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.40.4
23   0.0  0.0  0.0  0.0  0.0  0.0  0.0  1.2  0.0  0.01.2
24   0.0  0.0  0.0  0.0  1.0  0.0  0.0  2.0  8.2  0.0   11.2
25   3.2  0.0  0.0  0.6  0.6  0.0  0.0  0.0 11.8  0.0   16.2
26   0.0  0.0 26.2  0.0 12.6  0.0  0.0  2.2  0.0  0.0   41.0
27   0.2  0.0 10.6  0.0  1.2  0.0  0.0  1.8  0.0  0.0   13.8
28   0.0  4.0  0.0  5.8  0.0  0.0  0.0  0.0  0.0  0.09.8
29   0.2 12.4  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0   12.6
30   0.0  2.6  0.0  0.0  2.2  0.0  0.0  0.0  0.0  0.04.8
31   0.0  0.0  0.6  0.0  0.8  0.0  0.0  0.0  4.8  0.06.2
Year Totals 62.8 40.2 51.2 29.8 70.0 12.8 30.8 57.2 55.2 47.6  457.6


  

Be a better friend, newshound, and

__
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 commented, minimal, self-contained, reproducible code.


[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 http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[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
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


<    1   2