The data presented here was obtained through various tests performed on different Shinjei PPD42 units. Each test will be described in more detail below.
# Load libraries
require('openair')
## Loading required package: openair
## Loading required package: lazyeval
## Loading required package: dplyr
##
## Attaching package: 'dplyr'
##
## The following objects are masked from 'package:stats':
##
## filter, lag
##
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
##
## Loading required package: maps
This test was performed outdoors during May 2013 at a residence in Auckland’s North Shore. The data indicate that the sensor did respond to moderate concentrations during the test but did not report above noise levels during expected low concentration periods. There was no reference device during this test so the estimates of low and moderate were made refering to the general conditions throughout the city and may not necessarily reflect local sources impacting the device.
# Loading the data
shinyei_201305 <- read.csv("/data/GUSrepos/R_work/data/shinyei_201305.txt", sep="")
shinyei_201305$date=as.POSIXct(paste(shinyei_201305$Date,shinyei_201305$Time),tz='NZST')
shinyei_201305$Date<-NULL
shinyei_201305$Time<-NULL
summary(shinyei_201305)
## Conc Humdity Temperature
## Min. : 0.0000 Min. :24.10 Min. :18.50
## 1st Qu.: 0.0000 1st Qu.:27.50 1st Qu.:20.40
## Median : 0.0000 Median :30.10 Median :20.70
## Mean : 0.4835 Mean :32.41 Mean :20.68
## 3rd Qu.: 0.0000 3rd Qu.:33.30 3rd Qu.:21.10
## Max. :62.4800 Max. :99.90 Max. :22.10
## date
## Min. :2013-05-17 17:08:24
## 1st Qu.:2013-05-18 08:56:00
## Median :2013-05-19 00:44:00
## Mean :2013-05-19 00:44:00
## 3rd Qu.:2013-05-19 16:32:00
## Max. :2013-05-20 08:20:00
# Plot summaries
timePlot(shinyei_201305,pollutant = c('Conc'),normalise = 'mean', main = 'Raw data\nnormalised to the mean')
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
timePlot(shinyei_201305,pollutant = c('Conc'),normalise = 'mean', avg.time = '10 min',main = '10 min average\nnormalised to the mean')
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
timeVariation(shinyei_201305,pollutant = c('Conc'),normalise = 'TRUE',main = 'Temporal variation\nnormalised to the mean')
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
This test was performed indoors during November 2015 at my desk in a well sealed and air conditioned central Auckland office The data indicate that the sensor did not respond to the relatively low concentrations during the test. There was no reference device during this test so the estimates of low is refered to the expected conditions in this office space.
# Loading the data
shinyei_201511 <- read.csv("/data/GUSrepos/R_work/data/shinyei_201511.txt", sep=",")
shinyei_201511$date=as.POSIXct(shinyei_201511$date)
summary(shinyei_201511)
## date low_pulse ratio concentration
## Min. :2015-11-26 01:22:22 Min. :0 Min. :0 Min. :0.62
## 1st Qu.:2015-11-26 22:47:12 1st Qu.:0 1st Qu.:0 1st Qu.:0.62
## Median :2015-11-27 22:06:33 Median :0 Median :0 Median :0.62
## Mean :2015-11-27 22:07:45 Mean :0 Mean :0 Mean :0.62
## 3rd Qu.:2015-11-28 21:25:54 3rd Qu.:0 3rd Qu.:0 3rd Qu.:0.62
## Max. :2015-11-29 20:45:15 Max. :0 Max. :0 Max. :0.62
# Plot summaries
timePlot(shinyei_201511,pollutant = c('low_pulse','ratio','concentration'),normalise = 'mean', main = 'Raw data\nnormalised to the mean')
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
timePlot(shinyei_201511,pollutant = c('low_pulse','ratio','concentration'),normalise = 'mean', avg.time = '10 min',main = '10 min average\nnormalised to the mean')
## Warning in checkPrep(mydata, vars, type, remove.calm = FALSE): Detected
## data with Daylight Saving Time, converting to UTC/GMT
As there signals from the Shinjei device were constant, the time variation plot is irrelevant.
# timeVariation(shinyei_201511,pollutant = c('low_pulse','ratio','concentration'),normalise = 'TRUE',main = 'Temporal variation\nnormalised to the mean')