Analysis of dataset 0 (first measurement in lab/testing)

Dataset

  • simultaneous operation of the two scionix sensors (joint measurement test)
  • measurement starting in 2nd December 2016
  • output extracted from datalogger using PC200W
  • header and first two records (zero) are discarded (# before the corresponding line)

File: dataset-0

In [6]:
data0=read.table(file="Dataset-0files/CR800Series_Table1.dat",sep=",")

str(data0)
'data.frame':	3602 obs. of  8 variables:
 $ V1: Factor w/ 3602 levels "2016-12-02 16:17:00",..: 1 2 3 4 5 6 7 8 9 10 ...
 $ V2: int  2 3 4 5 6 7 8 9 10 11 ...
 $ V3: num  12 12 12 12 12 ...
 $ V4: num  12 12 12 12 12 ...
 $ V5: num  22.6 22.6 22.6 22.6 22.6 ...
 $ V6: num  22.4 22.2 22.2 22.2 22.2 ...
 $ V7: int  4374 4456 4339 4461 4529 4478 4278 4448 4436 4479 ...
 $ V8: int  4534 4679 4590 4708 4515 4596 4591 4649 4545 4398 ...

Exploratory data analysis

In [29]:
library(chron)

dat=as.POSIXct(as.character(data0$V1), format="%Y-%m-%d %H:%M:%S", tz="UTC")

sci1=data0$V7
sci2=data0$V8

options(repr.plot.width=8, repr.plot.height=5)

plot(dat,sci1,type="b",pch=20,xlab="time",ylab="counts", main="all measurements")
lines(dat,sci2,col="darkgrey",pch=20,type="b")

Zoom-in

From record number 39 (line 44) to 3586 (line 3591)
(12-05 10:47 to 12-07 21:54)

In [57]:
system("head -3591 Dataset-0files/CR800Series_Table1.dat | tail -3548 > data0Zoom.txt", intern=TRUE)
In [61]:
dataZoom=read.table(file="Dataset-0files/data0Zoom.txt",sep=",")
In [62]:
datz=as.POSIXct(as.character(dataZoom$V1), format="%Y-%m-%d %H:%M:%S", tz="UTC")

sci1=dataZoom$V7
sci2=dataZoom$V8

options(repr.plot.width=8, repr.plot.height=5)

plot(datz,sci1,type="b",pch=20,xlab="time",ylab="counts", main="scionix 1")
plot(datz,sci2,type="b",pch=20,xlab="time",ylab="counts", main="scionix 2")