Last data update: 2014.03.03

R: Plotting method for 'lab.fqcs' objects
plot.lab.fqcsR Documentation

Plotting method for 'lab.fqcs' objects

Description

Generic function to plot objects of 'lab.fqcs' class. Results of functional ILS studies are graphically shown.

Usage

## S3 method for class 'lab.fqcs'
plot(x, main = NULL, xlab = NULL, ylab = NULL,
  ylim = NULL, x.co = NULL, y.co = NULL, legend = TRUE, col = NULL,
  ...)

Arguments

x

Object functional data or a list with objects of functional data type

main

Main title for the plot

xlab

Title for the x axis

ylab

Title for the y axis

ylim

The y limits of the plot

x.co

It speficies the x co-ordinates to be used to place a legend.

y.co

It specifies the y co-ordinates to be used to place a legend.

legend

Logical argument. Default is TRUE then The legend default is used.

col

Color specifications

...

arguments to be passed to or from methods

References

Febrero-Bande, M. and Oviedo, M. (2012), "Statistical computing in functional data analysis: the R package fda.usc". Journal of Statistical Software 51 (4), 1-28.

Naya, S., Tarrio-Saavedra. J., Lopez- Beceiro, J., Francisco Fernandez, M., Flores, M. and Artiaga, R. (2014), "Statistical functional approach for interlaboratory studies with thermal data". Journal of Thermal Analysis and Calorimetry, 118,1229-1243.

Examples

library(ILS)
data(TG)
delta <- seq(from = 40 ,to = 850 ,length.out = 1000 )
curves.fqcd <- lab.fqcd(TG, argvals = delta)
curves.fqcs <- lab.fqcs(curves.fqcd)
summary(curves.fqcs)
names(curves.fqcs)
class(curves.fqcs$mean.i)
xlab <- "Temperature/ C"
ylab <- "Mass/ %"
main <- "Functional Mean Estimation by Laboratory"
p <- dim(curves.fqcd$curves)[3]
col <- 1:p
plot(curves.fqcs$mean.i,main = main, xlab = xlab, ylab = ylab, col = col,legend = FALSE)
legend(45,70,paste("Lab",1:p),
      col = col,lty = 1,lwd = 2,cex = 0.7)

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(ILS)
Loading required package: multcomp
Loading required package: mvtnorm
Loading required package: survival
Loading required package: TH.data
Loading required package: MASS

Attaching package: 'TH.data'

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

    geyser

Loading required package: depthTools
Loading required package: fda.usc
Loading required package: fda
Loading required package: splines
Loading required package: Matrix

Attaching package: 'fda'

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

    matplot

Loading required package: mgcv
Loading required package: nlme
This is mgcv 1.8-12. For overview type 'help("mgcv-package")'.
Loading required package: rpart

 Package ILS: Interlaboratory Study 
 version 0.1.0 (built on 2016-05-22).
 Copyright Miguel A. Flores Sanchez 2016. 


Attaching package: 'ILS'

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

    Glucose

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ILS/plot.lab.fqcs.Rd_%03d_medium.png", width=480, height=480)
> ### Name: plot.lab.fqcs
> ### Title: Plotting method for 'lab.fqcs' objects
> ### Aliases: plot.lab.fqcs
> 
> ### ** Examples
> 
> library(ILS)
> data(TG)
> delta <- seq(from = 40 ,to = 850 ,length.out = 1000 )
> curves.fqcd <- lab.fqcd(TG, argvals = delta)
> curves.fqcs <- lab.fqcs(curves.fqcd)
> summary(curves.fqcs)

Number of laboratories:  7
Number of replicates:  15
Summary of fdata list:
             Length Class      Mode   
curves       105000 array      numeric
curves.fdata      7 list.fdata list   
p                 1 -none-     numeric
n                 1 -none-     numeric
m                 1 -none-     numeric

Summary of mean fdata list:
             Length Class      Mode   
curves       7000   matrix     numeric
curves.fdata    7   list.fdata list   
p               1   -none-     numeric
n               1   -none-     numeric
m               1   -none-     numeric

Summary of var fdata list:
             Length Class      Mode   
curves       7000   matrix     numeric
curves.fdata    7   list.fdata list   
p               1   -none-     numeric
n               1   -none-     numeric
m               1   -none-     numeric
> names(curves.fqcs)
[1] "lab.fqcd" "mean.i"   "s2.i"     "mean"     "S2"       "p"        "n"       
[8] "m"       
> class(curves.fqcs$mean.i)
[1] "lab.fqcd"
> xlab <- "Temperature/ C"
> ylab <- "Mass/ %"
> main <- "Functional Mean Estimation by Laboratory"
> p <- dim(curves.fqcd$curves)[3]
> col <- 1:p
> plot(curves.fqcs$mean.i,main = main, xlab = xlab, ylab = ylab, col = col,legend = FALSE)
> legend(45,70,paste("Lab",1:p),
+       col = col,lty = 1,lwd = 2,cex = 0.7)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>