topic Re: Count Data in Statistical Procedures

Count Data

vdfdd — Thu, 12 Jul 2018 14:28:29 GMT

Please I need suggestions on how to simulate count data from a Poisson distribution and fit (using PROC GLIMMIX) Normal, lognormal and negbin to the resulting data.Thank you

Re: Count Data

Reeza — Thu, 12 Jul 2018 14:41:22 GMT

https://www.amazon.ca/Simulating-Data-SAS-Rick-Wicklin-ebook/dp/B00CC0DQ16

Re: Count Data

StatDave — Thu, 12 Jul 2018 14:53:35 GMT

Use the RAND('POISSON', mean) function in the DATA step to generate Poisson distributed data. Just set the desired value of the Poisson mean to produce a set of values. For example, these statements produce 100 observations in which Y has Poisson distributed values with mean=2.

data s;
do i=1 to 100;
y=rand("poisson",2);
output;
end;
run;

Re: Count Data

PeterClemmensen — Thu, 12 Jul 2018 14:59:29 GMT

Are you sure PROC GLIMMIX is the right tool for fitting this data?

Re: Count Data

vdfdd — Thu, 12 Jul 2018 15:31:08 GMT

Thank you for your swift response.

Re: Count Data

vdfdd — Thu, 12 Jul 2018 15:32:26 GMT

I am not completely sure. I was looking at the example in this link: https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_glimmix_sect026.htm

Re: Count Data

Reeza — Thu, 12 Jul 2018 15:38:20 GMT

What are you trying to do overall? Why the need to simulate data?

@vdfdd wrote:

I am not completely sure. I was looking at the example in this link: https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_glimmix_sect026.htm

Re: Count Data

vdfdd — Thu, 12 Jul 2018 15:47:31 GMT

To test which distribution best fit the (non-normal, discrete) data

Re: Count Data

Reeza — Thu, 12 Jul 2018 16:06:30 GMT

PROC UNIVARIATE can test this with histogram statements.

Re: Count Data

Ksharp — Fri, 13 Jul 2018 12:47:54 GMT

You could find answer at @Rick_SAS blog.

http://blogs.sas.com/content/iml/

Re: Count Data

vdfdd — Fri, 13 Jul 2018 14:51:02 GMT

Thank you for your suggestion.