Google Groups no longer supports new Usenet posts or subscriptions. Historical content remains viewable.
Dismiss

geometric mean month by month

4 views
Skip to first unread message

Saad

unread,
May 23, 2013, 1:56:11 PM5/23/13
to
Dear All,

I need some help please.
I have two columns in a matrix B.
B=[A, X];
X is a column vector that contains different daily observations.
A is a column vector that contains the number of the month (for example for January it has "1" for february "2" etc till December where it has "12").

Column X has daily data for January, Feb till december. Thus Column A has numbers that change according to the month number. For instance, for january, Column A will have thirty one "1" corresponding to 31 days in january and twenty nine "2" corresponding to 29 days in February until the end of the year then i will start again with january with "1" and february with "2" etc till december .

What I would like to do is to do the gemetric mean of the values of X for similar values of A. Basically, it is like if Iam doing a geometric mean month by month untill the end of my data series. How can I do that please? I am stuck in this. ANy help is much appreciated.
Best

S

dpb

unread,
May 24, 2013, 10:11:49 AM5/24/13
to
On 5/23/2013 12:56 PM, Saad wrote:
...
> I need some help please.
> I have two columns in a matrix B.
> B=[A, X];
> X is a column vector that contains different daily observations.
> A is a column vector that contains the number of the month (for example
> for January it has "1" for february "2" etc till December where it has
> "12").
>
...

> What I would like to do is to do the gemetric mean of the values of X
> for similar values of A. Basically, it is like if Iam doing a geometric
> mean month by month untill the end of my data series. How can I do that
> please? I am stuck in this. ANy help is much appreciated.

"Deadahead" solution...

g=zeros(12,1);
for i=1:12
g(i)=geomean(B(B(1,:)==i,2);
end

I've only very recently received my later version of ML that has
arrayfun and such modern niceties so my fluency therein isn't yet up to
writing the above on the fly via it but that's the "neat" way...

--

0 new messages