Last data update: 2014.03.03

 R: Discriptives Follow Groups
DgroupR Documentation

Discriptives Follow Groups

Description

This function descriptive statistics in group for a continous variable. Usual using to statistic a time series following dates in week or months in year.

Usage

Dgroup(x, follow = NULL, r = 2, answer = 1, statistic = "ALL")

Arguments

x

A continous variable or a time series.

follow

A factor or a list factor which contain not too two factors.

r

Rounds the answer to the specified number of decimal places (default 2).

answer

Form of answers are returned. Let answer=1 or answer=2 (default 1).

statistic

A list contain descriptive statistic values that user want R print screen (default ALL).

Details

Data is divided into groups by follow and then every group are calculated by Descriptives function.

Value

N

Length sample

NaN

Number NA values

Min

Min value

1sq QU

Value in 25% of interval probabilities

Median

Median value

Mean

Mean value

3rd QU

Value in 75% of interval probabilities

Max

Max value

VAR

Variance value

SD

Standard Deviation

SE

Standard Deviation of the Estimated Means

Note

The function just maximum calculated for two factors.

Author(s)

Hong Viet Minh <hongvietminh@gmail.com>

References

Theory of base statistic.

See Also

Descriptives,Frequencies

Examples

#Factor date
date<-as.factor(c("Tue","Wed","Thu","Fri","Mon","Tue","Wed","Thu","Fri","Mon",
"Tue","Wed","Thu","Fri","Mon","Tue","Wed","Thu","Fri","Mon","Tue","Wed","Thu",
"Fri","Mon","Tue","Wed","Thu","Fri","Mon","Tue","Wed","Thu"))

#Factor hk
hk<-as.factor(c("hk1","hk2","hk1","hk3","hk3","hk1","hk1","hk1","hk2","hk2","hk2",
"hk1","hk2","hk1","hk1","hk1","hk2","hk1","hk1","hk1","hk1","hk2","hk1",
"hk1","hk1","hk1","hk3","hk1","hk3","hk3","hk2","hk3","hk1"))

#A continous variable
coffee<-c(5,6,8,4,3,7,6,0,3,2,3,4,9,1,3,8,7,8,2,3,8,6,4,4,6,7,6,5,2,3,8,4,4)

#Descriptive statistics in group
Dgroup(coffee,r=4,answer=2)
Dgroup(coffee,follow=list(date),r=4)
Dgroup(coffee,follow=date,r=4,answer=2)
Dgroup(coffee,follow=date,r=4,statistic=list("Mean","Max"))
Dgroup(coffee,r=4,follow=list(date,hk),answer=1)
Dgroup(coffee,r=4,follow=list(date,hk),answer=2)
Dgroup(coffee,r=4,follow=list(hk,date),answer=1)
Dgroup(coffee,r=4,follow=list(hk,date),answer=2)

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(AnalyzeTS)
Loading required package: MASS
Loading required package: TSA
Loading required package: leaps
Loading required package: locfit
locfit 1.5-9.1 	 2013-03-22
Loading required package: mgcv
Loading required package: nlme
This is mgcv 1.8-12. For overview type 'help("mgcv-package")'.
Loading required package: tseries

Attaching package: 'TSA'

The following objects are masked from 'package:stats':

    acf, arima

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

    tar

Loading required package: TTR
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/AnalyzeTS/Dgroup.Rd_%03d_medium.png", width=480, height=480)
> ### Name: Dgroup
> ### Title: Discriptives Follow Groups
> ### Aliases: Dgroup
> ### Keywords: Dgroup
> 
> ### ** Examples
> 
> #Factor date
> date<-as.factor(c("Tue","Wed","Thu","Fri","Mon","Tue","Wed","Thu","Fri","Mon",
+ "Tue","Wed","Thu","Fri","Mon","Tue","Wed","Thu","Fri","Mon","Tue","Wed","Thu",
+ "Fri","Mon","Tue","Wed","Thu","Fri","Mon","Tue","Wed","Thu"))
> 
> #Factor hk
> hk<-as.factor(c("hk1","hk2","hk1","hk3","hk3","hk1","hk1","hk1","hk2","hk2","hk2",
+ "hk1","hk2","hk1","hk1","hk1","hk2","hk1","hk1","hk1","hk1","hk2","hk1",
+ "hk1","hk1","hk1","hk3","hk1","hk3","hk3","hk2","hk3","hk1"))
> 
> #A continous variable
> coffee<-c(5,6,8,4,3,7,6,0,3,2,3,4,9,1,3,8,7,8,2,3,8,6,4,4,6,7,6,5,2,3,8,4,4)
> 
> #Descriptive statistics in group
> Dgroup(coffee,r=4,answer=2)
  N: NaN: Min: 1sq QU: Median:  Mean: 3rd QU: Max:   VAR:    SD:   SE:
x 33    0    0       3       4 4.8182       7    9 5.4659 2.3379 0.407
> Dgroup(coffee,follow=list(date),r=4)
           Tue    Wed    Thu    Fri    Mon
