Last data update: 2014.03.03

R: Summarized magnitude distribution
mag.distrR Documentation

Summarized magnitude distribution

Description

Table and graphical representation of summarized magnitude distribution for a given magnitude dataset, specified meteor shower and time period.

Usage

mag.distr(data,year,month.beg,month.end=month.beg,day.beg,day.end=day.beg,
time.beg=0,time.end=2359,shw)

Arguments

data

data frame consisting of visual meteor magnitude data.

year

numeric vector of length 4 specifying year.

month.beg

numeric vector specifying the beginning month.

month.end

numeric vector specifying the ending month. By default, month.end is set to be equal to month.beg.

day.beg

numeric vector specifying the beginning day.

day.end

numeric vector specifying the ending day. By default, day.end is set to be equal to day.beg.

time.beg

numeric vector (0-2359) specifying lower boundary of time in hours and minutes, corresponding to day.beg. By default, time.beg is set to be equal to 0.

time.end

numeric vector(0-2359) specifying upper boundary of time in hours and minutes, corresponding to day.end. By default, time.end is set to be equal to 2359.

shw

character string consisting of three capital letters which represent meteor shower code.

Details

Summarized magnitude distribution is formed by summing frequencies of all observers for each magnitude value.

Value

Table and plot of summarized magnitude distribution consisting of histogram and boxplot.

The histogram cells are intervals of length 1, with midpoints at magnitude values.

Author(s)

Kristina Veljkovic

References

http://www.imo.net/data/visual

See Also

pop.index

Examples

## select data for observations of Perseids, period 12-14th August  2007 
## and make graphics of magnitude distribution
data(magn07)
magnPer<-filter(magn07,shw="PER", year=2007, month.beg=8, day.beg=12, day.end=14)
mag.distr(magnPer,year=2007, month.beg=8, day.beg=12, day.end=14, shw="PER")

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(MetFns)
Loading required package: astroFns
Loading required package: fBasics
Loading required package: timeDate
Loading required package: timeSeries


Rmetrics Package fBasics
Analysing Markets and calculating Basic Statistics
Copyright (C) 2005-2014 Rmetrics Association Zurich
Educational Software for Financial Engineering and Computational Science
Rmetrics is free software and comes with ABSOLUTELY NO WARRANTY.
https://www.rmetrics.org --- Mail to: info@rmetrics.org
Loading required package: lubridate

Attaching package: 'lubridate'

The following object is masked from 'package:base':

    date

Loading required package: plotrix
Loading required package: spatial

Attaching package: 'MetFns'

The following object is masked from 'package:timeSeries':

    filter

The following object is masked from 'package:stats':

    filter

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/MetFns/mag.distr.Rd_%03d_medium.png", width=480, height=480)
> ### Name: mag.distr
> ### Title: Summarized magnitude distribution
> ### Aliases: mag.distr
> 
> ### ** Examples
> 
> ## select data for observations of Perseids, period 12-14th August  2007 
> ## and make graphics of magnitude distribution
> data(magn07)
> magnPer<-filter(magn07,shw="PER", year=2007, month.beg=8, day.beg=12, day.end=14)
> mag.distr(magnPer,year=2007, month.beg=8, day.beg=12, day.end=14, shw="PER")
     m6 m5   m4 m3   m2    m1 zero  p1  p2     p3    p4  p5   p6 p7
[1,]  4  5 15.5 36 98.5 203.5  414 633 942 1007.5 734.5 380 65.5  0
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>