R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Panel's performance according to its capabilities to...

panelperf

R Documentation

Panel's performance according to its capabilities to dicriminate between products

Description

Computes automatically P-values associated with the F-test as well as the residual term for a given analysis of variance model.

Usage

panelperf(donnee, formul, subset = NULL, firstvar, 
    lastvar = ncol(donnee), random = TRUE)

Arguments

`donnee`	a data frame
`formul`	the model that is to be tested
`subset`	cf. function `lm` in the stats package
`firstvar`	the position of the first endogenous variable
`lastvar`	the position of the last endogenous variable (by default the last column of `donnee`
`random`	boolean, effect should be possible as fixed or random (default as random)

Details

The formul parameter must be filled in by an analysis of variance model and must begin with the categorical variable of interest (e.g. the product effect) followed by the different other factors of interest (and their combinations). E.g.:formul = "~Product+Session".

Value

A list containing the following components:

`p.value`	a matrix of dimension (k,m) of P-values associated with the F-test for the k descriptors and the m factors and their combinations considered in the analysis of variance model of interest
`variability`	a matrix of dimension (k,m) where the entries correspond to the percentages of variability due to the effects introduced in the analysis of variance model of interest
`res`	a vector of dimension k of residual terms for the analysis of variance model of interest
`r2`	a vector of dimension k of r-squared for the analysis of variance model of interest

Author(s)

F Husson, S Le

References

P. Lea, T. Naes, M. Rodbotten. Analysis of variance for sensory data.
H. Sahai, M. I. Ageel. The analysis of variance.

Examples

data(chocolates)
res=panelperf(sensochoc, firstvar = 5, formul = "~Product+Panelist+
    Session+Product:Panelist+Session:Product+Panelist:Session")
## Sort results by product p.values.
coltable(magicsort(res$p.value, sort.mat = res$p.value[,1], bycol = FALSE,
    method = "median"), main.title = "Panel performance (sorted by product P-value)")

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(SensoMineR)
Loading required package: FactoMineR
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/SensoMineR/panelperf.Rd_%03d_medium.png", width=480, height=480)
> ### Name: panelperf
> ### Title: Panel's performance according to its capabilities to dicriminate
> ###   between products
> ### Aliases: panelperf
> ### Keywords: models
> 
> ### ** Examples
> 
> data(chocolates)
> res=panelperf(sensochoc, firstvar = 5, formul = "~Product+Panelist+
+     Session+Product:Panelist+Session:Product+Panelist:Session")
> ## Sort results by product p.values.
> coltable(magicsort(res$p.value, sort.mat = res$p.value[,1], bycol = FALSE,
+     method = "median"), main.title = "Panel performance (sorted by product P-value)")
Error in dev.new() : no suitable unused file name for pdf()
Calls: coltable -> police -> dev.new
Execution halted