FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_cloud.wasp
Title produced by softwareTrivariate Scatterplots
Date of computationThu, 01 May 2008 20:03:09 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/May/02/t1209693918fib3g37kci70kvi.htm/, Retrieved Wed, 15 May 2024 16:05:51 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=11146, Retrieved Wed, 15 May 2024 16:05:51 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact205

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Trivariate Scatterplots] [Venus 0 SM] [2008-05-02 02:03:09] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
0.73
0.7
0.61
0.47
0.29
0.09
-0.12
-0.32
-0.5
-0.63
-0.7
-0.72
-0.67
-0.57
-0.42
-0.23
-0.03
0.18
0.37
0.54
0.66
0.72
0.73
0.68
0.57
0.41
0.22
0.02
-0.19
-0.38
-0.54
-0.66
-0.72
-0.72
-0.66
-0.54
-0.38
-0.19
0.02
0.22
0.41
0.57
0.68
0.73
0.73
0.67
0.56
0.4
0.21
0
-0.21
-0.4
-0.56
-0.67
-0.72
-0.72
-0.66
-0.54
-0.38
-0.19
0.01
0.22
0.41
0.56
0.68
0.74
0.75
0.74
0.68
0.58
0.42
0.23
0.03
-0.18
-0.37
-0.53
-0.65
-0.72
-0.73
-0.68
-0.58
-0.43
-0.24
-0.04
0.17
0.36
0.53
0.65
0.73
0.75
0.72
0.62
0.49
0.31
0.11
-0.1
-0.3
-0.48
-0.61
-0.7
-0.74
-0.71
-0.63
-0.5
-0.33
-0.14
0.07
0.27
0.45
0.6
0.7
0.75
0.75
0.69
0.58
0.42
0.24
0.03
-0.18
-0.37
-0.53
-0.65
-0.72
-0.74
-0.7
-0.6
-0.46
-0.28
-0.08
0.12
0.32
0.5
0.63
0.67
0.73
0.77
0.75
0.68
0.56
0.4
0.21
0
-0.2
-0.39
-0.55
-0.67
-0.73
-0.74
-0.69
-0.6
-0.45
-0.27
-0.07
0.13
0.33
0.5
0.64
0.73
0.77
0.76
0.69
0.57
0.41
0.22
0.02
-0.18
-0.37
-0.54
-0.66
-0.73
-0.75
-0.71
-0.62
-0.48
-0.31
-0.11
0.09
0.29
0.47
0.62
0.72
0.77
0.77
0.72
0.61
0.47
0.29
0.09
-0.12
-0.31
-0.49
-0.62
-0.71
-0.75
-0.74
-0.66
-0.54
-0.39
-0.2
0
0.01
0.21
0.39
0.56
0.68
0.76
0.78
0.76
0.68
0.55
0.39
0.2
0
-0.21
-0.39
-0.55
-0.67
-0.74
-0.76
-0.72
-0.63
-0.5
-0.33
-0.13
0.07
0.27
0.45
0.6
0.71
0.78
0.79
0.75
0.66
0.52
0.35
0.16
-0.05
-0.25
-0.43
-0.58
-0.69
-0.75
-0.76
-0.72
-0.62
-0.48
-0.31
-0.12
0.09
0.28
0.46
0.61
0.72
0.78
0.79
0.75
0.66
0.52
0.35
0.16
-0.04
-0.24
-0.42
-0.58
-0.69
-0.76
-0.77
-0.77
-0.73
-0.64
-0.51
-0.34
-0.15
0.05
0.25
0.44
0.59
0.71
0.78
0.8
0.77
0.69
0.56
0.4
0.21
0.01
-0.19
-0.38
-0.54
-0.67
-0.75
-0.78
-0.75
-0.68
-0.56
-0.41
-0.22
-0.02
0.18
0.37
0.53
0.66
0.75
0.8
0.79
0.73
0.63
0.48
0.31
0.11
-0.09
-0.29
-0.46
-0.61
-0.72
-0.77
-0.78
-0.74
-0.64
-0.51
-0.34
-0.15
0.05
0.25
0.43
0.58
0.7
0.78
0.8
0.78
0.71
0.59
0.44
0.36
0.25
0.06
-0.14
-0.34
-0.51
-0.64
-0.74
-0.79
-0.78
-0.73
-0.63
-0.49
-0.31
-0.12
0.08
0.27
0.45
0.6
0.72
0.79
0.81
0.78
0.7
0.58
0.43
0.25
0.05
-0.15
-0.34
-0.51
-0.65
-0.74
-0.79
-0.79
-0.74
-0.64
-0.5
-0.33
-0.14
0.06
0.25
0.43
0.59
0.71
0.78
0.81
0.79
0.72
0.61
0.46
0.29
0.09
-0.11
-0.3
-0.48
-0.62
-0.73
-0.79
-0.8
-0.76
-0.68
-0.55
-0.39
-0.2
-0.01
0.19
0.21
0.38
0.54
0.67
0.76
0.81
0.81
0.76
0.67
0.53
0.37
0.18
-0.02
-0.22
-0.4
-0.56
-0.69
-0.77
-0.81
-0.79
-0.73
-0.63
-0.48
-0.31
-0.12
0.08
0.27
0.45
0.6
0.72
0.79
0.82
0.8
0.74
0.63
0.48
0.31
0.12
-0.08
-0.28
-0.45
-0.6
-0.72
-0.79
-0.81
-0.79
-0.71
-0.6
-0.45
-0.27
-0.08
0.12
0.31
0.48
0.63
0.74
0.81
0.83
0.8
0.73
0.62
0.47
0.29
0.1
-0.1
-0.29
-0.43
-0.47
-0.61
-0.73
-0.79
-0.82
-0.79
-0.72
-0.6
-0.45
-0.28
-0.08
0.11
0.3
0.48
0.63
0.74
0.81
0.83
0.81
0.75
0.64
0.49
0.32
0.13
-0.07
-0.26
-0.44
-0.59
-0.71
-0.79
-0.82
-0.8
-0.74
-0.63
-0.49
-0.32
-0.13
0.06
0.26
0.43
0.59
0.71
0.8
0.84
0.83
0.78
0.69
0.56
0.39
0.21
0.01
-0.18
-0.37
-0.53
-0.67
-0.76
-0.81
-0.82
-0.78
-0.69
-0.57
-0.41
-0.23
-0.04
0.16
0.35
0.43
0.52
0.66
0.76
0.83
0.85
0.82
0.75
0.64
0.5
0.33
0.14
-0.06
-0.25
-0.43
-0.58
-0.71
-0.79
-0.82
-0.81
-0.76
-0.66
-0.53
-0.36
-0.18
0.02
0.21
0.39
0.56
0.69
0.79
0.84
0.85
0.82
0.74
0.63
0.48
0.3
0.11
-0.08
-0.28
-0.45
-0.6
-0.72
-0.79
-0.83
-0.81
-0.76
-0.66
-0.52
-0.35
-0.17
0.02
0.22
0.4
0.56
0.69
0.79
0.85
0.86
0.83
0.75
0.64
0.49
0.32
0.13
-0.06
-0.1
-0.26
-0.43
-0.59
-0.71
-0.79
-0.83
-0.82
-0.77
-0.68
-0.55
-0.39
-0.21
-0.02
0.18
0.36
0.53
0.67
0.77
0.84
0.87
0.85
0.78
0.68
0.54
0.38
0.2
0
-0.19
-0.37
-0.54
-0.67
-0.77
-0.82
-0.84
-0.8
-0.73
-0.61
-0.46
-0.29
-0.11
0.09
0.28
0.45
0.61
0.73
0.82
0.86
0.87
0.83
0.74
0.63
0.47
0.3
0.11
-0.08
-0.27
-0.45
-0.6
-0.72
-0.8
-0.84
-0.83
-0.79
-0.7
-0.57
Dataseries Y:
Download CSV
Histogram 
0
-0.21
-0.4
-0.55
-0.66
-0.72
-0.71
-0.64
-0.52
-0.35
-0.15
0.06
0.26
0.45
0.59
0.69
0.73
0.71
0.63
0.5
0.33
0.13
-0.08
-0.28
-0.46
-0.6
-0.69
-0.73
-0.7
-0.61
-0.48
-0.3
-0.1
0.11
0.32
0.49
0.63
0.71
0.74
0.7
0.61
0.48
0.3
0.1
-0.11
-0.31
-0.48
-0.62
-0.7
-0.73
-0.7
-0.61
-0.47
-0.29
-0.09
0.12
0.32
0.5
0.63
0.72
0.74
0.71
0.62
0.49
0.31
0.12
-0.03
-0.09
-0.29
-0.47
-0.61
-0.7
-0.73
-0.71
-0.63
-0.5
-0.33
-0.13
0.08
0.28
0.47
0.61
0.71
0.75
0.73
0.66
5.37E-001
0.37
0.18
-0.02
-0.23
-0.41
-0.57
-0.68
-0.73
-0.73
-0.67
-0.56
-0.4
-0.22
-0.01
0.2
0.39
0.55
0.67
0.74
0.75
0.71
0.61
0.47
0.29
0.09
-0.12
-0.31
-0.48
-0.62
-0.71
-0.74
-0.72
-0.64
-0.51
-0.34
-0.15
0.06
0.26
0.45
0.6
0.7
0.75
0.75
0.69
0.58
0.43
0.37
0.25
0.05
-0.16
-0.35
-0.52
-0.64
-0.72
-0.75
-0.72
-0.63
-0.5
-0.33
-0.13
0.08
0.28
0.46
0.61
0.71
0.76
0.76
0.7
0.59
0.43
0.25
0.05
-0.15
-0.35
-0.51
-0.64
-0.73
-0.76
-0.73
-0.65
-0.52
-0.35
-0.16
0.05
0.25
0.43
0.59
0.7
0.76
0.76
0.72
0.62
0.48
0.3
0.1
-0.1
-0.3
-0.47
-0.61
-0.71
-0.76
-0.75
-0.69
-0.58
-0.42
-0.24
-0.03
0.17
0.36
0.53
0.66
0.74
0.77
0.77
0.75
0.67
0.55
0.39
0.2
-9.32E-004
-0.2
-0.39
-0.55
-0.67
-0.75
-0.77
-0.74
-0.65
-0.52
-0.35
-0.16
0.04
0.24
0.43
0.58
0.7
0.76
0.77
0.73
0.64
0.51
0.34
0.15
-0.06
-0.25
-0.44
-0.59
-0.7
-0.76
-0.77
-0.73
-0.64
-0.5
-0.33
-0.13
0.07
0.27
0.45
0.6
0.71
0.77
0.78
0.73
0.64
0.5
0.33
0.14
-0.06
-0.26
-0.44
-0.59
-0.7
-0.77
-0.78
-0.74
-0.65
-0.51
-0.34
-0.15
-9.14E-004
0.05
0.25
0.43
0.59
0.7
0.77
0.78
0.74
0.66
0.53
0.37
0.18
-0.02
-0.22
-0.41
-0.56
-0.69
-0.76
-0.79
-0.76
-0.69
-0.56
-0.4
-0.22
-0.01
0.19
0.38
0.54
0.67
0.75
0.79
0.77
0.7
0.59
0.44
0.26
0.07
-0.13
-0.33
-0.5
-0.64
-0.74
-0.79
-0.79
-0.74
-0.64
-0.5
-0.32
-0.13
0.07
0.27
0.45
0.6
0.71
0.78
0.79
0.76
0.67
0.55
0.38
0.2
0
-0.2
-0.38
-0.55
-0.67
-0.72
-0.76
-0.8
-0.78
-0.72
-0.61
-0.46
-0.28
-0.09
0.12
0.31
0.48
0.63
0.73
0.79
0.8
0.75
0.66
0.53
0.37
0.18
-0.01
-0.21
-0.4
-0.56
-0.68
-0.77
-0.8
-0.79
-0.72
-0.61
-0.46
-0.29
-0.09
0.11
0.31
0.48
0.63
0.73
0.79
0.8
0.77
0.68
0.56
0.4
0.21
0.02
-0.18
-0.37
-0.53
-0.66
-0.76
-0.8
-0.8
-0.74
-0.64
-0.5
-0.33
-0.14
0.06
0.26
0.44
0.59
0.71
0.78
0.81
0.79
0.79
0.72
0.61
0.46
0.29
0.1
-0.1
-0.29
-0.47
-0.61
-0.72
-0.79
-0.81
-0.78
-0.7
-0.58
-0.42
-0.24
-0.04
0.16
0.35
0.52
0.66
0.76
0.81
0.82
0.77
0.69
0.56
0.4
0.22
0.02
-0.17
-0.36
-0.52
-0.66
-0.75
-0.8
-0.81
-0.76
-0.67
-0.54
-0.37
-0.19
0.01
0.21
0.4
0.56
0.69
0.78
0.82
0.82
0.77
0.67
0.54
0.38
0.19
0
-0.2
-0.38
-0.54
-0.67
-0.76
-0.81
-0.81
-0.76
-0.69
-0.66
-0.53
-0.37
-0.18
0.02
0.22
0.4
0.56
0.69
0.78
0.83
0.83
0.78
0.68
0.55
0.39
0.21
0.02
-0.18
-0.36
-0.53
-0.66
-0.76
-0.81
-0.82
-0.77
-0.69
-0.56
-0.4
-0.22
-0.02
0.18
0.36
0.53
0.67
0.77
0.83
0.84
0.8
0.72
0.6
0.45
0.27
0.08
-0.11
-0.3
-0.48
-0.62
-0.73
-0.8
-0.82
-0.8
-0.73
-0.62
-0.47
-0.3
-0.11
0.09
0.28
0.46
0.61
0.73
0.81
0.84
0.83
0.77
0.73
0.67
0.54
0.37
0.19
0
-0.2
-0.38
-0.54
-0.67
-0.77
-0.82
-0.83
-0.79
-0.7
-0.58
-0.42
-0.24
-0.05
0.15
0.34
0.51
0.65
0.76
0.83
0.85
0.82
0.76
0.65
0.51
0.34
0.15
-0.04
-0.24
-0.41
-0.57
-0.7
-0.78
-0.83
-0.83
-0.78
-0.69
-0.57
-0.41
-0.23
-0.03
0.16
0.35
0.52
0.66
0.77
0.83
0.85
0.83
0.76
0.65
0.51
0.34
0.16
-0.04
-0.23
-0.41
-0.57
-0.69
-0.78
-0.83
-0.84
-0.83
-0.8
-0.71
-0.59
-0.44
-0.26
-0.07
0.13
0.32
0.49
0.64
0.75
0.82
0.85
0.84
0.78
0.68
0.55
0.39
0.21
0.01
-0.18
-0.36
-0.53
-0.66
-0.77
-8.29E-001
-0.85
-0.82
-0.75
-0.64
-0.5
-0.33
-0.14
0.05
0.24
0.42
0.58
0.71
0.8
0.85
0.85
0.81
0.73
0.62
0.47
0.3
0.11
-0.08
-0.27
-0.45
-0.6
-0.72
-0.81
-0.85
-0.85
-8.05E-001
-0.72
-0.6
-0.44
-0.26
-0.07
0.12
0.31
0.48
0.63
Dataseries Z:
Download CSV
Histogram 
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.03
0.03
0.02
0
-0.01
-0.02
-0.03
-0.03
-0.04
-0.04
-0.04
-0.03
-0.03
-0.02
0
0.01
0.02
0.03
0.03
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-2.93E-002
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.04
-0.03
-2.78E-002
-0.03
-0.02
0
0.01
0.02
0.03
0.03
0.04
0.04
0.04
0.03
0.03
0.02
0
-0.01
-0.02
-0.03
-0.03
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
8.15E-003
0.02
0.03
0.04
0.04
0.04
0.04
0.03
0.03
0.02
6.47E-003
0
-0.02
-0.03
-0.03
-0.04
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
3.63E-002
0.04
0.04
0.04
0.03
0.03
0.02
0.01
-0.01
-0.02
-0.03
-0.03
-0.04
-0.04
-4.13E-002
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.04
0.03
1.84E-002
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.04
-0.03
-0.03
-0.02
-0.01
0.01
0.02
0.03
0.03
0.04
0.04
0.04
0.04
0.03
0.02
0.01
9.87E-004
-0.01
-0.02
-0.03
-0.03
-0.04
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
2.48E-004
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.03
0.03
0.02
0.01
0
-0.02
-0.03
-0.03
-0.04
-0.04
-0.04
-0.04
-0.04
-2.81E-002
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.04
0.03
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-4.16E-002
-0.04
-0.03
-0.02
-0.01
2.26E-004
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.04
0.03
0.02
0.01
9.32E-004
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.04
-0.04
-0.03
-0.03
-0.02
-5.01E-003
0.01
0.02
0.03
0.03
0.04
0.04
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.01
0.02
0.03
0.04
0.04
0.04
4.30E-002
0.04
0.03
0.03
0.02
0.01
-0.01
-0.02
-0.03
-0.03
-0.04
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.04
0.03
0.03
0.02
0.01
0
-0.02
-0.03
-0.03
-0.04
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.04
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.03
-4.07E-002
-0.04
-0.05
-0.04
-3.92E-002
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.03
0.04
0.04
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.05
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.04
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.05
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.03
0.04
4.42E-002
0.05
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.03
-0.04
-0.04
-0.05
-0.04
-0.04
-0.04
-0.03
-0.02
-0.01
0
1.36E-002
0.02
0.03
0.04
0.04
0.05
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.04
-0.05
-0.05
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
2.70E-002
0.03
0.04
0.04
0.05
0.04
0.04
0.03
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.05
-0.05
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
4.51E-002
0.05
0.04
0.04
0.03
0.03
0.02
0.01
0
-0.02
-0.03
-0.03
-0.04
-0.04
-0.05
-0.05
-0.04
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.04
0.04
0.05
0.05
0.04
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.05
-0.05
-0.04
-0.04
-0.03
-0.03
-0.02
-0.01
0.01
0.02
0.03
0.03
0.04
0.04
0.05
0.05
0.04
0.04
0.03
0.02
0.01
0
-0.01
-0.02
-0.03
-0.04
-0.04
-0.05
-0.05
-0.05
-0.04
-0.04
-0.03
-0.02
-0.01
-7.05E-004
0.01
0.02
0.03
0.04
0.04
0.05
0.05
0.05
0.04
0.04
0.03
0.02
0.01
0

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time5 seconds
R Server'George Udny Yule' @ 72.249.76.132

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 5 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=11146&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]5 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=11146&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=11146&T=0

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time5 seconds
R Server'George Udny Yule' @ 72.249.76.132


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 50 ; par2 = 50 ; par3 = Y ; par4 = Y ; par5 = X ; par6 = Y ; par7 = Z ;
Parameters (R input):
par1 = 50 ; par2 = 50 ; par3 = Y ; par4 = Y ; par5 = X ; par6 = Y ; par7 = Z ;
R code (references can be found in the software module):
x <- array(x,dim=c(length(x),1))
colnames(x) <- par5
y <- array(y,dim=c(length(y),1))
colnames(y) <- par6
z <- array(z,dim=c(length(z),1))
colnames(z) <- par7
d <- data.frame(cbind(z,y,x))
colnames(d) <- list(par7,par6,par5)
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
if (par1>500) par1 <- 500
if (par2>500) par2 <- 500
if (par1<10) par1 <- 10
if (par2<10) par2 <- 10
library(GenKern)
library(lattice)
panel.hist <- function(x, ...)
{
usr <- par('usr'); on.exit(par(usr))
par(usr = c(usr[1:2], 0, 1.5) )
h <- hist(x, plot = FALSE)
breaks <- h$breaks; nB <- length(breaks)
y <- h$counts; y <- y/max(y)
rect(breaks[-nB], 0, breaks[-1], y, col='black', ...)
}
bitmap(file='cloud1.png')
cloud(z~x*y, screen = list(x=-45, y=45, z=35),xlab=par5,ylab=par6,zlab=par7)
dev.off()
bitmap(file='cloud2.png')
cloud(z~x*y, screen = list(x=35, y=45, z=25),xlab=par5,ylab=par6,zlab=par7)
dev.off()
bitmap(file='cloud3.png')
cloud(z~x*y, screen = list(x=35, y=-25, z=90),xlab=par5,ylab=par6,zlab=par7)
dev.off()
bitmap(file='pairs.png')
pairs(d,diag.panel=panel.hist)
dev.off()
x <- as.vector(x)
y <- as.vector(y)
z <- as.vector(z)
bitmap(file='bidensity1.png')
op <- KernSur(x,y, xgridsize=par1, ygridsize=par2, correlation=cor(x,y), xbandwidth=dpik(x), ybandwidth=dpik(y))
image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main='Bivariate Kernel Density Plot (x,y)',xlab=par5,ylab=par6)
if (par3=='Y') contour(op$xords, op$yords, op$zden, add=TRUE)
if (par4=='Y') points(x,y)
(r<-lm(y ~ x))
abline(r)
box()
dev.off()
bitmap(file='bidensity2.png')
op <- KernSur(y,z, xgridsize=par1, ygridsize=par2, correlation=cor(y,z), xbandwidth=dpik(y), ybandwidth=dpik(z))
op
image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main='Bivariate Kernel Density Plot (y,z)',xlab=par6,ylab=par7)
if (par3=='Y') contour(op$xords, op$yords, op$zden, add=TRUE)
if (par4=='Y') points(y,z)
(r<-lm(z ~ y))
abline(r)
box()
dev.off()
bitmap(file='bidensity3.png')
op <- KernSur(x,z, xgridsize=par1, ygridsize=par2, correlation=cor(x,z), xbandwidth=dpik(x), ybandwidth=dpik(z))
op
image(op$xords, op$yords, op$zden, col=terrain.colors(100), axes=TRUE,main='Bivariate Kernel Density Plot (x,z)',xlab=par5,ylab=par7)
if (par3=='Y') contour(op$xords, op$yords, op$zden, add=TRUE)
if (par4=='Y') points(x,z)
(r<-lm(z ~ x))
abline(r)
box()
dev.off()