N:      7.0000 7.0000 7.0000 6.0000 6.0000
NaN:    0.0000 0.0000 0.0000 0.0000 0.0000
Min:    3.0000 4.0000 0.0000 1.0000 2.0000
1sq QU: 6.0000 5.0000 4.0000 2.0000 3.0000
Median: 7.0000 6.0000 5.0000 2.5000 3.0000
Mean:   6.5714 5.5714 5.4286 2.6667 3.3333
3rd QU: 8.0000 6.0000 8.0000 3.7500 3.0000
Max:    8.0000 7.0000 9.0000 4.0000 6.0000
VAR:    3.6190 1.2857 9.9524 1.4667 1.8667
SD:     1.9024 1.1339 3.1547 1.2111 1.3663
SE:     0.7190 0.4286 1.1924 0.4944 0.5578
> Dgroup(coffee,follow=date,r=4,answer=2)
    N: NaN: Min: 1sq QU: Median:  Mean: 3rd QU: Max:   VAR:    SD:    SE:
Tue  7    0    3       6     7.0 6.5714    8.00    8 3.6190 1.9024 0.7190
Wed  7    0    4       5     6.0 5.5714    6.00    7 1.2857 1.1339 0.4286
Thu  7    0    0       4     5.0 5.4286    8.00    9 9.9524 3.1547 1.1924
Fri  6    0    1       2     2.5 2.6667    3.75    4 1.4667 1.2111 0.4944
Mon  6    0    2       3     3.0 3.3333    3.00    6 1.8667 1.3663 0.5578
> Dgroup(coffee,follow=date,r=4,statistic=list("Mean","Max"))
         Tue    Wed    Thu    Fri    Mon
Mean: 6.5714 5.5714 5.4286 2.6667 3.3333
Max:  8.0000 7.0000 9.0000 4.0000 6.0000
> Dgroup(coffee,r=4,follow=list(date,hk),answer=1)
$Fri
           hk3 hk2    hk1
N:      2.0000   1 3.0000
NaN:    0.0000   0 0.0000
Min:    2.0000   3 1.0000
1sq QU: 2.5000   3 1.5000
Median: 3.0000   3 2.0000
Mean:   3.0000   3 2.3333
3rd QU: 3.5000   3 3.0000
Max:    4.0000   3 4.0000
VAR:    2.0000  NA 2.3333
SD:     1.4142  NA 1.5275
SE:     1.0000  NA 0.8819

$Mon
        hk3 hk2    hk1
N:        2   1 3.0000
NaN:      0   0 0.0000
Min:      3   2 3.0000
1sq QU:   3   2 3.0000
Median:   3   2 3.0000
Mean:     3   2 4.0000
3rd QU:   3   2 4.5000
Max:      3   2 6.0000
VAR:      0  NA 3.0000
SD:       0  NA 1.7321
SE:       0  NA 1.0000

$Thu
           hk1 hk2
N:      6.0000   1
NaN:    0.0000   0
Min:    0.0000   9
1sq QU: 4.0000   9
Median: 4.5000   9
Mean:   4.8333   9
3rd QU: 7.2500   9
Max:    8.0000   9
VAR:    8.9667  NA
SD:     2.9944  NA
SE:     1.2225  NA

$Tue
           hk1     hk2
N:      5.0000  2.0000
NaN:    0.0000  0.0000
Min:    5.0000  3.0000
1sq QU: 7.0000  4.2500
Median: 7.0000  5.5000
Mean:   7.0000  5.5000
3rd QU: 8.0000  6.7500
Max:    8.0000  8.0000
VAR:    1.5000 12.5000
SD:     1.2247  3.5355
SE:     0.5477  2.5000

$Wed
           hk2    hk1    hk3
N:      3.0000 2.0000 2.0000
NaN:    0.0000 0.0000 0.0000
Min:    6.0000 4.0000 4.0000
1sq QU: 6.0000 4.5000 4.5000
Median: 6.0000 5.0000 5.0000
Mean:   6.3333 5.0000 5.0000
3rd QU: 6.5000 5.5000 5.5000
Max:    7.0000 6.0000 6.0000
VAR:    0.3333 2.0000 2.0000
SD:     0.5774 1.4142 1.4142
SE:     0.3333 1.0000 1.0000

> Dgroup(coffee,r=4,follow=list(date,hk),answer=2)
$Fri
    N: NaN: Min: 1sq QU: Median:  Mean: 3rd QU: Max:   VAR:    SD:    SE:
hk3  2    0    2     2.5       3 3.0000     3.5    4 2.0000 1.4142 1.0000
hk2  1    0    3     3.0       3 3.0000     3.0    3     NA     NA     NA
hk1  3    0    1     1.5       2 2.3333     3.0    4 2.3333 1.5275 0.8819

$Mon
    N: NaN: Min: 1sq QU: Median: Mean: 3rd QU: Max: VAR:    SD: SE:
hk3  2    0    3       3       3     3     3.0    3    0 0.0000   0
hk2  1    0    2       2       2     2     2.0    2   NA     NA  NA
hk1  3    0    3       3       3     4     4.5    6    3 1.7321   1

