Thanks Robin.  I used the following line, which got my table looking
correct.

Wind_freq_speed <- t(table(cut(Wind_Dir_vec, 0:36), cut(Wind_Speed_vec,
seq(0, to=60, by=10))))

Unfortunately, the rosavent function produced a wind rose diagram that
repeats my frequency data 6 times throughout the 360 degrees and I get the
following error:

Error in legend(-fmaxi - 2 * fint, fmaxi, fill = col, legend = attr(frec,  :

  'legend' is of length 0

Not sure what to do now that my table is set up like the windfreq.dat table
used in the rosavent example.  I think I may check out some of these
auxillary controls for rose diagrams.  Any thoughts are welcome.

Thanks.


On Sat, Apr 17, 2010 at 2:24 AM, Robin Evans <rje42@stat.washington.edu>wrote:

> > but I need to incorporate the speed values and have a table like the
> > following:
> >                 1     2     3     4     5     6 .... 36
> > 0-9
> > 10-19
> > 20-29
> > 30-39
> > 40-49
> > 50+
> > The final table will be used in the rosavent function to produce a wind
> rose
> > diagram just like in this example:
>
> The table() command will accept an arbitrary number of factors as its
> first arguments, and cross tabulate the frequencies for all of them.
> So
>
> table(cut(Wind_Dir_vec, 0:36), cut(Wind_Speed, seq(9.5, to=49.5, by=10)))
>
> should do more or less what you want.
>
> > --
> > Kevin Turner
> > Department of Geography
> > and Environmental Studies
> > Wilfrid Laurier University
> > Waterloo, Ontario
> > [[alternative HTML version deleted]]
>
> --
> Robin Evans
> Statistics Department
> University of Washington
> www.stat.washington.edu/~rje42 <http://www.stat.washington.edu/%7Erje42>
>



-- 
Kevin Turner
Department of Geography
and Environmental Studies
Wilfrid Laurier University
Waterloo, Ontario

	[[alternative HTML version deleted]]

