pop.index.plot: Plot of population index

Description

Plots population index of a meteor shower for data containing population index values and given limits on x and y axes.

Usage

pop.index.plot(rdata,xlim1,xlim2,xinc,ylim1,ylim2,yinc,dlim1=NULL,dlim2=NULL,
dinc=NULL,dunit=NULL)

Arguments

rdata

data frame consisting of population index values.

xlim1

numeric vector specifying minimum value on x-axis.

xlim2

numeric vector specifying maximum value on x-axis.

xinc

numeric vector specifying increment between labels on x-axis.

ylim1

numeric vector specifying minimum value on y-axis.

ylim2

numeric vector specifying maximum value on y-axis.

yinc

numeric vector specifying increment between labels on y-axis.

dlim1

character vector specifying start date on secondary x-axis, given in "%Y-%m-%d %H:%M:%S" format (UTC timezone).

dlim2

character vector specifying end date on secondary x-axis, given in "%Y-%m-%d %H:%M:%S" format (UTC timezone).

dinc

numeric vector specifying increment between labels on secondary x-axis.

dunit

character vector specifying unit of increment dinc. Accepted values are "min" for minute, "h" for hour and "day" for day.

Value

xy plot of population index is made. Solar longitude is on x-axis and population index on y-axis. Population index is represented with black filled circles and 68% confidence intervals. Values of limits (minimum and maximum values) on x and y axis (xlim1,xlim2,ylim1,ylim2), as well as increments between the axis labels (xinc,yinc) should be provided to function call.

Details

rdata represents table of population index values calculated using functions pop.index or pop.index2.

References

Arlt, R. (2003). Bulletin 19 of the International Leonid Watch: Population index study of the 2002 Leonid meteors. WGN, Journal of the IMO,31:3, 77-87.

Examples

Run this code

# NOT RUN {
## calculate and make plot of population index for observations of 2015 Orionids,
## time period around maximum, 21th to 25th October, 
## min bin size 0.16 degrees (4 hours), max bin size 1 degree, 
## number of meteors equals 100 
## First select magnitude data for Orionids activity from 21th to 25th October, 
## limiting magnitudes above 5, radiant elevation above 20 degrees, 
## percentage of clouds below 20
ori2015mag<-filter(magn2015,date.start="2015-10-21",date.end="2015-10-25",
shw="ORI",P.up=20,mag.low=5,h.low=20)
oripop<-pop.index2(ori2015mag,date.start="2015-10-21",date.end="2015-10-25",
shw="ORI",kmin=0.16,kmax=1,num=100)
pop.index.plot(oripop,xlim1=207,xlim2=211,xinc=1,ylim1=1.6,ylim2=2.4,yinc=0.2)

## change date tick marks on secondary x-axis, to show date at midnight
pop.index.plot(oripop,xlim1=207,xlim2=211,xinc=1,ylim1=1.5,ylim2=2.5,
yinc=0.2,dlim1="2015-10-21",dlim2="2015-10-25",dinc=1,dunit="day")
# }

Run the code above in your browser using DataLab

Data Engineering and BI courses are free this week!