$Thu
    N: NaN: Min: 1sq QU: Median:  Mean: 3rd QU: Max:   VAR:    SD:    SE:
hk1  6    0    0       4     4.5 4.8333    7.25    8 8.9667 2.9944 1.2225
hk2  1    0    9       9     9.0 9.0000    9.00    9     NA     NA     NA

$Tue
    N: NaN: Min: 1sq QU: Median: Mean: 3rd QU: Max: VAR:    SD:    SE:
hk1  5    0    5    7.00     7.0   7.0    8.00    8  1.5 1.2247 0.5477
hk2  2    0    3    4.25     5.5   5.5    6.75    8 12.5 3.5355 2.5000

$Wed
    N: NaN: Min: 1sq QU: Median:  Mean: 3rd QU: Max:   VAR:    SD:    SE:
hk2  3    0    6     6.0       6 6.3333     6.5    7 0.3333 0.5774 0.3333
hk1  2    0    4     4.5       5 5.0000     5.5    6 2.0000 1.4142 1.0000
hk3  2    0    4     4.5       5 5.0000     5.5    6 2.0000 1.4142 1.0000

> Dgroup(coffee,r=4,follow=list(hk,date),answer=1)
$hk1
           Tue    Thu    Wed    Fri    Mon
N:      5.0000 6.0000 2.0000 3.0000 3.0000
NaN:    0.0000 0.0000 0.0000 0.0000 0.0000
Min:    5.0000 0.0000 4.0000 1.0000 3.0000
1sq QU: 7.0000 4.0000 4.5000 1.5000 3.0000
Median: 7.0000 4.5000 5.0000 2.0000 3.0000
Mean:   7.0000 4.8333 5.0000 2.3333 4.0000
3rd QU: 8.0000 7.2500 5.5000 3.0000 4.5000
Max:    8.0000 8.0000 6.0000 4.0000 6.0000
VAR:    1.5000 8.9667 2.0000 2.3333 3.0000
SD:     1.2247 2.9944 1.4142 1.5275 1.7321
SE:     0.5477 1.2225 1.0000 0.8819 1.0000

$hk2
           Wed Fri Mon     Tue Thu
N:      3.0000   1   1  2.0000   1
NaN:    0.0000   0   0  0.0000   0
Min:    6.0000   3   2  3.0000   9
1sq QU: 6.0000   3   2  4.2500   9
Median: 6.0000   3   2  5.5000   9
Mean:   6.3333   3   2  5.5000   9
3rd QU: 6.5000   3   2  6.7500   9
Max:    7.0000   3   2  8.0000   9
VAR:    0.3333  NA  NA 12.5000  NA
SD:     0.5774  NA  NA  3.5355  NA
SE:     0.3333  NA  NA  2.5000  NA

$hk3
           Fri Mon    Wed
N:      2.0000   2 2.0000
NaN:    0.0000   0 0.0000
Min:    2.0000   3 4.0000
1sq QU: 2.5000   3 4.5000
Median: 3.0000   3 5.0000
Mean:   3.0000   3 5.0000
3rd QU: 3.5000   3 5.5000
Max:    4.0000   3 6.0000
VAR:    2.0000   0 2.0000
SD:     1.4142   0 1.4142
SE:     1.0000   0 1.0000

> Dgroup(coffee,r=4,follow=list(hk,date),answer=2)
$hk1
    N: NaN: Min: 1sq QU: Median:  Mean: 3rd QU: Max:   VAR:    SD:    SE:
Tue  5    0    5     7.0     7.0 7.0000    8.00    8 1.5000 1.2247 0.5477
Thu  6    0    0     4.0     4.5 4.8333    7.25    8 8.9667 2.9944 1.2225
Wed  2    0    4     4.5     5.0 5.0000    5.50    6 2.0000 1.4142 1.0000
Fri  3    0    1     1.5     2.0 2.3333    3.00    4 2.3333 1.5275 0.8819
Mon  3    0    3     3.0     3.0 4.0000    4.50    6 3.0000 1.7321 1.0000

$hk2
    N: NaN: Min: 1sq QU: Median:  Mean: 3rd QU: Max:    VAR:    SD:    SE:
Wed  3    0    6    6.00     6.0 6.3333    6.50    7  0.3333 0.5774 0.3333
Fri  1    0    3    3.00     3.0 3.0000    3.00    3      NA     NA     NA
Mon  1    0    2    2.00     2.0 2.0000    2.00    2      NA     NA     NA
Tue  2    0    3    4.25     5.5 5.5000    6.75    8 12.5000 3.5355 2.5000
Thu  1    0    9    9.00     9.0 9.0000    9.00    9      NA     NA     NA

$hk3
    N: NaN: Min: 1sq QU: Median: Mean: 3rd QU: Max: VAR:    SD: SE:
Fri  2    0    2     2.5       3     3     3.5    4    2 1.4142   1
Mon  2    0    3     3.0       3     3     3.0    3    0 0.0000   0
Wed  2    0    4     4.5       5     5     5.5    6    2 1.4142   1

> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>


Created & Maintained by Osamu Ogasawara (osamu.ogasawara@gmail.com) and IMSBIO