[R] HPD credible sets

2005-08-10 Thread pantd
Hi R users Is there a function in R that gives HPD credible sets. i googled it but was in vain! - dev __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] gamma distribution

2005-07-28 Thread pantd
thanks for your response. btw i am calculating the power of the wilcoxon test. i divide the total no. of rejections by the no. of simulations. so for 1000 simulations, at 0.05 level of significance if the no. of rejections are 50 then the power will be 50/1000 = 0.05. thats y im importing in excel

Re: [R] gamma distribution

2005-07-28 Thread pantd
Hi Christopher and Uwe. thanks for your time and guidance. I deeply appreciate it. -dev Quoting Christoph Buser [EMAIL PROTECTED]: Hi As Uwe mentioned be careful about the difference the significance level alpha and the power of a test. To do power calculations you should specify and

Re: [R] gamma distribution

2005-07-27 Thread pantd
Hi You are right but here I am taking into account the p values I get from the tests on the raw and the transformed samples. And then I calculate the power of the tests based on the # of rejections of the p values. DO you think its a good way to determine the power of a test? thanks -dev

Re: [R] gamma distribution

2005-07-27 Thread pantd
Hi I ran your code. I think it should give me the number of p values below 0.05 significance level (thats what i could understand from your code), but after running your code there is neither any error that shows up nor any value that the console displays. thanks in advance -dev. Quoting Uwe

[R] gamma distribution

2005-07-26 Thread pantd
Hi R Users This is a code I wrote and just want to confirm if the first 1000 values are raw gamma (z) and the next 1000 values are transformed gamma (k) or not. As I get 2000 rows once I import into excel, the p - values beyond 1000 dont look that good, they are very high. -- sink(a1.txt);

[R] A Question About Inverse Gamma

2005-07-21 Thread pantd
Hi R users, I am having a little problem finding the the solution to this problem in R: 1. I need to generate normal distribution of sample size 30, mean = 50, sd = 5. 2. From the statistics obtained in step 1, I need to generate the Inverse Gamma distribution. Your views and help will be

[R] Code Verification

2005-07-19 Thread pantd
Hi R Users I have a code which I am running for my thesis work. Just want to make sure that its ok. Its a t test I am conducting between two gamma distributions with different shape parameters. the code looks like: sink(a1.txt); for (i in 1:1000) { x-rgamma(40, 2.5, 10) # n = 40, shape = 2.5,

[R] question from environmental statistics

2005-07-15 Thread pantd
thanks Fran. that was useful but Im still in a fix. its a real life data which looks like this: 0.9 10.9 24.0 6.7 0.6 1.0 2.4 12.4 7.9 15.8 1.4 7.9 11000.0 (benzene conc. taken after WTC attacks)..its just a small chunk of data i pasted for you to look at. its neither normal nor lognormal.

[R] question from environmental statistics

2005-07-14 Thread pantd
Dear R users I want to knw if there is a way in which a raw dataset can be modelled by some distribution. besides the gof test is there any test involving gamma or lognormal that would fit the data. thank you -dev __ R-help@stat.math.ethz.ch mailing