## Dump S Version 4 Dump ##
.Last.value
character
character
53
A1
A1inv
I.0
I.1
I.p
change.pt
circ.aov
circ.cor
circ.disp
circ.lin.reg
circ.mean
circ.plot
circ.range
circ.reg
circ.summary
dcard
deg
dmixedvm
dtri
dvm
dwrpcauchy
dwrpnorm
eqscplot
equal.kappa.test
est.kappa
est.rho
kuiper
last.dump
ls.circle
plot.edf
plot.ls.circle
pp.plot
pvm
r.test
rad
rao.homogeneity
rao.spacing
rao.table
rcard
rmixedvm
rose.diag
rtri
rvm
rwrpcauchy
rwrpnorm
rwrpstab
trig.moment
v0.test
vm.bootstrap.ci
vm.ml
watson
watson.two
wrpcauchy.ml
A1
.named_I
structure
2
function
function
2
missing
missing
0
{
{
2
<-
<-
2
name
name
1
result
call
call
3
name
name
1
/
call
call
2
name
name
1
I.1
name
name
1
kappa
call
call
2
name
name
1
I.0
name
name
1
kappa
name
name
1
result
character
character
2
kappa

A1inv
.named_I
structure
2
function
function
2
missing
missing
0
{
{
4
if
if
2
call
call
3
name
name
1
&&
call
call
3
name
name
1
<=
numeric
numeric
1
0
name
name
1
x
call
call
3
name
name
1
<
name
name
1
x
numeric
numeric
1
0.53000000000000003
{
{
1
<-
<-
2
name
name
1
result
call
call
3
name
name
1
+
call
call
3
name
name
1
+
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
x
call
call
3
name
name
1
^
name
name
1
x
numeric
numeric
1
3
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
5
call
call
3
name
name
1
^
name
name
1
x
numeric
numeric
1
5
numeric
numeric
1
6
if
if
2
call
call
3
name
name
1
&&
call
call
3
name
name
1
<=
numeric
numeric
1
0.53000000000000003
name
name
1
x
call
call
3
name
name
1
<
name
name
1
x
numeric
numeric
1
0.84999999999999998
{
{
1
<-
<-
2
name
name
1
result
call
call
3
name
name
1
+
call
call
3
name
name
1
+
numeric
numeric
1
-0.40000000000000002
call
call
3
name
name
1
*
numeric
numeric
1
1.3899999999999999
name
name
1
x
call
call
3
name
name
1
/
numeric
numeric
1
0.42999999999999999
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
1
name
name
1
x
if
if
2
call
call
3
name
name
1
>=
name
name
1
x
numeric
numeric
1
0.84999999999999998
{
{
1
<-
<-
2
name
name
1
result
call
call
3
name
name
1
/
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
^
name
name
1
x
numeric
numeric
1
3
call
call
3
name
name
1
*
numeric
numeric
1
4
call
call
3
name
name
1
^
name
name
1
x
numeric
numeric
1
2
call
call
3
name
name
1
*
numeric
numeric
1
3
name
name
1
x
name
name
1
result
character
character
2
x

I.0
.named_I
structure
2
function
function
2
missing
missing
0
{
{
20
<-
<-
2
name
name
1
p1
numeric
numeric
1
1
<-
<-
2
name
name
1
p2
numeric
numeric
1
3.5156228999999999
<-
<-
2
name
name
1
p3
numeric
numeric
1
3.0899424
<-
<-
2
name
name
1
p4
numeric
numeric
1
1.2067492
<-
<-
2
name
name
1
p5
numeric
numeric
1
0.26597320000000002
<-
<-
2
name
name
1
p6
call
call
3
name
name
1
/
numeric
numeric
1
0.36076799999999998
numeric
numeric
1
10
<-
<-
2
name
name
1
p7
call
call
3
name
name
1
/
numeric
numeric
1
0.45812999999999998
numeric
numeric
1
100
<-
<-
2
name
name
1
q1
numeric
numeric
1
0.39894227999999998
<-
<-
2
name
name
1
q2
call
call
3
name
name
1
/
numeric
numeric
1
0.13285920000000001
numeric
numeric
1
10
<-
<-
2
name
name
1
q3
call
call
3
name
name
1
/
numeric
numeric
1
0.22531899999999999
numeric
numeric
1
100
<-
<-
2
name
name
1
q4
call
call
3
name
name
1
/
numeric
numeric
1
-0.15756500000000001
numeric
numeric
1
100
<-
<-
2
name
name
1
q5
call
call
3
name
name
1
/
numeric
numeric
1
0.91628100000000001
numeric
numeric
1
100
<-
<-
2
name
name
1
q6
call
call
3
name
name
1
/
numeric
numeric
1
-0.2057706
numeric
numeric
1
10
<-
<-
2
name
name
1
q7
call
call
3
name
name
1
/
numeric
numeric
1
0.2635537
numeric
numeric
1
10
<-
<-
2
name
name
1
q8
call
call
3
name
name
1
/
numeric
numeric
1
-0.1647633
numeric
numeric
1
10
<-
<-
2
name
name
1
q9
call
call
3
name
name
1
/
numeric
numeric
1
0.39237699999999998
numeric
numeric
1
100
<-
<-
2
name
name
1
y
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
/
name
name
1
x
numeric
numeric
1
3.75
numeric
numeric
1
2
<-
<-
2
name
name
1
ax
call
call
2
name
name
1
abs
name
name
1
x
<-
<-
2
name
name
1
z
call
call
3
name
name
1
/
numeric
numeric
1
3.75
name
name
1
ax
call
call
4
name
name
1
ifelse
call
call
3
name
name
1
<
call
call
2
name
name
1
abs
name
name
1
x
numeric
numeric
1
3.75
call
call
3
name
name
1
+
name
name
1
p1
call
call
3
name
name
1
*
name
name
1
y
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
p2
call
call
3
name
name
1
*
name
name
1
y
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
p3
call
call
3
name
name
1
*
name
name
1
y
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
p4
call
call
3
name
name
1
*
name
name
1
y
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
p5
call
call
3
name
name
1
*
name
name
1
y
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
p6
call
call
3
name
name
1
*
name
name
1
y
name
name
1
p7
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
/
call
call
2
name
name
1
exp
name
name
1
ax
call
call
2
name
name
1
sqrt
name
name
1
ax
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
q1
call
call
3
name
name
1
*
name
name
1
z
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
q2
call
call
3
name
name
1
*
name
name
1
z
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
q3
call
call
3
name
name
1
*
name
name
1
z
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
q4
call
call
3
name
name
1
*
name
name
1
z
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
q5
call
call
3
name
name
1
*
name
name
1
z
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
q6
call
call
3
name
name
1
*
name
name
1
z
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
q7
call
call
3
name
name
1
*
name
name
1
z
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
q8
call
call
3
name
name
1
*
name
name
1
z
name
name
1
q9
character
character
2
x

I.1
.named_I
structure
2
function
function
2
missing
missing
0
{
{
4
<-
<-
2
name
name
1
t
call
call
3
name
name
1
/
name
name
1
x
numeric
numeric
1
3.75
if
if
2
call
call
3
name
name
1
<
name
name
1
x
numeric
numeric
1
3.75
{
{
1
<-
<-
2
name
name
1
I1
call
call
3
name
name
1
*
name
name
1
x
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
+
call
call
3
name
name
1
+
call
call
3
name
name
1
+
call
call
3
name
name
1
+
call
call
3
name
name
1
+
numeric
numeric
1
0.5
call
call
3
name
name
1
*
numeric
numeric
1
0.87890594
call
call
3
name
name
1
^
name
name
1
t
numeric
numeric
1
2
call
call
3
name
name
1
*
numeric
numeric
1
0.51498869000000003
call
call
3
name
name
1
^
name
name
1
t
numeric
numeric
1
4
call
call
3
name
name
1
*
numeric
numeric
1
0.15084934
call
call
3
name
name
1
^
name
name
1
t
numeric
numeric
1
6
call
call
3
name
name
1
*
numeric
numeric
1
0.026587329999999999
call
call
3
name
name
1
^
name
name
1
t
numeric
numeric
1
8
call
call
3
name
name
1
*
numeric
numeric
1
0.0030153200000000002
call
call
3
name
name
1
^
name
name
1
t
numeric
numeric
1
10
call
call
3
name
name
1
*
numeric
numeric
1
0.00032411000000000001
call
call
3
name
name
1
^
name
name
1
t
numeric
numeric
1
12
if
if
2
call
call
3
name
name
1
>=
name
name
1
x
numeric
numeric
1
3.75
{
{
1
<-
<-
2
name
name
1
I1
call
call
3
name
name
1
*
call
call
3
name
name
1
*
call
call
3
name
name
1
^
name
name
1
x
(
(
2
name
name
1
(
numeric
numeric
1
-0.5
call
call
2
name
name
1
exp
name
name
1
x
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
-
numeric
numeric
1
0.39894227999999998
call
call
3
name
name
1
*
numeric
numeric
1
0.039880239999999997
call
call
3
name
name
1
^
name
name
1
t
(
(
2
name
name
1
(
numeric
numeric
1
-1
call
call
3
name
name
1
*
numeric
numeric
1
0.0036201800000000002
call
call
3
name
name
1
^
name
name
1
t
(
(
2
name
name
1
(
numeric
numeric
1
-2
call
call
3
name
name
1
*
numeric
numeric
1
0.0016380100000000001
call
call
3
name
name
1
^
name
name
1
t
(
(
2
name
name
1
(
numeric
numeric
1
-3
call
call
3
name
name
1
*
numeric
numeric
1
0.01031555
call
call
3
name
name
1
^
name
name
1
t
(
(
2
name
name
1
(
numeric
numeric
1
-4
call
call
3
name
name
1
*
numeric
numeric
1
0.02282967
call
call
3
name
name
1
^
name
name
1
t
(
(
2
name
name
1
(
numeric
numeric
1
-5
call
call
3
name
name
1
*
numeric
numeric
1
0.028953119999999999
call
call
3
name
name
1
^
name
name
1
t
(
(
2
name
name
1
(
numeric
numeric
1
-6
call
call
3
name
name
1
*
numeric
numeric
1
0.01787654
call
call
3
name
name
1
^
name
name
1
t
(
(
2
name
name
1
(
numeric
numeric
1
-7
call
call
3
name
name
1
*
numeric
numeric
1
0.0042005899999999997
call
call
3
name
name
1
^
name
name
1
t
(
(
2
name
name
1
(
numeric
numeric
1
-8
name
name
1
I1
character
character
2
x

I.p
.named_I
structure
2
function
function
3
missing
missing
0
missing
missing
0
{
{
6
<-
<-
2
name
name
1
I.before
call
call
2
name
name
1
I.0
name
name
1
x
<-
<-
2
name
name
1
I.curr
call
call
2
name
name
1
I.1
name
name
1
x
if
if
2
call
call
3
name
name
1
==
name
name
1
p
numeric
numeric
1
0
<-
<-
2
name
name
1
result
name
name
1
I.before
if
if
2
call
call
3
name
name
1
==
name
name
1
p
numeric
numeric
1
1
<-
<-
2
name
name
1
result
name
name
1
I.curr
if
if
2
call
call
3
name
name
1
&&
call
call
3
name
name
1
!=
name
name
1
p
numeric
numeric
1
0
call
call
3
name
name
1
!=
name
name
1
p
numeric
numeric
1
1
{
{
4
<-
<-
2
name
name
1
n
numeric
numeric
1
1
<-
<-
2
name
name
1
I.next
numeric
numeric
1
0
while
while
2
call
call
3
name
name
1
&&
call
call
3
name
name
1
<
name
name
1
n
name
name
1
p
call
call
3
name
name
1
>=
name
name
1
I.next
numeric
numeric
1
0
{
{
4
<-
<-
2
name
name
1
I.next
call
call
3
name
name
1
-
name
name
1
I.before
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
n
name
name
1
I.curr
name
name
1
x
<-
<-
2
name
name
1
I.before
name
name
1
I.curr
<-
<-
2
name
name
1
I.curr
name
name
1
I.next
<-
<-
2
name
name
1
n
call
call
3
name
name
1
+
name
name
1
n
numeric
numeric
1
1
if
if
3
call
call
3
name
name
1
<
name
name
1
I.next
numeric
numeric
1
0
<-
<-
2
name
name
1
result
numeric
numeric
1
0
<-
<-
2
name
name
1
result
name
name
1
I.next
name
name
1
result
character
character
3
p
x

change.pt
.named_I
structure
2
function
function
2
missing
missing
0
{
{
16
<-
<-
2
name
name
1
phi
.named_I
structure
2
function
function
2
missing
missing
0
{
{
3
<-
<-
2
name
name
1
arg
call
call
2
name
name
1
A1inv
name
name
1
x
if
if
3
call
call
3
name
name
1
!=
call
call
2
name
name
1
I.0
name
name
1
arg
name
name
1
Inf
<-
<-
2
name
name
1
result
call
call
3
name
name
1
-
call
call
3
name
name
1
*
name
name
1
x
call
call
2
name
name
1
A1inv
name
name
1
x
call
call
2
name
name
1
log
call
call
2
name
name
1
I.0
name
name
1
arg
<-
<-
2
name
name
1
result
call
call
3
name
name
1
-
call
call
3
name
name
1
*
name
name
1
x
call
call
2
name
name
1
A1inv
name
name
1
x
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
arg
call
call
2
name
name
1
log
call
call
3
name
name
1
*
call
call
3
name
name
1
/
numeric
numeric
1
1
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
name
name
1
arg
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
+
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
/
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
8
name
name
1
arg
call
call
3
name
name
1
/
numeric
numeric
1
9
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
128
call
call
3
name
name
1
^
name
name
1
arg
numeric
numeric
1
2
call
call
3
name
name
1
/
numeric
numeric
1
225
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
1024
call
call
3
name
name
1
^
name
name
1
arg
numeric
numeric
1
3
name
name
1
result
character
character
2
x

<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
rho
call
call
2
name
name
1
est.rho
name
name
1
data
<-
<-
2
name
name
1
R1
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
<-
<-
2
name
name
1
R2
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
<-
<-
2
name
name
1
V
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
numeric
numeric
1
1
{
{
3
<-
<-
2
call
call
3
name
name
1
[
name
name
1
R1
name
name
1
k
call
call
3
name
name
1
*
call
call
2
name
name
1
est.rho
call
call
3
name
name
1
[
name
name
1
data
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
k
name
name
1
k
<-
<-
2
call
call
3
name
name
1
[
name
name
1
R2
name
name
1
k
call
call
3
name
name
1
*
call
call
2
name
name
1
est.rho
call
call
3
name
name
1
[
name
name
1
data
call
call
3
name
name
1
:
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
k
numeric
numeric
1
1
name
name
1
n
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
name
name
1
k
if
if
2
call
call
3
name
name
1
&
call
call
3
name
name
1
>=
name
name
1
k
numeric
numeric
1
2
call
call
3
name
name
1
<=
name
name
1
k
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
numeric
numeric
1
2
{
{
1
<-
<-
2
call
call
3
name
name
1
[
name
name
1
V
name
name
1
k
call
call
3
name
name
1
+
call
call
3
name
name
1
*
call
call
3
name
name
1
/
name
name
1
k
name
name
1
n
call
call
2
name
name
1
phi
call
call
3
name
name
1
/
call
call
3
name
name
1
[
name
name
1
R1
name
name
1
k
name
name
1
k
call
call
3
name
name
1
*
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
name
name
1
k
name
name
1
n
call
call
2
name
name
1
phi
call
call
3
name
name
1
/
call
call
3
name
name
1
[
name
name
1
R2
name
name
1
k
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
name
name
1
k
character
character
3
k


<-
<-
2
call
call
3
name
name
1
[
name
name
1
R1
name
name
1
n
call
call
3
name
name
1
*
name
name
1
rho
name
name
1
n
<-
<-
2
call
call
3
name
name
1
[
name
name
1
R2
name
name
1
n
numeric
numeric
1
0
<-
<-
2
name
name
1
R.diff
call
call
3
name
name
1
-
call
call
3
name
name
1
+
name
name
1
R1
name
name
1
R2
call
call
3
name
name
1
*
name
name
1
rho
name
name
1
n
<-
<-
2
name
name
1
rmax
call
call
2
name
name
1
max
name
name
1
R.diff
<-
<-
2
name
name
1
rave
call
call
2
name
name
1
mean
name
name
1
R.diff
<-
<-
2
name
name
1
k.r
call
call
3
name
name
1
[
(
(
2
name
name
1
(
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
call
call
3
name
name
1
==
name
name
1
R.diff
call
call
2
name
name
1
max
name
name
1
R.diff
<-
<-
2
name
name
1
V
call
call
3
name
name
1
[
name
name
1
V
call
call
3
name
name
1
:
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
numeric
numeric
1
2
if
if
3
call
call
3
name
name
1
>
name
name
1
n
numeric
numeric
1
3
{
{
3
<-
<-
2
name
name
1
tmax
call
call
2
name
name
1
max
name
name
1
V
<-
<-
2
name
name
1
tave
call
call
2
name
name
1
mean
name
name
1
V
<-
<-
2
name
name
1
k.t
call
call
3
name
name
1
+
call
call
3
name
name
1
[
(
(
2
name
name
1
(
call
call
3
name
name
1
:
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
numeric
numeric
1
3
call
call
3
name
name
1
==
name
name
1
V
call
call
2
name
name
1
max
name
name
1
V
numeric
numeric
1
1
call
call
2
name
name
1
stop
character
character
1
Sample size must be at least 4
call
call
9
name
name
1
data.frame
name
name
1
n
name
name
1
rho
name
name
1
rmax
name
name
1
k.r
name
name
1
rave
name
name
1
tmax
name
name
1
k.t
name
name
1
tave
character
character
2
data

circ.aov
.named_I
structure
2
function
function
7
missing
missing
0
missing
missing
0
name
name
1
NULL
character
character
1
F.test
name
name
1
T
name
name
1
T
{
{
11
<-
<-
2
name
name
1
ns
.named_I
structure
2
call
call
4
name
name
1
tapply
name
name
1
theta
name
name
1
group
name
name
1
length
character
character
4



FUN
<-
<-
2
name
name
1
resultant
.named_I
structure
2
call
call
4
name
name
1
tapply
name
name
1
theta
name
name
1
group
.named_I
structure
2
function
function
2
missing
missing
0
call
call
3
name
name
1
*
call
call
2
name
name
1
est.rho
name
name
1
x
call
call
2
name
name
1
length
name
name
1
x
character
character
2
x

character
character
4



FUN
<-
<-
2
name
name
1
mean.dirs
.named_I
structure
2
call
call
4
name
name
1
tapply
name
name
1
theta
name
name
1
group
name
name
1
circ.mean
character
character
4



FUN
<-
<-
2
name
name
1
kappas
.named_I
structure
2
call
call
4
name
name
1
tapply
name
name
1
theta
name
name
1
group
name
name
1
est.kappa
character
character
4



FUN
<-
<-
2
name
name
1
grps
call
call
2
name
name
1
length
name
name
1
resultant
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
group
<-
<-
2
name
name
1
res.all
call
call
3
name
name
1
*
call
call
2
name
name
1
est.rho
name
name
1
theta
name
name
1
n
<-
<-
2
name
name
1
mean.dir.all
call
call
2
name
name
1
circ.mean
name
name
1
theta
<-
<-
2
name
name
1
kappa.all
call
call
2
name
name
1
est.kappa
name
name
1
theta
if
if
3
call
call
3
name
name
1
==
name
name
1
method
character
character
1
F.test
{
{
12
if
if
2
call
call
2
name
name
1
!
call
call
2
name
name
1
is.null
name
name
1
kappa
call
call
2
name
name
1
warning
character
character
1
Specified value of kappa is not used in the F-test
<-
<-
2
name
name
1
sum.res
call
call
2
name
name
1
sum
name
name
1
resultant
<-
<-
2
name
name
1
df
call
call
4
name
name
1
c
call
call
3
name
name
1
-
name
name
1
grps
integer
integer
1
1
call
call
3
name
name
1
-
name
name
1
n
name
name
1
grps
call
call
3
name
name
1
-
name
name
1
n
integer
integer
1
1
<-
<-
2
name
name
1
SS
call
call
4
name
name
1
c
call
call
3
name
name
1
-
name
name
1
sum.res
name
name
1
res.all
call
call
3
name
name
1
-
name
name
1
n
name
name
1
sum.res
call
call
3
name
name
1
-
name
name
1
n
name
name
1
res.all
<-
<-
2
name
name
1
MS
call
call
3
name
name
1
/
name
name
1
SS
name
name
1
df
if
if
2
call
call
3
name
name
1
==
name
name
1
F.mod
name
name
1
T
<-
<-
2
name
name
1
F.stat
call
call
3
name
name
1
/
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
integer
integer
1
1
call
call
3
name
name
1
/
integer
integer
1
3
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
8
call
call
2
name
name
1
est.kappa
name
name
1
theta
call
call
3
name
name
1
[
name
name
1
MS
integer
integer
1
1
call
call
3
name
name
1
[
name
name
1
MS
integer
integer
1
2
if
if
2
call
call
3
name
name
1
==
name
name
1
F.mod
name
name
1
F
<-
<-
2
name
name
1
F.stat
call
call
3
name
name
1
/
call
call
3
name
name
1
[
name
name
1
MS
integer
integer
1
1
call
call
3
name
name
1
[
name
name
1
MS
integer
integer
1
2
<-
<-
2
name
name
1
p.value
call
call
3
name
name
1
-
integer
integer
1
1
call
call
4
name
name
1
pf
name
name
1
F.stat
call
call
3
name
name
1
-
name
name
1
grps
integer
integer
1
1
call
call
3
name
name
1
-
name
name
1
n
name
name
1
grps
<-
<-
2
name
name
1
result.matrix
.named_I
structure
2
call
call
6
name
name
1
cbind
name
name
1
df
name
name
1
SS
name
name
1
MS
call
call
4
name
name
1
c
name
name
1
F.stat
name
name
1
NA
name
name
1
NA
call
call
4
name
name
1
c
name
name
1
p.value
name
name
1
NA
name
name
1
NA
character
character
6




F
p
<-
<-
2
call
call
2
name
name
1
dimnames
name
name
1
result.matrix
call
call
3
name
name
1
list
call
call
4
name
name
1
c
character
character
1
Between
character
character
1
Within
character
character
1
Total
call
call
6
name
name
1
c
character
character
1
df
character
character
1
SS
character
character
1
MS
character
character
1
F
character
character
1
p
<-
<-
2
name
name
1
result.list
.named_I
structure
2
call
call
12
name
name
1
list
name
name
1
mean.dirs
name
name
1
mean.dir.all
name
name
1
kappas
name
name
1
kappa.all
name
name
1
resultant
name
name
1
res.all
name
name
1
df
name
name
1
SS
name
name
1
MS
name
name
1
F.stat
name
name
1
p.value
character
character
12

mu
mu.all
kappa
kappa.all
r
r.all
df
SS
MS
F.stat
p.value
if
if
2
call
call
3
name
name
1
==
name
name
1
print
name
name
1
T
{
{
3
call
call
5
name
name
1
cat
character
character
1
\n
character
character
1
Circular Analysis of Variance: High Concentration F-Test
character
character
1
\n
character
character
1
\n
call
call
2
name
name
1
print
name
name
1
result.matrix
call
call
3
name
name
1
cat
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
==
name
name
1
method
character
character
1
LRT
{
{
8
if
if
2
call
call
2
name
name
1
is.null
name
name
1
kappa
<-
<-
2
name
name
1
kappa
name
name
1
kappa.all
<-
<-
2
name
name
1
stat1
call
call
3
name
name
1
-
integer
integer
1
1
call
call
3
name
name
1
*
call
call
3
name
name
1
*
call
call
3
name
name
1
/
integer
integer
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
4
name
name
1
kappa
call
call
2
name
name
1
A1
name
name
1
kappa
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
integer
integer
1
1
name
name
1
ns
call
call
3
name
name
1
/
integer
integer
1
1
name
name
1
n
<-
<-
2
name
name
1
stat2
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
kappa
call
call
2
name
name
1
sum
call
call
3
name
name
1
*
name
name
1
resultant
(
(
2
name
name
1
(
call
call
3
name
name
1
-
integer
integer
1
1
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
mean.dirs
name
name
1
mean.dir.all
<-
<-
2
name
name
1
chisq.stat
call
call
3
name
name
1
*
name
name
1
stat1
name
name
1
stat2
<-
<-
2
name
name
1
p.value
call
call
3
name
name
1
-
integer
integer
1
1
call
call
3
name
name
1
pchisq
name
name
1
chisq.stat
call
call
3
name
name
1
-
name
name
1
grps
integer
integer
1
1
<-
<-
2
name
name
1
result.vector
.named_I
structure
2
call
call
4
name
name
1
c
call
call
3
name
name
1
-
name
name
1
grps
integer
integer
1
1
name
name
1
chisq.stat
name
name
1
p.value
character
character
4

df
Chi.Sq
p
<-
<-
2
name
name
1
result.list
.named_I
structure
2
call
call
10
name
name
1
list
name
name
1
mean.dirs
name
name
1
mean.dir.all
name
name
1
kappas
name
name
1
kappa.all
name
name
1
resultant
name
name
1
res.all
call
call
3
name
name
1
-
name
name
1
grps
integer
integer
1
1
name
name
1
chisq.stat
name
name
1
p.value
character
character
10

mu
mu.all
kappa
kappa.all
r
r.all
df
chisq.stat
p.value
if
if
2
call
call
3
name
name
1
==
name
name
1
print
name
name
1
T
{
{
3
call
call
5
name
name
1
cat
character
character
1
\n
character
character
1
Circular Analysis of Variance: Likelihood Ratio Test
character
character
1
\n
character
character
1
\n
call
call
2
name
name
1
print
name
name
1
result.vector
call
call
3
name
name
1
cat
character
character
1
\n
character
character
1
\n
call
call
2
name
name
1
stop
character
character
1
Invalid method specified
name
name
1
result.list
character
character
7
group
theta
kappa
method
F.mod
print

circ.cor
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
name
name
1
F
{
{
9
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
alpha
<-
<-
2
name
name
1
alpha.bar
call
call
2
name
name
1
circ.mean
name
name
1
alpha
<-
<-
2
name
name
1
beta.bar
call
call
2
name
name
1
circ.mean
name
name
1
beta
<-
<-
2
name
name
1
num
call
call
2
name
name
1
sum
call
call
3
name
name
1
*
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
alpha
name
name
1
alpha.bar
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
beta
name
name
1
beta.bar
<-
<-
2
name
name
1
den
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
*
call
call
2
name
name
1
sum
call
call
3
name
name
1
^
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
alpha
name
name
1
alpha.bar
numeric
numeric
1
2
call
call
2
name
name
1
sum
call
call
3
name
name
1
^
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
beta
name
name
1
beta.bar
numeric
numeric
1
2
<-
<-
2
name
name
1
r
call
call
3
name
name
1
/
name
name
1
num
name
name
1
den
<-
<-
2
name
name
1
result
call
call
2
name
name
1
data.frame
name
name
1
r
if
if
2
call
call
3
name
name
1
==
name
name
1
test
name
name
1
T
{
{
6
<-
<-
2
name
name
1
l20
call
call
2
name
name
1
mean
call
call
3
name
name
1
^
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
alpha
name
name
1
alpha.bar
numeric
numeric
1
2
<-
<-
2
name
name
1
l02
call
call
2
name
name
1
mean
call
call
3
name
name
1
^
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
beta
name
name
1
beta.bar
numeric
numeric
1
2
<-
<-
2
name
name
1
l22
call
call
2
name
name
1
mean
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
^
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
alpha
name
name
1
alpha.bar
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
^
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
beta
name
name
1
beta.bar
numeric
numeric
1
2
<-
<-
2
name
name
1
test.stat
call
call
3
name
name
1
*
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
name
name
1
n
name
name
1
l20
name
name
1
l02
name
name
1
l22
name
name
1
r
<-
<-
2
name
name
1
p.value
call
call
3
name
name
1
*
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
2
name
name
1
pnorm
call
call
2
name
name
1
abs
name
name
1
test.stat
<-
<-
2
name
name
1
result
call
call
4
name
name
1
data.frame
name
name
1
r
name
name
1
test.stat
name
name
1
p.value
name
name
1
result
character
character
4
alpha
beta
test

circ.disp
.named_I
structure
2
function
function
2
missing
missing
0
{
{
7
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
c
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
name
name
1
data
<-
<-
2
name
name
1
s
call
call
2
name
name
1
sum
call
call
2
name
name
1
sin
name
name
1
data
<-
<-
2
name
name
1
r
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
c
numeric
numeric
1
2
call
call
3
name
name
1
^
name
name
1
s
numeric
numeric
1
2
<-
<-
2
name
name
1
rbar
call
call
3
name
name
1
/
name
name
1
r
name
name
1
n
<-
<-
2
name
name
1
var
call
call
3
name
name
1
-
numeric
numeric
1
1
name
name
1
rbar
call
call
5
name
name
1
data.frame
name
name
1
n
name
name
1
r
name
name
1
rbar
name
name
1
var
character
character
2
data

circ.lin.reg
.named_I
structure
2
function
function
7
missing
missing
0
missing
missing
0
missing
missing
0
name
name
1
F
name
name
1
T
numeric
numeric
1
1e-010
{
{
25
if
if
2
call
call
2
name
name
1
is.vector
name
name
1
x
<-
<-
2
name
name
1
x
call
call
2
name
name
1
cbind
name
name
1
x
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
theta
<-
<-
2
name
name
1
betaPrev
name
name
1
beta0
<-
<-
2
name
name
1
S
call
call
3
name
name
1
/
call
call
2
name
name
1
sum
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
theta
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
atan
call
call
3
name
name
1
%*%
name
name
1
x
name
name
1
betaPrev
name
name
1
n
<-
<-
2
name
name
1
C
call
call
3
name
name
1
/
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
theta
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
atan
call
call
3
name
name
1
%*%
name
name
1
x
name
name
1
betaPrev
name
name
1
n
<-
<-
2
name
name
1
R
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
S
integer
integer
1
2
call
call
3
name
name
1
^
name
name
1
C
integer
integer
1
2
<-
<-
2
name
name
1
mu
call
call
3
name
name
1
atan
name
name
1
S
name
name
1
C
<-
<-
2
name
name
1
k
call
call
2
name
name
1
A1inv
name
name
1
R
<-
<-
2
name
name
1
diff
call
call
3
name
name
1
+
name
name
1
tol
integer
integer
1
1
<-
<-
2
name
name
1
iter
integer
integer
1
0
while
while
2
call
call
3
name
name
1
>
name
name
1
diff
name
name
1
tol
{
{
14
<-
<-
2
name
name
1
iter
call
call
3
name
name
1
+
name
name
1
iter
integer
integer
1
1
<-
<-
2
name
name
1
u
call
call
3
name
name
1
*
name
name
1
k
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
atan
call
call
3
name
name
1
%*%
name
name
1
x
name
name
1
betaPrev
<-
<-
2
name
name
1
A
.named_I
structure
2
call
call
3
name
name
1
diag
call
call
3
name
name
1
*
name
name
1
k
call
call
2
name
name
1
A1
name
name
1
k
name
name
1
n
character
character
3


nrow
<-
<-
2
name
name
1
g.p
.named_I
structure
2
call
call
3
name
name
1
diag
.named_I
structure
2
call
call
5
name
name
1
apply
name
name
1
x
integer
integer
1
1
.named_I
structure
2
function
function
3
missing
missing
0
missing
missing
0
call
call
3
name
name
1
/
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
+
integer
integer
1
1
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
betaPrev
name
name
1
row
integer
integer
1
2
character
character
3
row
betaPrev

name
name
1
betaPrev
character
character
5




betaPrev
name
name
1
n
character
character
3


nrow
<-
<-
2
name
name
1
D
call
call
3
name
name
1
%*%
name
name
1
g.p
name
name
1
x
<-
<-
2
name
name
1
betaNew
call
call
3
name
name
1
$
call
call
2
name
name
1
lm
call
call
3
name
name
1
~
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
D
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
u
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
name
name
1
A
name
name
1
D
name
name
1
betaPrev
call
call
3
name
name
1
-
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
D
name
name
1
A
name
name
1
D
integer
integer
1
1
character
character
1
coefficients
<-
<-
2
name
name
1
diff
call
call
2
name
name
1
abs
call
call
2
name
name
1
max
call
call
3
name
name
1
-
name
name
1
betaNew
name
name
1
betaPrev
<-
<-
2
name
name
1
betaPrev
name
name
1
betaNew
<-
<-
2
name
name
1
S
call
call
3
name
name
1
/
call
call
2
name
name
1
sum
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
name
name
1
theta
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
atan
call
call
3
name
name
1
%*%
name
name
1
x
name
name
1
betaPrev
name
name
1
n
<-
<-
2
name
name
1
C
call
call
3
name
name
1
/
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
theta
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
atan
call
call
3
name
name
1
%*%
name
name
1
x
name
name
1
betaPrev
name
name
1
n
<-
<-
2
name
name
1
R
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
S
integer
integer
1
2
call
call
3
name
name
1
^
name
name
1
C
integer
integer
1
2
<-
<-
2
name
name
1
mu
call
call
3
name
name
1
atan
name
name
1
S
name
name
1
C
<-
<-
2
name
name
1
k
call
call
2
name
name
1
A1inv
name
name
1
R
if
if
2
call
call
3
name
name
1
==
name
name
1
trace
name
name
1
T
{
{
2
<-
<-
2
name
name
1
log.lik
call
call
3
name
name
1
+
call
call
3
name
name
1
*
call
call
2
name
name
1
-
name
name
1
n
call
call
2
name
name
1
log
call
call
2
name
name
1
I.0
name
name
1
k
call
call
3
name
name
1
*
name
name
1
k
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
atan
call
call
3
name
name
1
%*%
name
name
1
x
name
name
1
betaNew
call
call
6
name
name
1
cat
character
character
1
Iteration 
name
name
1
iter
character
character
1
:    Log-Likelihood = 
name
name
1
log.lik
character
character
1
\n
<-
<-
2
name
name
1
log.lik
call
call
3
name
name
1
+
call
call
3
name
name
1
*
call
call
2
name
name
1
-
name
name
1
n
call
call
2
name
name
1
log
call
call
2
name
name
1
I.0
name
name
1
k
call
call
3
name
name
1
*
name
name
1
k
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
atan
call
call
3
name
name
1
%*%
name
name
1
x
name
name
1
betaNew
<-
<-
2
name
name
1
cov.beta
call
call
2
name
name
1
solve
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
D
name
name
1
A
name
name
1
D
<-
<-
2
name
name
1
se.beta
call
call
2
name
name
1
sqrt
call
call
2
name
name
1
diag
name
name
1
cov.beta
<-
<-
2
name
name
1
se.kappa
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
/
integer
integer
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
name
name
1
n
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
-
integer
integer
1
1
call
call
3
name
name
1
^
call
call
2
name
name
1
A1
name
name
1
k
integer
integer
1
2
call
call
3
name
name
1
/
call
call
2
name
name
1
A1
name
name
1
k
name
name
1
k
<-
<-
2
name
name
1
circ.se.mu
call
call
3
name
name
1
/
integer
integer
1
1
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
*
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
call
call
2
name
name
1
ncol
name
name
1
x
name
name
1
k
call
call
2
name
name
1
A1
name
name
1
k
<-
<-
2
name
name
1
z.values
call
call
2
name
name
1
abs
call
call
3
name
name
1
/
name
name
1
betaNew
name
name
1
se.beta
<-
<-
2
name
name
1
p.values
call
call
3
name
name
1
-
integer
integer
1
1
call
call
2
name
name
1
pnorm
name
name
1
z.values
<-
<-
2
name
name
1
result.matrix
.named_I
structure
2
call
call
5
name
name
1
cbind
name
name
1
betaNew
name
name
1
se.beta
name
name
1
z.values
name
name
1
p.values
character
character
5

Coef
SE
Z
p
<-
<-
2
call
call
2
name
name
1
dimnames
name
name
1
result.matrix
call
call
3
name
name
1
list
call
call
3
name
name
1
[[
call
call
2
name
name
1
dimnames
name
name
1
x
integer
integer
1
2
call
call
5
name
name
1
c
character
character
1
Coef
character
character
1
SE
character
character
1
|z|
character
character
1
p
call
call
5
name
name
1
cat
character
character
1
\n
character
character
1
Circular-Linear Regression
character
character
1
\n
character
character
1
\n
call
call
2
name
name
1
print
name
name
1
result.matrix
call
call
3
name
name
1
cat
character
character
1
\n
character
character
1
\n
<-
<-
2
call
call
2
name
name
1
dimnames
name
name
1
betaNew
call
call
3
name
name
1
list
call
call
3
name
name
1
[[
call
call
2
name
name
1
dimnames
name
name
1
x
integer
integer
1
2
call
call
2
name
name
1
c
character
character
1
Estimate
.named_I
structure
2
call
call
9
name
name
1
list
name
name
1
mu
name
name
1
k
name
name
1
betaNew
name
name
1
log.lik
name
name
1
circ.se.mu
name
name
1
se.kappa
name
name
1
cov.beta
name
name
1
se.beta
character
character
9

mu
kappa
beta
log.lik
circ.se.mu
se.kappa
cov.beta
se.beta
character
character
7
x
theta
beta0
trace
print
tol

circ.mean
.named_I
structure
2
function
function
2
missing
missing
0
{
{
4
<-
<-
2
name
name
1
sinr
call
call
2
name
name
1
sum
call
call
2
name
name
1
sin
name
name
1
data
<-
<-
2
name
name
1
cosr
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
name
name
1
data
<-
<-
2
name
name
1
circmean
call
call
3
name
name
1
atan
name
name
1
sinr
name
name
1
cosr
name
name
1
circmean
character
character
2
data

circ.plot
.named_I
structure
2
function
function
17
missing
missing
0
name
name
1
T
name
name
1
F
name
name
1
F
character
character
1
E
character
character
1
counter
character
character
1
deg
numeric
numeric
1
1.5
integer
integer
1
8
character
character
1

integer
integer
1
16
integer
integer
1
0
integer
integer
1
1
integer
integer
1
40
integer
integer
1
1
integer
integer
1
1
{
{
8
<-
<-
2
name
name
1
data
call
call
3
name
name
1
%%
name
name
1
data
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
.named_I
structure
2
call
call
10
name
name
1
eqscplot
call
call
2
name
name
1
cos
.named_I
structure
2
call
call
4
name
name
1
seq
integer
integer
1
0
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
integer
integer
1
1000
character
character
4



length
call
call
2
name
name
1
sin
.named_I
structure
2
call
call
4
name
name
1
seq
integer
integer
1
0
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
integer
integer
1
1000
character
character
4



length
name
name
1
F
character
character
1

character
character
1

name
name
1
main
character
character
1
l
call
call
3
name
name
1
*
name
name
1
shrink
call
call
3
name
name
1
c
integer
integer
1
-1
integer
integer
1
1
call
call
3
name
name
1
*
name
name
1
shrink
call
call
3
name
name
1
c
integer
integer
1
-1
integer
integer
1
1
character
character
10



axes
xlab
ylab
main
type
xlim
ylim
if
if
2
call
call
3
name
name
1
|
call
call
3
name
name
1
==
name
name
1
label
character
character
1
deg
call
call
3
name
name
1
==
name
name
1
label
character
character
1
dir
{
{
10
if
if
2
call
call
3
name
name
1
==
name
name
1
rose
name
name
1
F
.named_I
structure
2
call
call
6
name
name
1
text
integer
integer
1
0
numeric
numeric
1
0.014999999999999999
character
character
1
+
name
name
1
labelex
numeric
numeric
1
0.5
character
character
6




cex
adj
if
if
3
call
call
3
name
name
1
==
name
name
1
label
character
character
1
deg
{
{
4
if
if
2
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
E
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
counter
<-
<-
2
name
name
1
labels
call
call
5
name
name
1
c
character
character
1
0
character
character
1
90
character
character
1
180
character
character
1
270
if
if
2
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
E
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
clock
<-
<-
2
name
name
1
labels
call
call
5
name
name
1
c
character
character
1
0
character
character
1
270
character
character
1
180
character
character
1
90
if
if
2
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
N
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
counter
<-
<-
2
name
name
1
labels
call
call
5
name
name
1
c
character
character
1
270
character
character
1
0
character
character
1
90
character
character
1
180
if
if
2
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
N
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
clock
<-
<-
2
name
name
1
labels
call
call
5
name
name
1
c
character
character
1
90
character
character
1
0
character
character
1
270
character
character
1
180
if
if
2
call
call
3
name
name
1
==
name
name
1
label
character
character
1
dir
<-
<-
2
name
name
1
labels
call
call
5
name
name
1
c
character
character
1
E
character
character
1
N
character
character
1
W
character
character
1
S
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0.90000000000000002
integer
integer
1
1
call
call
3
name
name
1
c
integer
integer
1
0
integer
integer
1
0
.named_I
structure
2
call
call
6
name
name
1
text
numeric
numeric
1
0.875
numeric
numeric
1
0.0074999999999999997
call
call
3
name
name
1
[
name
name
1
labels
integer
integer
1
1
name
name
1
labelex
integer
integer
1
1
character
character
6




cex
adj
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
integer
integer
1
0
numeric
numeric
1
0
call
call
3
name
name
1
c
numeric
numeric
1
0.90000000000000002
integer
integer
1
1
.named_I
structure
2
call
call
6
name
name
1
text
integer
integer
1
0
numeric
numeric
1
0.84999999999999998
call
call
3
name
name
1
[
name
name
1
labels
integer
integer
1
2
name
name
1
labelex
numeric
numeric
1
0.5
character
character
6




cex
adj
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
integer
integer
1
-1
numeric
numeric
1
-0.90000000000000002
call
call
3
name
name
1
c
integer
integer
1
0
integer
integer
1
0
.named_I
structure
2
call
call
6
name
name
1
text
numeric
numeric
1
-0.875
numeric
numeric
1
0.0074999999999999997
call
call
3
name
name
1
[
name
name
1
labels
integer
integer
1
3
name
name
1
labelex
integer
integer
1
0
character
character
6




cex
adj
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
integer
integer
1
0
integer
integer
1
0
call
call
3
name
name
1
c
numeric
numeric
1
-0.90000000000000002
integer
integer
1
-1
.named_I
structure
2
call
call
6
name
name
1
text
integer
integer
1
0
numeric
numeric
1
-0.82499999999999996
call
call
3
name
name
1
[
name
name
1
labels
integer
integer
1
4
name
name
1
labelex
numeric
numeric
1
0.5
character
character
6




cex
adj
if
if
2
call
call
3
name
name
1
>
name
name
1
minor
integer
integer
1
0
{
{
2
<-
<-
2
name
name
1
increment
call
call
3
name
name
1
/
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
minor
integer
integer
1
1
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
minor
{
{
4
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
name
name
1
pi
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
name
name
1
pi
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
name
name
1
pi
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
name
name
1
pi
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
/
call
call
3
name
name
1
*
integer
integer
1
3
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
/
call
call
3
name
name
1
*
integer
integer
1
3
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
/
call
call
3
name
name
1
*
integer
integer
1
3
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
/
call
call
3
name
name
1
*
integer
integer
1
3
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
cos
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
cos
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
sin
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
sin
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
character
character
3
tick


if
if
2
call
call
3
name
name
1
==
name
name
1
pts
name
name
1
T
{
{
2
if
if
3
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
E
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
counter
<-
<-
2
name
name
1
pts.data
name
name
1
data
if
if
3
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
E
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
clock
<-
<-
2
name
name
1
pts.data
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
name
name
1
data
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
if
if
3
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
N
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
counter
<-
<-
2
name
name
1
pts.data
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
data
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
if
if
2
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
N
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
clock
<-
<-
2
name
name
1
pts.data
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
name
name
1
data
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
if
if
3
call
call
3
name
name
1
==
name
name
1
stack
name
name
1
F
{
{
3
<-
<-
2
name
name
1
x
call
call
2
name
name
1
cos
name
name
1
pts.data
<-
<-
2
name
name
1
y
call
call
2
name
name
1
sin
name
name
1
pts.data
.named_I
structure
2
call
call
5
name
name
1
points
name
name
1
x
name
name
1
y
name
name
1
cex
name
name
1
pch
character
character
5



cex
pch
{
{
7
<-
<-
2
name
name
1
bins.count
call
call
2
name
name
1
c
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
bins
<-
<-
2
name
name
1
mids
call
call
2
name
name
1
c
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
bins
<-
<-
2
name
name
1
arc
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
name
name
1
bins
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
bins
{
{
1
if
if
3
call
call
3
name
name
1
==
name
name
1
i
integer
integer
1
1
{
{
2
<-
<-
2
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
call
call
2
name
name
1
sum
call
call
3
name
name
1
|
call
call
3
name
name
1
<=
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
arc
integer
integer
1
2
call
call
3
name
name
1
$
name
name
1
.Machine
character
character
1
single.eps
call
call
3
name
name
1
>
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
call
call
3
name
name
1
/
name
name
1
arc
integer
integer
1
2
call
call
3
name
name
1
$
name
name
1
.Machine
character
character
1
single.eps
<-
<-
2
call
call
3
name
name
1
[
name
name
1
mids
name
name
1
i
integer
integer
1
0
{
{
2
<-
<-
2
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
call
call
2
name
name
1
sum
call
call
3
name
name
1
&
call
call
3
name
name
1
<=
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
call
call
3
name
name
1
/
integer
integer
1
1
integer
integer
1
2
name
name
1
arc
call
call
3
name
name
1
$
name
name
1
.Machine
character
character
1
single.eps
call
call
3
name
name
1
>
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
call
call
3
name
name
1
/
integer
integer
1
1
integer
integer
1
2
name
name
1
arc
call
call
3
name
name
1
$
name
name
1
.Machine
character
character
1
single.eps
<-
<-
2
call
call
3
name
name
1
[
name
name
1
mids
name
name
1
i
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
mids
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
name
name
1
arc
character
character
3
i


if
if
3
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
E
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
clock
<-
<-
2
name
name
1
mids
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
name
name
1
mids
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
if
if
3
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
N
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
counter
<-
<-
2
name
name
1
mids
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
mids
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
if
if
2
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
N
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
clock
<-
<-
2
name
name
1
mids
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
name
name
1
mids
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
<-
<-
2
name
name
1
index
call
call
3
name
name
1
/
name
name
1
cex
name
name
1
dotsep
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
bins
{
{
1
if
if
2
call
call
3
name
name
1
!=
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
numeric
numeric
1
0
{
{
1
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
0
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
numeric
numeric
1
1
{
{
4
<-
<-
2
name
name
1
r
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
*
name
name
1
j
name
name
1
index
<-
<-
2
name
name
1
x
call
call
3
name
name
1
*
name
name
1
r
call
call
2
name
name
1
cos
call
call
3
name
name
1
[
name
name
1
mids
name
name
1
i
<-
<-
2
name
name
1
y
call
call
3
name
name
1
*
name
name
1
r
call
call
2
name
name
1
sin
call
call
3
name
name
1
[
name
name
1
mids
name
name
1
i
.named_I
structure
2
call
call
5
name
name
1
points
name
name
1
x
name
name
1
y
name
name
1
cex
name
name
1
pch
character
character
5



cex
pch
character
character
3
j


character
character
3
i


if
if
2
call
call
3
name
name
1
==
name
name
1
rose
name
name
1
T
{
{
7
<-
<-
2
name
name
1
bins.count
call
call
2
name
name
1
c
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
bins
<-
<-
2
name
name
1
sectors
call
call
2
name
name
1
c
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
bins
<-
<-
2
name
name
1
arc
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
name
name
1
bins
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
bins
{
{
1
if
if
3
call
call
3
name
name
1
==
name
name
1
i
integer
integer
1
1
{
{
2
<-
<-
2
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
call
call
2
name
name
1
sum
call
call
3
name
name
1
|
call
call
3
name
name
1
<=
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
arc
integer
integer
1
2
call
call
3
name
name
1
$
name
name
1
.Machine
character
character
1
single.eps
call
call
3
name
name
1
>
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
call
call
3
name
name
1
/
name
name
1
arc
integer
integer
1
2
call
call
3
name
name
1
$
name
name
1
.Machine
character
character
1
single.eps
<-
<-
2
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
call
call
3
name
name
1
/
name
name
1
arc
integer
integer
1
2
{
{
2
<-
<-
2
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
call
call
2
name
name
1
sum
call
call
3
name
name
1
&
call
call
3
name
name
1
<=
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
call
call
3
name
name
1
/
integer
integer
1
1
integer
integer
1
2
name
name
1
arc
call
call
3
name
name
1
$
name
name
1
.Machine
character
character
1
single.eps
call
call
3
name
name
1
>
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
call
call
3
name
name
1
/
integer
integer
1
1
integer
integer
1
2
name
name
1
arc
call
call
3
name
name
1
$
name
name
1
.Machine
character
character
1
single.eps
<-
<-
2
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
name
name
1
arc
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
character
character
3
i


if
if
3
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
E
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
clock
<-
<-
2
name
name
1
sectors
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
name
name
1
sectors
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
if
if
3
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
N
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
counter
<-
<-
2
name
name
1
sectors
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
sectors
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
if
if
2
call
call
3
name
name
1
&
call
call
3
name
name
1
==
name
name
1
zero
character
character
1
N
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
clock
<-
<-
2
name
name
1
sectors
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
name
name
1
sectors
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
<-
<-
2
name
name
1
radius
call
call
3
name
name
1
*
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
/
name
name
1
bins.count
name
name
1
n
name
name
1
prop
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
bins
{
{
1
if
if
2
call
call
3
name
name
1
!=
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
numeric
numeric
1
0
{
{
1
if
if
3
call
call
3
name
name
1
==
name
name
1
rotation
character
character
1
counter
{
{
3
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
call
call
3
name
name
1
c
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
{
{
3
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
call
call
3
name
name
1
c
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
character
character
3
i


call
call
1
name
name
1
invisible
character
character
17
data
pts
stack
rose
zero
rotation
label
labelex
minor
main
pch
bins
cex
dotsep
shrink
prop

circ.range
.named_I
structure
2
function
function
3
missing
missing
0
name
name
1
F
{
{
7
<-
<-
2
name
name
1
data
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
data
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
spacings
call
call
3
name
name
1
c
call
call
2
name
name
1
diff
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
data
numeric
numeric
1
1
call
call
3
name
name
1
[
name
name
1
data
name
name
1
n
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
range
call
call
3
name
name
1
-
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
2
name
name
1
max
name
name
1
spacings
<-
<-
2
name
name
1
result
call
call
2
name
name
1
data.frame
name
name
1
range
if
if
2
call
call
3
name
name
1
==
name
name
1
test
name
name
1
T
{
{
5
<-
<-
2
name
name
1
stop
call
call
2
name
name
1
floor
call
call
3
name
name
1
/
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
/
name
name
1
range
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
index
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
stop
<-
<-
2
name
name
1
sequence
call
call
3
name
name
1
*
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
^
(
(
2
name
name
1
(
numeric
numeric
1
-1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
index
numeric
numeric
1
1
call
call
3
name
name
1
choose
name
name
1
n
name
name
1
index
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
*
name
name
1
index
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
/
name
name
1
range
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
numeric
numeric
1
1
<-
<-
2
name
name
1
p.value
call
call
2
name
name
1
sum
name
name
1
sequence
<-
<-
2
name
name
1
result
call
call
3
name
name
1
data.frame
name
name
1
range
name
name
1
p.value
name
name
1
result
character
character
3
data
test

circ.reg
.named_I
structure
2
function
function
5
missing
missing
0
missing
missing
0
numeric
numeric
1
1
numeric
numeric
1
0.050000000000000003
{
{
44
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
alpha
<-
<-
2
name
name
1
ctheta
call
call
2
name
name
1
cos
name
name
1
theta
<-
<-
2
name
name
1
stheta
call
call
2
name
name
1
sin
name
name
1
theta
<-
<-
2
name
name
1
order.matrix
call
call
2
name
name
1
t
.named_I
structure
2
call
call
3
name
name
1
matrix
call
call
3
name
name
1
rep
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
order
name
name
1
n
name
name
1
n
character
character
3


ncol
<-
<-
2
name
name
1
cos.alpha
call
call
2
name
name
1
cos
call
call
3
name
name
1
*
name
name
1
alpha
name
name
1
order.matrix
<-
<-
2
name
name
1
sin.alpha
call
call
2
name
name
1
sin
call
call
3
name
name
1
*
name
name
1
alpha
name
name
1
order.matrix
<-
<-
2
name
name
1
cos.lm
call
call
2
name
name
1
lm
call
call
3
name
name
1
~
name
name
1
ctheta
call
call
3
name
name
1
+
name
name
1
cos.alpha
name
name
1
sin.alpha
<-
<-
2
name
name
1
sin.lm
call
call
2
name
name
1
lm
call
call
3
name
name
1
~
name
name
1
stheta
call
call
3
name
name
1
+
name
name
1
cos.alpha
name
name
1
sin.alpha
<-
<-
2
name
name
1
cos.fit
call
call
3
name
name
1
$
name
name
1
cos.lm
character
character
1
fitted
<-
<-
2
name
name
1
sin.fit
call
call
3
name
name
1
$
name
name
1
sin.lm
character
character
1
fitted
<-
<-
2
name
name
1
g1.sq
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
cos.fit
name
name
1
cos.fit
<-
<-
2
name
name
1
g2.sq
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
sin.fit
name
name
1
sin.fit
<-
<-
2
name
name
1
rho
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
g1.sq
name
name
1
g2.sq
name
name
1
n
<-
<-
2
name
name
1
theta.fitted
call
call
3
name
name
1
atan
name
name
1
sin.fit
name
name
1
cos.fit
<-
<-
2
name
name
1
Y1
name
name
1
ctheta
<-
<-
2
name
name
1
Y2
name
name
1
stheta
<-
<-
2
name
name
1
ones
call
call
4
name
name
1
matrix
numeric
numeric
1
1
name
name
1
n
numeric
numeric
1
1
<-
<-
2
name
name
1
X
call
call
4
name
name
1
cbind
name
name
1
ones
name
name
1
cos.alpha
name
name
1
sin.alpha
<-
<-
2
name
name
1
W
call
call
3
name
name
1
cbind
call
call
2
name
name
1
cos
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
order
numeric
numeric
1
1
name
name
1
alpha
call
call
2
name
name
1
sin
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
order
numeric
numeric
1
1
name
name
1
alpha
<-
<-
2
name
name
1
M
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
name
name
1
X
call
call
2
name
name
1
solve
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
X
name
name
1
X
call
call
2
name
name
1
t
name
name
1
X
<-
<-
2
name
name
1
I
call
call
2
name
name
1
diag
name
name
1
n
<-
<-
2
name
name
1
H
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
W
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
I
name
name
1
M
name
name
1
W
<-
<-
2
name
name
1
N
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
name
name
1
W
call
call
2
name
name
1
solve
name
name
1
H
call
call
2
name
name
1
t
name
name
1
W
<-
<-
2
name
name
1
c
call
call
3
name
name
1
-
name
name
1
n
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
order
numeric
numeric
1
1
<-
<-
2
name
name
1
N1
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
Y1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
I
name
name
1
M
name
name
1
N
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
I
name
name
1
M
name
name
1
Y1
<-
<-
2
name
name
1
D1
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
Y1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
I
name
name
1
M
name
name
1
Y1
<-
<-
2
name
name
1
T1
call
call
3
name
name
1
*
name
name
1
c
(
(
2
name
name
1
(
call
call
3
name
name
1
/
name
name
1
N1
name
name
1
D1
<-
<-
2
name
name
1
N2
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
Y2
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
I
name
name
1
M
name
name
1
N
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
I
name
name
1
M
name
name
1
Y2
<-
<-
2
name
name
1
D2
call
call
3
name
name
1
%*%
call
call
3
name
name
1
%*%
call
call
2
name
name
1
t
name
name
1
Y2
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
I
name
name
1
M
name
name
1
Y2
<-
<-
2
name
name
1
T2
call
call
3
name
name
1
*
name
name
1
c
(
(
2
name
name
1
(
call
call
3
name
name
1
/
name
name
1
N2
name
name
1
D2
<-
<-
2
name
name
1
p1
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
pchisq
name
name
1
T1
numeric
numeric
1
2
<-
<-
2
name
name
1
p2
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
pchisq
name
name
1
T2
numeric
numeric
1
2
<-
<-
2
name
name
1
pvalues
call
call
3
name
name
1
cbind
name
name
1
p1
name
name
1
p2
<-
<-
2
name
name
1
circ.lm
call
call
1
name
name
1
list
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
rho
name
name
1
rho
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
fitted
call
call
3
name
name
1
%%
name
name
1
theta.fitted
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
data
call
call
3
name
name
1
cbind
name
name
1
alpha
name
name
1
theta
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
residuals
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
theta.fitted
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
coef
call
call
3
name
name
1
cbind
call
call
3
name
name
1
$
name
name
1
cos.lm
character
character
1
coefficients
call
call
3
name
name
1
$
name
name
1
sin.lm
character
character
1
coefficients
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
pvalues
name
name
1
pvalues
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
A.k
call
call
2
name
name
1
mean
call
call
2
name
name
1
cos
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
residuals
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
kappa
call
call
2
name
name
1
A1inv
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
A.k
if
if
3
call
call
3
name
name
1
&
call
call
3
name
name
1
>
call
call
3
name
name
1
[
name
name
1
pvalues
numeric
numeric
1
1
name
name
1
level
call
call
3
name
name
1
>
call
call
3
name
name
1
[
name
name
1
pvalues
numeric
numeric
1
2
name
name
1
level
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
message
.named_I
structure
2
call
call
5
name
name
1
paste
character
character
1
Higher order terms are not significant at the 
name
name
1
level
character
character
1
 level
character
character
1

character
character
5




sep
<-
<-
2
call
call
3
name
name
1
$
name
name
1
circ.lm
character
character
1
message
.named_I
structure
2
call
call
5
name
name
1
paste
character
character
1
Higher order terms are significant at the 
name
name
1
level
character
character
1
 level
character
character
1

character
character
5




sep
name
name
1
circ.lm
character
character
5
alpha
theta
order
level

circ.summary
.named_I
structure
2
function
function
2
missing
missing
0
{
{
7
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
sinr
call
call
2
name
name
1
sum
call
call
2
name
name
1
sin
name
name
1
data
<-
<-
2
name
name
1
cosr
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
name
name
1
data
<-
<-
2
name
name
1
rho
call
call
3
name
name
1
/
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
sinr
numeric
numeric
1
2
call
call
3
name
name
1
^
name
name
1
cosr
numeric
numeric
1
2
name
name
1
n
<-
<-
2
name
name
1
mean.dir
call
call
3
name
name
1
atan
name
name
1
sinr
name
name
1
cosr
<-
<-
2
name
name
1
result
call
call
4
name
name
1
data.frame
name
name
1
n
name
name
1
mean.dir
name
name
1
rho
name
name
1
result
character
character
2
data

dcard
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
1
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
r
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
character
character
4
theta
mu
r

deg
.named_I
structure
2
function
function
2
missing
missing
0
{
{
1
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
name
name
1
radian
numeric
numeric
1
180
name
name
1
pi
character
character
2
radian

dmixedvm
.named_I
structure
2
function
function
7
missing
missing
0
missing
missing
0
missing
missing
0
missing
missing
0
missing
missing
0
missing
missing
0
{
{
1
call
call
3
name
name
1
+
call
call
3
name
name
1
*
call
call
3
name
name
1
/
name
name
1
p
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
2
name
name
1
I.0
name
name
1
kappa1
call
call
2
name
name
1
exp
call
call
3
name
name
1
*
name
name
1
kappa1
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu1
call
call
3
name
name
1
*
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
1
name
name
1
p
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
2
name
name
1
I.0
name
name
1
kappa2
call
call
2
name
name
1
exp
call
call
3
name
name
1
*
name
name
1
kappa2
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu2
character
character
7
theta
mu1
mu2
kappa1
kappa2
p

dtri
.named_I
structure
2
function
function
3
missing
missing
0
missing
missing
0
{
{
1
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
numeric
numeric
1
4
call
call
3
name
name
1
*
call
call
3
name
name
1
^
name
name
1
pi
numeric
numeric
1
2
name
name
1
r
call
call
3
name
name
1
*
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
name
name
1
r
call
call
2
name
name
1
abs
call
call
3
name
name
1
-
name
name
1
pi
name
name
1
theta
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
8
name
name
1
pi
character
character
3
theta
r

dvm
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
1
call
call
3
name
name
1
*
call
call
3
name
name
1
/
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
2
name
name
1
I.0
name
name
1
kappa
call
call
2
name
name
1
exp
call
call
3
name
name
1
*
name
name
1
kappa
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu
character
character
4
theta
mu
kappa

dwrpcauchy
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
1
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
^
name
name
1
rho
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
^
name
name
1
rho
numeric
numeric
1
2
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
rho
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu
character
character
4
theta
mu
rho

dwrpnorm
.named_I
structure
2
function
function
5
missing
missing
0
missing
missing
0
missing
missing
0
numeric
numeric
1
1.0000000000000001e-005
{
{
7
<-
<-
2
name
name
1
var
call
call
3
name
name
1
*
numeric
numeric
1
-2
call
call
2
name
name
1
log
name
name
1
rho
<-
<-
2
name
name
1
term
.named_I
structure
2
function
function
5
missing
missing
0
missing
missing
0
missing
missing
0
missing
missing
0
{
{
1
call
call
3
name
name
1
*
call
call
3
name
name
1
/
numeric
numeric
1
1
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
*
call
call
3
name
name
1
*
name
name
1
var
numeric
numeric
1
2
name
name
1
pi
call
call
2
name
name
1
exp
call
call
3
name
name
1
/
call
call
2
name
name
1
-
(
(
2
name
name
1
(
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
name
name
1
k
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
var
character
character
5
theta
mu
var
k

<-
<-
2
name
name
1
k
numeric
numeric
1
0
<-
<-
2
name
name
1
Next
call
call
5
name
name
1
term
name
name
1
theta
name
name
1
mu
name
name
1
var
name
name
1
k
<-
<-
2
name
name
1
delta
numeric
numeric
1
1
while
while
2
call
call
3
name
name
1
>
name
name
1
delta
name
name
1
acc
{
{
4
<-
<-
2
name
name
1
k
call
call
3
name
name
1
+
name
name
1
k
numeric
numeric
1
1
<-
<-
2
name
name
1
Last
name
name
1
Next
<-
<-
2
name
name
1
Next
call
call
3
name
name
1
+
call
call
3
name
name
1
+
name
name
1
Last
call
call
5
name
name
1
term
name
name
1
theta
name
name
1
mu
name
name
1
var
name
name
1
k
call
call
5
name
name
1
term
name
name
1
theta
name
name
1
mu
name
name
1
var
call
call
2
name
name
1
-
name
name
1
k
<-
<-
2
name
name
1
delta
call
call
2
name
name
1
abs
call
call
3
name
name
1
-
name
name
1
Next
name
name
1
Last
name
name
1
Next
character
character
5
theta
mu
rho
acc

eqscplot
.named_I
structure
2
function
function
7
missing
missing
0
missing
missing
0
numeric
numeric
1
0.040000000000000001
call
call
2
name
name
1
range
name
name
1
x
call
call
2
name
name
1
range
name
name
1
y
missing
missing
0
{
{
11
if
if
2
call
call
2
name
name
1
is.matrix
name
name
1
x
{
{
2
<-
<-
2
name
name
1
y
call
call
4
name
name
1
[
name
name
1
x
missing
missing
0
numeric
numeric
1
2
<-
<-
2
name
name
1
x
call
call
4
name
name
1
[
name
name
1
x
missing
missing
0
numeric
numeric
1
1
if
if
2
call
call
2
name
name
1
is.list
name
name
1
x
{
{
2
<-
<-
2
name
name
1
y
call
call
3
name
name
1
$
name
name
1
x
character
character
1
y
<-
<-
2
name
name
1
x
call
call
3
name
name
1
$
name
name
1
x
character
character
1
x
<-
<-
2
name
name
1
oldpin
call
call
2
name
name
1
par
character
character
1
pin
<-
<-
2
name
name
1
midx
call
call
3
name
name
1
*
numeric
numeric
1
0.5
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
xlim
numeric
numeric
1
2
call
call
3
name
name
1
[
name
name
1
xlim
numeric
numeric
1
1
<-
<-
2
name
name
1
xlim
call
call
3
name
name
1
+
name
name
1
midx
call
call
3
name
name
1
*
call
call
3
name
name
1
*
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
name
name
1
tol
numeric
numeric
1
0.5
call
call
3
name
name
1
c
numeric
numeric
1
-1
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
xlim
numeric
numeric
1
2
call
call
3
name
name
1
[
name
name
1
xlim
numeric
numeric
1
1
<-
<-
2
name
name
1
midy
call
call
3
name
name
1
*
numeric
numeric
1
0.5
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
ylim
numeric
numeric
1
2
call
call
3
name
name
1
[
name
name
1
ylim
numeric
numeric
1
1
<-
<-
2
name
name
1
ylim
call
call
3
name
name
1
+
name
name
1
midy
call
call
3
name
name
1
*
call
call
3
name
name
1
*
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
name
name
1
tol
numeric
numeric
1
0.5
call
call
3
name
name
1
c
numeric
numeric
1
-1
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
ylim
numeric
numeric
1
2
call
call
3
name
name
1
[
name
name
1
ylim
numeric
numeric
1
1
<-
<-
2
name
name
1
xr
call
call
3
name
name
1
/
call
call
3
name
name
1
[
name
name
1
oldpin
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
xlim
numeric
numeric
1
2
call
call
3
name
name
1
[
name
name
1
xlim
numeric
numeric
1
1
<-
<-
2
name
name
1
yr
call
call
3
name
name
1
/
call
call
3
name
name
1
[
name
name
1
oldpin
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
ylim
numeric
numeric
1
2
call
call
3
name
name
1
[
name
name
1
ylim
numeric
numeric
1
1
if
if
3
call
call
3
name
name
1
>
name
name
1
yr
name
name
1
xr
{
{
1
<-
<-
2
name
name
1
ylim
call
call
3
name
name
1
+
name
name
1
midy
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
name
name
1
yr
call
call
3
name
name
1
c
numeric
numeric
1
-1
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
ylim
numeric
numeric
1
2
call
call
3
name
name
1
[
name
name
1
ylim
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
xr
{
{
1
<-
<-
2
name
name
1
xlim
call
call
3
name
name
1
+
name
name
1
midx
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
name
name
1
xr
call
call
3
name
name
1
c
numeric
numeric
1
-1
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
xlim
numeric
numeric
1
2
call
call
3
name
name
1
[
name
name
1
xlim
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
yr
call with ...
call with ...
1
.named_I
structure
2
call
call
8
name
name
1
plot
name
name
1
x
name
name
1
y
name
name
1
xlim
name
name
1
ylim
character
character
1
i
character
character
1
i
name
name
1
...
character
character
8



xlim
ylim
xaxs
yaxs

character
character
7
x
y
tol
xlim
ylim
...

equal.kappa.test
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
name
name
1
T
{
{
17
<-
<-
2
name
name
1
ns
.named_I
structure
2
call
call
4
name
name
1
tapply
name
name
1
theta
name
name
1
group
name
name
1
length
character
character
4



FUN
<-
<-
2
name
name
1
r.bars
.named_I
structure
2
call
call
4
name
name
1
tapply
name
name
1
theta
name
name
1
group
name
name
1
est.rho
character
character
4



FUN
<-
<-
2
name
name
1
rs
call
call
3
name
name
1
*
name
name
1
r.bars
name
name
1
ns
<-
<-
2
name
name
1
kappas
.named_I
structure
2
call
call
4
name
name
1
tapply
name
name
1
theta
name
name
1
group
name
name
1
est.kappa
character
character
4



FUN
<-
<-
2
name
name
1
grps
call
call
2
name
name
1
length
name
name
1
r.bars
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
group
<-
<-
2
name
name
1
r.bar.all
call
call
2
name
name
1
est.rho
name
name
1
theta
<-
<-
2
name
name
1
kappa.all
call
call
2
name
name
1
est.kappa
name
name
1
theta
<-
<-
2
name
name
1
warn1
integer
integer
1
0
if
if
2
call
call
3
name
name
1
<
name
name
1
r.bar.all
numeric
numeric
1
0.45000000000000001
{
{
5
<-
<-
2
name
name
1
g1
.named_I
structure
2
function
function
2
missing
missing
0
{
{
1
call
call
2
name
name
1
asin
call
call
3
name
name
1
*
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
/
integer
integer
1
3
integer
integer
1
8
name
name
1
x
character
character
2
x

<-
<-
2
name
name
1
ws
call
call
3
name
name
1
/
call
call
3
name
name
1
*
integer
integer
1
4
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
ns
integer
integer
1
4
integer
integer
1
3
<-
<-
2
name
name
1
g1s
call
call
2
name
name
1
g1
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
r.bars
<-
<-
2
name
name
1
U
call
call
3
name
name
1
-
call
call
2
name
name
1
sum
call
call
3
name
name
1
*
name
name
1
ws
call
call
3
name
name
1
^
name
name
1
g1s
integer
integer
1
2
call
call
3
name
name
1
/
call
call
3
name
name
1
^
call
call
2
name
name
1
sum
call
call
3
name
name
1
*
name
name
1
ws
name
name
1
g1s
integer
integer
1
2
call
call
2
name
name
1
sum
name
name
1
ws
if
if
2
call
call
2
name
name
1
sum
call
call
3
name
name
1
>
call
call
2
name
name
1
is.na
name
name
1
g1s
integer
integer
1
0
{
{
2
<-
<-
2
name
name
1
warn1
integer
integer
1
1
call
call
2
name
name
1
warning
character
character
1
An argument outide of [-1,1] was passed to asin function in calculation of approximate chi-squared test statistic.  \n            Bartlett's test of homogeneity was used instead of the approximation using asin.
if
if
2
call
call
3
name
name
1
&
call
call
3
name
name
1
>=
name
name
1
r.bar.all
numeric
numeric
1
0.45000000000000001
call
call
3
name
name
1
<=
name
name
1
r.bar.all
numeric
numeric
1
0.69999999999999996
{
{
4
<-
<-
2
name
name
1
g2
.named_I
structure
2
function
function
2
missing
missing
0
{
{
3
<-
<-
2
name
name
1
c1
numeric
numeric
1
1.089
<-
<-
2
name
name
1
c2
numeric
numeric
1
0.25800000000000001
call
call
2
name
name
1
asinh
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
x
name
name
1
c1
name
name
1
c2
character
character
2
x

<-
<-
2
name
name
1
ws
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
ns
integer
integer
1
3
numeric
numeric
1
0.79800000000000004
<-
<-
2
name
name
1
g2s
call
call
2
name
name
1
g2
name
name
1
r.bars
<-
<-
2
name
name
1
U
call
call
3
name
name
1
-
call
call
2
name
name
1
sum
call
call
3
name
name
1
*
name
name
1
ws
call
call
3
name
name
1
^
name
name
1
g2s
integer
integer
1
2
call
call
3
name
name
1
/
call
call
3
name
name
1
^
call
call
2
name
name
1
sum
call
call
3
name
name
1
*
name
name
1
ws
name
name
1
g2s
integer
integer
1
2
call
call
2
name
name
1
sum
name
name
1
ws
if
if
2
call
call
3
name
name
1
|
call
call
3
name
name
1
>
name
name
1
r.bar.all
numeric
numeric
1
0.69999999999999996
call
call
3
name
name
1
==
name
name
1
warn1
integer
integer
1
1
{
{
4
<-
<-
2
name
name
1
vs
call
call
3
name
name
1
-
name
name
1
ns
integer
integer
1
1
<-
<-
2
name
name
1
v
call
call
3
name
name
1
-
name
name
1
n
name
name
1
grps
<-
<-
2
name
name
1
d
call
call
3
name
name
1
*
call
call
3
name
name
1
/
integer
integer
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
3
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
grps
integer
integer
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
integer
integer
1
1
name
name
1
vs
call
call
3
name
name
1
/
integer
integer
1
1
name
name
1
v
<-
<-
2
name
name
1
U
call
call
3
name
name
1
*
call
call
3
name
name
1
/
integer
integer
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
+
integer
integer
1
1
name
name
1
d
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
*
name
name
1
v
call
call
2
name
name
1
log
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
call
call
2
name
name
1
sum
name
name
1
rs
name
name
1
v
call
call
2
name
name
1
sum
call
call
3
name
name
1
*
name
name
1
vs
call
call
2
name
name
1
log
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
ns
name
name
1
rs
name
name
1
vs
<-
<-
2
name
name
1
p.value
call
call
3
name
name
1
-
integer
integer
1
1
call
call
3
name
name
1
pchisq
name
name
1
U
call
call
3
name
name
1
-
name
name
1
grps
integer
integer
1
1
<-
<-
2
name
name
1
result.vector
.named_I
structure
2
call
call
4
name
name
1
c
call
call
3
name
name
1
-
name
name
1
grps
integer
integer
1
1
name
name
1
U
name
name
1
p.value
character
character
4

df
ChiSq
p
<-
<-
2
name
name
1
result.list
.named_I
structure
2
call
call
8
name
name
1
list
name
name
1
kappas
name
name
1
kappa.all
name
name
1
r.bars
name
name
1
r.bar.all
call
call
3
name
name
1
-
name
name
1
grps
integer
integer
1
1
name
name
1
U
name
name
1
p.value
character
character
8

kappa
kappa.all
r.mean
r.mean.all
df
chisq.stat
p.value
if
if
2
call
call
3
name
name
1
==
name
name
1
print
name
name
1
T
{
{
3
call
call
5
name
name
1
cat
character
character
1
\n
character
character
1
Test for Homogeneity of Concentration Parameters
character
character
1
\n
character
character
1
\n
call
call
2
name
name
1
print
name
name
1
result.vector
call
call
3
name
name
1
cat
character
character
1
\n
character
character
1
\n
name
name
1
result.list
character
character
4
group
theta
print

est.kappa
.named_I
structure
2
function
function
3
missing
missing
0
name
name
1
F
{
{
4
<-
<-
2
name
name
1
mean.dir
call
call
2
name
name
1
circ.mean
name
name
1
data
<-
<-
2
name
name
1
kappa
call
call
2
name
name
1
A1inv
call
call
2
name
name
1
mean
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
data
name
name
1
mean.dir
if
if
2
call
call
3
name
name
1
==
name
name
1
bias
name
name
1
T
{
{
4
<-
<-
2
name
name
1
kappa.ml
name
name
1
kappa
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
if
if
2
call
call
3
name
name
1
<
name
name
1
kappa.ml
numeric
numeric
1
2
<-
<-
2
name
name
1
kappa
call
call
3
name
name
1
max
call
call
3
name
name
1
-
name
name
1
kappa.ml
call
call
3
name
name
1
*
numeric
numeric
1
2
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
*
name
name
1
n
name
name
1
kappa.ml
numeric
numeric
1
-1
numeric
numeric
1
0
if
if
2
call
call
3
name
name
1
>=
name
name
1
kappa.ml
numeric
numeric
1
2
<-
<-
2
name
name
1
kappa
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
numeric
numeric
1
1
numeric
numeric
1
3
name
name
1
kappa.ml
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
n
numeric
numeric
1
3
name
name
1
n
name
name
1
kappa
character
character
3
data
bias

est.rho
.named_I
structure
2
function
function
2
missing
missing
0
{
{
4
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
sinr
call
call
2
name
name
1
sum
call
call
2
name
name
1
sin
name
name
1
data
<-
<-
2
name
name
1
cosr
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
name
name
1
data
call
call
3
name
name
1
/
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
sinr
numeric
numeric
1
2
call
call
3
name
name
1
^
name
name
1
cosr
numeric
numeric
1
2
name
name
1
n
character
character
2
data

kuiper
.named_I
structure
2
function
function
3
missing
missing
0
numeric
numeric
1
0
{
{
10
call
call
5
name
name
1
cat
character
character
1
\n
character
character
1
      Kuiper's Test of Uniformity
character
character
1
\n
character
character
1
\n
<-
<-
2
name
name
1
kuiper.crits
call
call
3
name
name
1
cbind
call
call
6
name
name
1
c
numeric
numeric
1
0.14999999999999999
numeric
numeric
1
0.10000000000000001
numeric
numeric
1
0.050000000000000003
numeric
numeric
1
0.025000000000000001
numeric
numeric
1
0.01
call
call
6
name
name
1
c
numeric
numeric
1
1.5369999999999999
numeric
numeric
1
1.6200000000000001
numeric
numeric
1
1.7470000000000001
numeric
numeric
1
1.8620000000000001
numeric
numeric
1
2.0009999999999999
<-
<-
2
name
name
1
x
call
call
3
name
name
1
/
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
data
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
i
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
<-
<-
2
name
name
1
D.P
call
call
2
name
name
1
max
call
call
3
name
name
1
-
call
call
3
name
name
1
/
name
name
1
i
name
name
1
n
name
name
1
x
<-
<-
2
name
name
1
D.M
call
call
2
name
name
1
max
call
call
3
name
name
1
-
name
name
1
x
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
i
numeric
numeric
1
1
name
name
1
n
<-
<-
2
name
name
1
V
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
D.P
name
name
1
D.M
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
+
call
call
2
name
name
1
sqrt
name
name
1
n
numeric
numeric
1
0.155
call
call
3
name
name
1
/
numeric
numeric
1
0.23999999999999999
call
call
2
name
name
1
sqrt
name
name
1
n
call
call
4
name
name
1
cat
character
character
1
Test Statistic:
call
call
3
name
name
1
round
name
name
1
V
numeric
numeric
1
4
character
character
1
\n
if
if
3
call
call
3
name
name
1
==
name
name
1
alpha
numeric
numeric
1
0
{
{
1
if
if
3
call
call
3
name
name
1
<
name
name
1
V
numeric
numeric
1
1.5369999999999999
call
call
4
name
name
1
cat
character
character
1
P-value > 0.15
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
<
name
name
1
V
numeric
numeric
1
1.6200000000000001
call
call
4
name
name
1
cat
character
character
1
0.10 < P-value < 0.15
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
<
name
name
1
V
numeric
numeric
1
1.7470000000000001
call
call
4
name
name
1
cat
character
character
1
0.05 < P-value < 0.10
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
<
name
name
1
V
numeric
numeric
1
1.8620000000000001
call
call
4
name
name
1
cat
character
character
1
0.025 < P-value < 0.05
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
<
name
name
1
V
numeric
numeric
1
2.0009999999999999
call
call
4
name
name
1
cat
character
character
1
0.01 < P-value < 0.025
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
P-value < 0.01
character
character
1
\n
character
character
1
\n
{
{
3
<-
<-
2
name
name
1
Critical
call
call
4
name
name
1
[
name
name
1
kuiper.crits
call
call
3
name
name
1
[
(
(
2
name
name
1
(
call
call
3
name
name
1
:
numeric
numeric
1
1
numeric
numeric
1
5
call
call
3
name
name
1
==
name
name
1
alpha
call
call
2
name
name
1
c
call
call
4
name
name
1
[
name
name
1
kuiper.crits
missing
missing
0
numeric
numeric
1
1
numeric
numeric
1
2
call
call
6
name
name
1
cat
character
character
1
Level
name
name
1
alpha
character
character
1
Critical Value:
call
call
3
name
name
1
round
name
name
1
Critical
numeric
numeric
1
4
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
V
name
name
1
Critical
call
call
4
name
name
1
cat
character
character
1
Reject Null Hypothesis
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Do Not Reject Null Hypothesis
character
character
1
\n
character
character
1
\n
character
character
3
data
alpha

last.dump
.named_I
structure
2
structure
structure
3
list
list
4
character
character
1
No Frame Available
character
character
1
No Frame Available
character
character
1
No Frame Available
character
character
1
No Frame Available
character
character
4

eval(expression(dir))
doErrorAction("Problem: Object \\"dir\\" not found", 1000)
eval(action, sys.parent())
character
character
1
Problem: Object "dir" not found
character
character
3
.Data
.Names
message
ls.circle
.named_I
structure
2
function
function
5
missing
missing
0
missing
missing
0
name
name
1
NULL
missing
missing
0
{
{
4
if
if
2
call
call
2
name
name
1
is.null
name
name
1
p.0
{
{
4
<-
<-
2
name
name
1
x.mean
call
call
2
name
name
1
mean
name
name
1
x
<-
<-
2
name
name
1
y.mean
call
call
2
name
name
1
mean
name
name
1
y
<-
<-
2
name
name
1
est.r
call
call
2
name
name
1
max
call
call
3
name
name
1
c
call
call
2
name
name
1
abs
call
call
3
name
name
1
-
name
name
1
x
name
name
1
x.mean
call
call
2
name
name
1
abs
call
call
3
name
name
1
-
name
name
1
y
name
name
1
y.mean
<-
<-
2
name
name
1
p.0
call
call
4
name
name
1
c
name
name
1
est.r
name
name
1
x.mean
name
name
1
y.mean
<-
<-
2
name
name
1
obj.fun
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
1
call
call
2
name
name
1
sum
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
p
integer
integer
1
1
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
x
call
call
3
name
name
1
[
name
name
1
p
integer
integer
1
2
integer
integer
1
2
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
y
call
call
3
name
name
1
[
name
name
1
p
integer
integer
1
3
integer
integer
1
2
integer
integer
1
2
character
character
4
p
x
y

<-
<-
2
name
name
1
grad.fun
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
15
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
x
<-
<-
2
name
name
1
r
call
call
3
name
name
1
[
name
name
1
p
integer
integer
1
1
<-
<-
2
name
name
1
a
call
call
3
name
name
1
[
name
name
1
p
integer
integer
1
2
<-
<-
2
name
name
1
b
call
call
3
name
name
1
[
name
name
1
p
integer
integer
1
3
<-
<-
2
name
name
1
common
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
x
name
name
1
a
integer
integer
1
2
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
y
name
name
1
b
integer
integer
1
2
<-
<-
2
name
name
1
g.e1
call
call
3
name
name
1
-
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
n
name
name
1
r
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
sum
name
name
1
common
<-
<-
2
name
name
1
g.e2
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
n
name
name
1
a
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
n
call
call
2
name
name
1
mean
name
name
1
x
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
r
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
x
name
name
1
a
name
name
1
common
<-
<-
2
name
name
1
g.e3
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
n
name
name
1
b
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
n
call
call
2
name
name
1
mean
name
name
1
y
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
r
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
y
name
name
1
b
name
name
1
common
<-
<-
2
name
name
1
h.e11
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
n
<-
<-
2
name
name
1
h.e21
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
x
name
name
1
a
name
name
1
common
<-
<-
2
name
name
1
h.e22
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
n
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
r
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
y
name
name
1
b
integer
integer
1
2
call
call
3
name
name
1
^
name
name
1
common
integer
integer
1
3
<-
<-
2
name
name
1
h.e31
call
call
3
name
name
1
*
integer
integer
1
2
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
y
name
name
1
b
name
name
1
common
<-
<-
2
name
name
1
h.e32
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
r
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
x
name
name
1
a
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
y
name
name
1
b
call
call
3
name
name
1
^
name
name
1
common
integer
integer
1
3
<-
<-
2
name
name
1
h.e33
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
n
call
call
3
name
name
1
*
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
r
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
x
name
name
1
a
integer
integer
1
2
call
call
3
name
name
1
^
name
name
1
common
integer
integer
1
3
.named_I
structure
2
call
call
3
name
name
1
list
call
call
4
name
name
1
c
name
name
1
g.e1
name
name
1
g.e2
name
name
1
g.e3
call
call
7
name
name
1
c
name
name
1
h.e11
name
name
1
h.e21
name
name
1
h.e22
name
name
1
h.e31
name
name
1
h.e32
name
name
1
h.e33
character
character
3

gradient
hessian
character
character
4
p
x
y

call with ...
call with ...
1
.named_I
structure
2
call
call
8
name
name
1
nlminb
name
name
1
p.0
name
name
1
obj.fun
name
name
1
grad.fun
name
name
1
T
name
name
1
x
name
name
1
y
name
name
1
...
character
character
8

start
obj
gradient
hessian
x
y

character
character
5
x
y
p.0
...

plot.edf
.named_I
structure
2
function
function
3
missing
missing
0
missing
missing
0
{
{
4
<-
<-
2
name
name
1
data
call
call
3
name
name
1
%%
name
name
1
data
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
data
call
call
2
name
name
1
sort
name
name
1
data
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
call with ...
call with ...
1
.named_I
structure
2
call
call
7
name
name
1
plot
call
call
4
name
name
1
c
numeric
numeric
1
0
name
name
1
data
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
4
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
/
call
call
2
name
name
1
seq
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
name
name
1
n
numeric
numeric
1
1
character
character
1
s
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
3
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
1
name
name
1
...
character
character
7



type
xlim
ylim

character
character
3
data
...

plot.ls.circle
.named_I
structure
2
function
function
10
missing
missing
0
missing
missing
0
missing
missing
0
missing
missing
0
missing
missing
0
integer
integer
1
16
integer
integer
1
1
integer
integer
1
1
missing
missing
0
{
{
6
<-
<-
2
name
name
1
xvals
call
call
3
name
name
1
c
call
call
3
name
name
1
+
name
name
1
a
call
call
3
name
name
1
*
name
name
1
r
call
call
2
name
name
1
cos
.named_I
structure
2
call
call
4
name
name
1
seq
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
numeric
numeric
1
1000
character
character
4



length
name
name
1
x
<-
<-
2
name
name
1
yvals
call
call
3
name
name
1
c
call
call
3
name
name
1
+
name
name
1
b
call
call
3
name
name
1
*
name
name
1
r
call
call
2
name
name
1
sin
.named_I
structure
2
call
call
4
name
name
1
seq
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
numeric
numeric
1
1000
character
character
4



length
name
name
1
y
call with ...
call with ...
1
.named_I
structure
2
call
call
5
name
name
1
eqscplot
call
call
3
name
name
1
c
call
call
2
name
name
1
min
call
call
3
name
name
1
c
name
name
1
xvals
name
name
1
yvals
call
call
2
name
name
1
max
call
call
3
name
name
1
c
name
name
1
xvals
name
name
1
yvals
call
call
3
name
name
1
c
call
call
2
name
name
1
min
call
call
3
name
name
1
c
name
name
1
xvals
name
name
1
yvals
call
call
2
name
name
1
max
call
call
3
name
name
1
c
name
name
1
xvals
name
name
1
yvals
character
character
1
n
name
name
1
...
character
character
5



type

.named_I
structure
2
call
call
4
name
name
1
points
call
call
3
name
name
1
+
name
name
1
a
call
call
3
name
name
1
*
name
name
1
r
call
call
2
name
name
1
cos
.named_I
structure
2
call
call
4
name
name
1
seq
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
numeric
numeric
1
1000
character
character
4



length
call
call
3
name
name
1
+
name
name
1
b
call
call
3
name
name
1
*
name
name
1
r
call
call
2
name
name
1
sin
.named_I
structure
2
call
call
4
name
name
1
seq
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
numeric
numeric
1
1000
character
character
4



length
character
character
1
l
character
character
4



type
.named_I
structure
2
call
call
5
name
name
1
points
name
name
1
x
name
name
1
y
name
name
1
points.pch
name
name
1
points.cex
character
character
5



pch
cex
.named_I
structure
2
call
call
5
name
name
1
text
name
name
1
a
name
name
1
b
character
character
1
+
name
name
1
plus.cex
character
character
5




cex
character
character
10
x
y
r
a
b
points.pch
points.cex
plus.cex
...

pp.plot
.named_I
structure
2
function
function
3
missing
missing
0
name
name
1
T
{
{
10
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
data
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
data
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
x
call
call
3
name
name
1
/
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
n
numeric
numeric
1
1
<-
<-
2
name
name
1
mu
call
call
3
name
name
1
%%
call
call
2
name
name
1
circ.mean
name
name
1
data
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
kappa
call
call
2
name
name
1
est.kappa
name
name
1
data
<-
<-
2
name
name
1
y
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
{
{
1
<-
<-
2
call
call
3
name
name
1
[
name
name
1
y
name
name
1
i
call
call
4
name
name
1
pvm
call
call
3
name
name
1
[
name
name
1
data
name
name
1
i
name
name
1
mu
name
name
1
kappa
character
character
3
i


.named_I
structure
2
call
call
5
name
name
1
plot
name
name
1
x
name
name
1
y
character
character
1
von Mises Distribution
character
character
1
Empirical Distribution
character
character
5



xlab
ylab
if
if
2
name
name
1
ref.line
call
call
3
name
name
1
abline
numeric
numeric
1
0
numeric
numeric
1
1
call
call
3
name
name
1
data.frame
name
name
1
mu
name
name
1
kappa
character
character
3
data
ref.line

pvm
.named_I
structure
2
function
function
5
missing
missing
0
missing
missing
0
missing
missing
0
numeric
numeric
1
9.9999999999999995e-021
{
{
5
<-
<-
2
name
name
1
theta
call
call
3
name
name
1
%%
name
name
1
theta
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
mu
call
call
3
name
name
1
%%
name
name
1
mu
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
pvm.mu0
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
5
<-
<-
2
name
name
1
flag
character
character
1
true
<-
<-
2
name
name
1
p
numeric
numeric
1
1
<-
<-
2
name
name
1
sum
numeric
numeric
1
0
while
while
2
call
call
3
name
name
1
==
name
name
1
flag
character
character
1
true
{
{
4
<-
<-
2
name
name
1
term
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
I.p
name
name
1
p
name
name
1
kappa
call
call
2
name
name
1
sin
call
call
3
name
name
1
*
name
name
1
p
name
name
1
theta
name
name
1
p
<-
<-
2
name
name
1
sum
call
call
3
name
name
1
+
name
name
1
sum
name
name
1
term
<-
<-
2
name
name
1
p
call
call
3
name
name
1
+
name
name
1
p
numeric
numeric
1
1
if
if
2
call
call
3
name
name
1
<
call
call
2
name
name
1
abs
name
name
1
term
name
name
1
acc
<-
<-
2
name
name
1
flag
character
character
1
false
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
theta
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
3
name
name
1
/
name
name
1
sum
(
(
2
name
name
1
(
call
call
3
name
name
1
*
name
name
1
pi
call
call
2
name
name
1
I.0
name
name
1
kappa
character
character
4
theta
kappa
acc

if
if
3
call
call
3
name
name
1
==
name
name
1
mu
numeric
numeric
1
0
{
{
1
<-
<-
2
name
name
1
result
call
call
4
name
name
1
pvm.mu0
name
name
1
theta
name
name
1
kappa
name
name
1
acc
{
{
1
if
if
3
call
call
3
name
name
1
<=
name
name
1
theta
name
name
1
mu
{
{
4
<-
<-
2
name
name
1
upper
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
if
if
2
call
call
3
name
name
1
==
name
name
1
upper
numeric
numeric
1
0
<-
<-
2
name
name
1
upper
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
lower
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
2
name
name
1
-
name
name
1
mu
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
result
call
call
3
name
name
1
-
call
call
4
name
name
1
pvm.mu0
name
name
1
upper
name
name
1
kappa
name
name
1
acc
call
call
4
name
name
1
pvm.mu0
name
name
1
lower
name
name
1
kappa
name
name
1
acc
{
{
3
<-
<-
2
name
name
1
upper
call
call
3
name
name
1
-
name
name
1
theta
name
name
1
mu
<-
<-
2
name
name
1
lower
call
call
3
name
name
1
%%
name
name
1
mu
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
result
call
call
3
name
name
1
+
call
call
4
name
name
1
pvm.mu0
name
name
1
upper
name
name
1
kappa
name
name
1
acc
call
call
4
name
name
1
pvm.mu0
name
name
1
lower
name
name
1
kappa
name
name
1
acc
name
name
1
result
character
character
5
theta
mu
kappa
acc

r.test
.named_I
structure
2
function
function
2
missing
missing
0
{
{
9
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
ss
call
call
2
name
name
1
sum
call
call
2
name
name
1
sin
name
name
1
data
<-
<-
2
name
name
1
cc
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
name
name
1
data
<-
<-
2
name
name
1
rbar
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
ss
numeric
numeric
1
2
call
call
3
name
name
1
^
name
name
1
cc
numeric
numeric
1
2
name
name
1
n
<-
<-
2
name
name
1
z
(
(
2
name
name
1
(
call
call
3
name
name
1
*
name
name
1
n
call
call
3
name
name
1
^
name
name
1
rbar
numeric
numeric
1
2
<-
<-
2
name
name
1
p.value
call
call
2
name
name
1
exp
call
call
2
name
name
1
-
name
name
1
z
if
if
3
call
call
3
name
name
1
<
name
name
1
n
numeric
numeric
1
50
<-
<-
2
name
name
1
temp
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
z
call
call
3
name
name
1
^
name
name
1
z
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
4
name
name
1
n
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
numeric
numeric
1
24
name
name
1
z
call
call
3
name
name
1
*
numeric
numeric
1
132
call
call
3
name
name
1
^
name
name
1
z
numeric
numeric
1
2
call
call
3
name
name
1
*
numeric
numeric
1
76
call
call
3
name
name
1
^
name
name
1
z
numeric
numeric
1
3
call
call
3
name
name
1
*
numeric
numeric
1
9
call
call
3
name
name
1
^
name
name
1
z
numeric
numeric
1
4
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
288
call
call
3
name
name
1
^
name
name
1
n
numeric
numeric
1
2
<-
<-
2
name
name
1
temp
numeric
numeric
1
1
<-
<-
2
name
name
1
result
.named_I
structure
2
call
call
3
name
name
1
list
name
name
1
rbar
call
call
3
name
name
1
*
name
name
1
p.value
name
name
1
temp
character
character
3

r.bar
p.value
name
name
1
result
character
character
2
data

rad
.named_I
structure
2
function
function
2
missing
missing
0
{
{
1
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
name
name
1
degree
name
name
1
pi
numeric
numeric
1
180
character
character
2
degree

rao.homogeneity
.named_I
structure
2
function
function
3
missing
missing
0
numeric
numeric
1
0
{
{
20
if
if
2
call
call
2
name
name
1
!
call
call
2
name
name
1
is.list
name
name
1
data
call
call
2
name
name
1
stop
character
character
1
Data must be of mode list
<-
<-
2
name
name
1
n
call
call
2
name
name
1
unlist
call
call
3
name
name
1
lapply
name
name
1
data
name
name
1
length
<-
<-
2
name
name
1
k
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
c.data
call
call
3
name
name
1
lapply
name
name
1
data
name
name
1
cos
<-
<-
2
name
name
1
s.data
call
call
3
name
name
1
lapply
name
name
1
data
name
name
1
sin
<-
<-
2
name
name
1
x
call
call
2
name
name
1
unlist
call
call
3
name
name
1
lapply
name
name
1
c.data
name
name
1
mean
<-
<-
2
name
name
1
y
call
call
2
name
name
1
unlist
call
call
3
name
name
1
lapply
name
name
1
s.data
name
name
1
mean
<-
<-
2
name
name
1
s.co
call
call
2
name
name
1
unlist
call
call
3
name
name
1
lapply
name
name
1
c.data
name
name
1
var
<-
<-
2
name
name
1
s.ss
call
call
2
name
name
1
unlist
call
call
3
name
name
1
lapply
name
name
1
s.data
name
name
1
var
<-
<-
2
name
name
1
s.cs
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
k
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
k
{
{
1
<-
<-
2
call
call
3
name
name
1
[
name
name
1
s.cs
name
name
1
i
call
call
3
name
name
1
var
call
call
3
name
name
1
[[
name
name
1
c.data
name
name
1
i
call
call
3
name
name
1
[[
name
name
1
s.data
name
name
1
i
character
character
3
i


<-
<-
2
name
name
1
s.polar
call
call
3
name
name
1
*
call
call
3
name
name
1
/
numeric
numeric
1
1
name
name
1
n
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
s.ss
call
call
3
name
name
1
^
name
name
1
x
numeric
numeric
1
2
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
^
name
name
1
y
numeric
numeric
1
2
name
name
1
s.co
call
call
3
name
name
1
^
name
name
1
x
numeric
numeric
1
4
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
y
name
name
1
s.cs
call
call
3
name
name
1
^
name
name
1
x
numeric
numeric
1
3
<-
<-
2
name
name
1
tan
call
call
3
name
name
1
/
name
name
1
y
name
name
1
x
<-
<-
2
name
name
1
H.polar
call
call
3
name
name
1
-
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
call
call
3
name
name
1
^
name
name
1
tan
numeric
numeric
1
2
name
name
1
s.polar
call
call
3
name
name
1
/
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
name
name
1
tan
name
name
1
s.polar
numeric
numeric
1
2
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
numeric
numeric
1
1
name
name
1
s.polar
<-
<-
2
name
name
1
U
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
x
numeric
numeric
1
2
call
call
3
name
name
1
^
name
name
1
y
numeric
numeric
1
2
<-
<-
2
name
name
1
s.disp
call
call
3
name
name
1
*
call
call
3
name
name
1
/
numeric
numeric
1
4
name
name
1
n
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
+
call
call
3
name
name
1
*
call
call
3
name
name
1
^
name
name
1
x
numeric
numeric
1
2
name
name
1
s.co
call
call
3
name
name
1
*
call
call
3
name
name
1
^
name
name
1
y
numeric
numeric
1
2
name
name
1
s.ss
call
call
3
name
name
1
*
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
x
name
name
1
y
name
name
1
s.cs
<-
<-
2
name
name
1
H.disp
call
call
3
name
name
1
-
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
call
call
3
name
name
1
^
name
name
1
U
numeric
numeric
1
2
name
name
1
s.disp
call
call
3
name
name
1
/
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
name
name
1
U
name
name
1
s.disp
numeric
numeric
1
2
call
call
2
name
name
1
sum
call
call
3
name
name
1
/
numeric
numeric
1
1
name
name
1
s.disp
call
call
2
name
name
1
cat
character
character
1
\n
call
call
3
name
name
1
cat
character
character
1
Rao's Tests for Homogeneity
character
character
1
\n
if
if
3
call
call
3
name
name
1
==
name
name
1
alpha
numeric
numeric
1
0
{
{
9
call
call
2
name
name
1
cat
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
       Test for Equality of Polar Vectors:
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Test Statistic =
call
call
3
name
name
1
round
name
name
1
H.polar
numeric
numeric
1
5
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Degrees of Freedom =
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
character
character
1
\n
call
call
5
name
name
1
cat
character
character
1
P-value of test =
call
call
3
name
name
1
round
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
pchisq
name
name
1
H.polar
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
numeric
numeric
1
5
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
       Test for Equality of Dispersions:
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Test Statistic =
call
call
3
name
name
1
round
name
name
1
H.disp
numeric
numeric
1
5
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Degrees of Freedom =
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
character
character
1
\n
call
call
5
name
name
1
cat
character
character
1
P-value of test =
call
call
3
name
name
1
round
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
pchisq
name
name
1
H.disp
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
numeric
numeric
1
5
character
character
1
\n
character
character
1
\n
{
{
11
call
call
2
name
name
1
cat
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
       Test for Equality of Polar Vectors:
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Test Statistic =
call
call
3
name
name
1
round
name
name
1
H.polar
numeric
numeric
1
5
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Degrees of Freedom =
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
character
character
1
\n
call
call
6
name
name
1
cat
character
character
1
Level
name
name
1
alpha
character
character
1
critical value =
call
call
3
name
name
1
round
call
call
3
name
name
1
qchisq
call
call
3
name
name
1
-
numeric
numeric
1
1
name
name
1
alpha
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
numeric
numeric
1
5
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
H.polar
call
call
3
name
name
1
qchisq
call
call
3
name
name
1
-
numeric
numeric
1
1
name
name
1
alpha
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
call
call
4
name
name
1
cat
character
character
1
Reject null hypothesis of equal polar vectors
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Do not reject null hypothesis of equal polar vectors
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
       Test for Equality of Dispersions:
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Test Statistic =
call
call
3
name
name
1
round
name
name
1
H.disp
numeric
numeric
1
5
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Degrees of Freedom =
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
character
character
1
\n
call
call
6
name
name
1
cat
character
character
1
Level
name
name
1
alpha
character
character
1
critical value =
call
call
3
name
name
1
round
call
call
3
name
name
1
qchisq
call
call
3
name
name
1
-
numeric
numeric
1
1
name
name
1
alpha
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
numeric
numeric
1
5
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
H.disp
call
call
3
name
name
1
qchisq
call
call
3
name
name
1
-
numeric
numeric
1
1
name
name
1
alpha
call
call
3
name
name
1
-
name
name
1
k
numeric
numeric
1
1
call
call
4
name
name
1
cat
character
character
1
Reject null hypothesis of equal dispersions
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Do not reject null hypothesis of equal dispersions
character
character
1
\n
character
character
1
\n
character
character
3
data
alpha

rao.spacing
.named_I
structure
2
function
function
3
missing
missing
0
integer
integer
1
0
{
{
8
<-
<-
2
name
name
1
data
call
call
2
name
name
1
deg
name
name
1
data
<-
<-
2
name
name
1
data
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
data
numeric
numeric
1
360
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
spacings
call
call
3
name
name
1
c
call
call
2
name
name
1
diff
name
name
1
data
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
data
numeric
numeric
1
1
call
call
3
name
name
1
[
name
name
1
data
name
name
1
n
numeric
numeric
1
360
<-
<-
2
name
name
1
U
call
call
3
name
name
1
*
call
call
3
name
name
1
/
numeric
numeric
1
1
numeric
numeric
1
2
call
call
2
name
name
1
sum
call
call
2
name
name
1
abs
call
call
3
name
name
1
-
name
name
1
spacings
call
call
3
name
name
1
/
numeric
numeric
1
360
name
name
1
n
if
if
2
call
call
3
name
name
1
<
name
name
1
n
numeric
numeric
1
4
call
call
2
name
name
1
stop
character
character
1
Sample size too small
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
30
<-
<-
2
name
name
1
table.row
call
call
3
name
name
1
-
name
name
1
n
numeric
numeric
1
3
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
32
<-
<-
2
name
name
1
table.row
numeric
numeric
1
27
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
37
<-
<-
2
name
name
1
table.row
numeric
numeric
1
28
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
42
<-
<-
2
name
name
1
table.row
numeric
numeric
1
29
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
47
<-
<-
2
name
name
1
table.row
numeric
numeric
1
30
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
62
<-
<-
2
name
name
1
table.row
numeric
numeric
1
31
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
87
<-
<-
2
name
name
1
table.row
numeric
numeric
1
32
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
125
<-
<-
2
name
name
1
table.row
numeric
numeric
1
33
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
175
<-
<-
2
name
name
1
table.row
numeric
numeric
1
34
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
250
<-
<-
2
name
name
1
table.row
numeric
numeric
1
35
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
350
<-
<-
2
name
name
1
table.row
numeric
numeric
1
36
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
450
<-
<-
2
name
name
1
table.row
numeric
numeric
1
37
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
550
<-
<-
2
name
name
1
table.row
numeric
numeric
1
38
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
650
<-
<-
2
name
name
1
table.row
numeric
numeric
1
39
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
750
<-
<-
2
name
name
1
table.row
numeric
numeric
1
40
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
850
<-
<-
2
name
name
1
table.row
numeric
numeric
1
41
if
if
3
call
call
3
name
name
1
<=
name
name
1
n
numeric
numeric
1
950
<-
<-
2
name
name
1
table.row
numeric
numeric
1
42
<-
<-
2
name
name
1
table.row
numeric
numeric
1
43
if
if
3
call
call
3
name
name
1
==
name
name
1
alpha
numeric
numeric
1
0
{
{
4
call
call
2
name
name
1
cat
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
       Rao's Spacing Test of Uniformity
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Test Statistic =
call
call
3
name
name
1
round
name
name
1
U
numeric
numeric
1
5
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
U
call
call
4
name
name
1
[
name
name
1
rao.table
name
name
1
table.row
numeric
numeric
1
1
call
call
4
name
name
1
cat
character
character
1
P-value < 0.001
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
U
call
call
4
name
name
1
[
name
name
1
rao.table
name
name
1
table.row
numeric
numeric
1
2
call
call
4
name
name
1
cat
character
character
1
0.001 < P-value < 0.01
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
U
call
call
4
name
name
1
[
name
name
1
rao.table
name
name
1
table.row
numeric
numeric
1
3
call
call
4
name
name
1
cat
character
character
1
0.01 < P-value < 0.05
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
U
call
call
4
name
name
1
[
name
name
1
rao.table
name
name
1
table.row
numeric
numeric
1
4
call
call
4
name
name
1
cat
character
character
1
0.05 < P-value < 0.10
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
P-value > 0.10
character
character
1
\n
character
character
1
\n
{
{
8
call
call
2
name
name
1
cat
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
       Rao's Spacing Test of Uniformity
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Test Statistic =
call
call
3
name
name
1
round
name
name
1
U
numeric
numeric
1
5
character
character
1
\n
if
if
2
call
call
3
name
name
1
==
call
call
2
name
name
1
sum
call
call
3
name
name
1
==
name
name
1
alpha
call
call
5
name
name
1
c
numeric
numeric
1
0.001
numeric
numeric
1
0.01
numeric
numeric
1
0.050000000000000003
numeric
numeric
1
0.10000000000000001
numeric
numeric
1
0
call
call
2
name
name
1
stop
character
character
1
Invalid significance level
<-
<-
2
name
name
1
table.col
call
call
3
name
name
1
[
(
(
2
name
name
1
(
call
call
3
name
name
1
:
numeric
numeric
1
1
numeric
numeric
1
4
call
call
3
name
name
1
==
name
name
1
alpha
call
call
5
name
name
1
c
numeric
numeric
1
0.001
numeric
numeric
1
0.01
numeric
numeric
1
0.050000000000000003
numeric
numeric
1
0.10000000000000001
<-
<-
2
name
name
1
critical
call
call
4
name
name
1
[
name
name
1
rao.table
name
name
1
table.row
name
name
1
table.col
call
call
6
name
name
1
cat
character
character
1
Level
name
name
1
alpha
character
character
1
critical value =
name
name
1
critical
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
U
name
name
1
critical
call
call
4
name
name
1
cat
character
character
1
Reject null hypothesis of uniformity
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Do not reject null hypothesis of uniformity
character
character
1
\n
character
character
1
\n
character
character
3
data
alpha

rao.table
matrix
structure
3
numeric
numeric
172
247.31999999999999
245.19
236.81
229.46000000000001
224.41
219.52000000000001
215.44
211.87
208.69
205.87
203.33000000000001
201.03999999999999
198.96000000000001
197.05000000000001
195.28999999999999
193.66999999999999
192.16999999999999
190.78
189.47
188.25
187.11000000000001
186.03
185.00999999999999
184.05000000000001
183.13999999999999
182.28
181.44999999999999
177.88
174.99000000000001
172.58000000000001
170.53999999999999
163.59999999999999
159.44999999999999
154.50999999999999
151.56
148.06
145.96000000000001
144.53999999999999
143.47999999999999
142.66
142
141.44999999999999
140.99000000000001
221.13999999999999
211.93000000000001
206.78999999999999
202.55000000000001
198.46000000000001
195.27000000000001
192.37
189.88
187.66
185.68000000000001
183.90000000000001
182.28
180.81
179.46000000000001
178.22
177.08000000000001
176.00999999999999
175.02000000000001
174.09999999999999
173.22999999999999
172.41
171.63999999999999
170.91999999999999
170.22999999999999
169.58000000000001
168.96000000000001
168.38
165.81
163.72999999999999
162
160.53
155.49000000000001
152.46000000000001
148.84
146.66999999999999
144.09
142.53999999999999
141.47999999999999
140.69999999999999
140.09
139.59999999999999
139.19
138.84
186.44999999999999
183.44
180.65000000000001
177.83000000000001
175.68000000000001
173.68000000000001
171.97999999999999
170.44999999999999
169.09
167.87
166.75999999999999
165.75
164.83000000000001
163.97999999999999
163.19999999999999
162.47
161.78999999999999
161.16
160.56
160.00999999999999
159.47999999999999
158.99000000000001
158.52000000000001
158.06999999999999
157.65000000000001
157.25
156.87
155.19
153.81999999999999
152.68000000000001
151.69999999999999
148.34
146.28999999999999
143.83000000000001
142.34999999999999
140.56999999999999
139.5
138.77000000000001
138.22999999999999
137.80000000000001
137.46000000000001
137.18000000000001
136.94
168.02000000000001
168.66
166.30000000000001
165.05000000000001
163.56
162.36000000000001
161.22999999999999
160.24000000000001
159.33000000000001
158.5
157.75
157.06
156.43000000000001
155.84
155.28999999999999
154.78
154.31
153.86000000000001
153.44
153.05000000000001
152.68000000000001
152.31999999999999
151.99000000000001
151.66999999999999
151.37
151.08000000000001
150.80000000000001
149.59
148.59999999999999
147.75999999999999
147.05000000000001
144.56
143.03
141.18000000000001
140.06
138.71000000000001
137.88999999999999
137.33000000000001
136.91
136.59
136.33000000000001
136.11000000000001
135.91999999999999
integer
integer
2
43
4
NULL
NULL
0
rcard
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
4
<-
<-
2
name
name
1
i
numeric
numeric
1
1
<-
<-
2
name
name
1
result
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
while
while
2
call
call
3
name
name
1
<=
name
name
1
i
name
name
1
n
{
{
4
<-
<-
2
name
name
1
x
call
call
4
name
name
1
runif
numeric
numeric
1
1
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
y
call
call
4
name
name
1
runif
numeric
numeric
1
1
numeric
numeric
1
0
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
r
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
f
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
r
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
x
name
name
1
mu
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
if
if
2
call
call
3
name
name
1
<=
name
name
1
y
name
name
1
f
{
{
2
<-
<-
2
call
call
3
name
name
1
[
name
name
1
result
name
name
1
i
name
name
1
x
<-
<-
2
name
name
1
i
call
call
3
name
name
1
+
name
name
1
i
numeric
numeric
1
1
name
name
1
result
character
character
4
n
mu
r

rmixedvm
.named_I
structure
2
function
function
7
missing
missing
0
missing
missing
0
missing
missing
0
missing
missing
0
missing
missing
0
missing
missing
0
{
{
3
<-
<-
2
name
name
1
result
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
{
{
2
<-
<-
2
name
name
1
test
call
call
2
name
name
1
runif
numeric
numeric
1
1
if
if
3
call
call
3
name
name
1
<
name
name
1
test
name
name
1
p
<-
<-
2
call
call
3
name
name
1
[
name
name
1
result
name
name
1
i
call
call
4
name
name
1
rvm
numeric
numeric
1
1
name
name
1
mu1
name
name
1
kappa1
<-
<-
2
call
call
3
name
name
1
[
name
name
1
result
name
name
1
i
call
call
4
name
name
1
rvm
numeric
numeric
1
1
name
name
1
mu2
name
name
1
kappa2
character
character
3
i


name
name
1
result
character
character
7
n
mu1
mu2
kappa1
kappa2
p

rose.diag
.named_I
structure
2
function
function
15
missing
missing
0
missing
missing
0
integer
integer
1
1
character
character
1
angle
character
character
1
counter
name
name
1
T
numeric
numeric
1
1.5
integer
integer
1
8
character
character
1

integer
integer
1
16
name
name
1
F
integer
integer
1
1
integer
integer
1
40
integer
integer
1
1
{
{
16
if
if
2
call
call
3
name
name
1
==
name
name
1
type
character
character
1
angle
if
if
2
call
call
3
name
name
1
==
name
name
1
orientation
character
character
1
clock
call
call
2
name
name
1
stop
character
character
1
orientation must equal "counter" when type="angle"
<-
<-
2
name
name
1
data
call
call
3
name
name
1
%%
name
name
1
data
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
.named_I
structure
2
call
call
10
name
name
1
eqscplot
call
call
2
name
name
1
cos
.named_I
structure
2
call
call
4
name
name
1
seq
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
numeric
numeric
1
1000
character
character
4



length
call
call
2
name
name
1
sin
.named_I
structure
2
call
call
4
name
name
1
seq
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
numeric
numeric
1
1000
character
character
4



length
name
name
1
F
character
character
1

character
character
1

name
name
1
main
character
character
1
l
call
call
3
name
name
1
*
name
name
1
shrink
call
call
3
name
name
1
c
numeric
numeric
1
-1
numeric
numeric
1
1
call
call
3
name
name
1
*
name
name
1
shrink
call
call
3
name
name
1
c
numeric
numeric
1
-1
numeric
numeric
1
1
character
character
10



axes
xlab
ylab
main
type
xlim
ylim
if
if
2
call
call
3
name
name
1
==
name
name
1
label
name
name
1
T
{
{
1
if
if
3
call
call
3
name
name
1
==
name
name
1
type
character
character
1
direction
{
{
9
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0
call
call
3
name
name
1
c
numeric
numeric
1
0.90000000000000002
numeric
numeric
1
1
.named_I
structure
2
call
call
6
name
name
1
text
numeric
numeric
1
0
numeric
numeric
1
0.84999999999999998
character
character
1
N
name
name
1
labelex
numeric
numeric
1
0.5
character
character
6




cex
adj
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0
call
call
3
name
name
1
c
numeric
numeric
1
-0.90000000000000002
numeric
numeric
1
-1
.named_I
structure
2
call
call
6
name
name
1
text
numeric
numeric
1
0
numeric
numeric
1
-0.82499999999999996
character
character
1
S
name
name
1
labelex
numeric
numeric
1
0.5
character
character
6




cex
adj
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
-1
numeric
numeric
1
-0.90000000000000002
call
call
3
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0
.named_I
structure
2
call
call
5
name
name
1
text
numeric
numeric
1
-0.82499999999999996
numeric
numeric
1
0.0074999999999999997
character
character
1
W
name
name
1
labelex
character
character
5




cex
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0.90000000000000002
numeric
numeric
1
1
call
call
3
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0
.named_I
structure
2
call
call
5
name
name
1
text
numeric
numeric
1
0.84999999999999998
numeric
numeric
1
0.0074999999999999997
character
character
1
E
name
name
1
labelex
character
character
5




cex
.named_I
structure
2
call
call
6
name
name
1
text
numeric
numeric
1
0
numeric
numeric
1
0.014999999999999999
character
character
1
+
name
name
1
labelex
numeric
numeric
1
0.5
character
character
6




cex
adj
if
if
2
call
call
3
name
name
1
==
name
name
1
type
character
character
1
angle
{
{
8
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0
call
call
3
name
name
1
c
numeric
numeric
1
0.90000000000000002
numeric
numeric
1
1
.named_I
structure
2
call
call
6
name
name
1
text
numeric
numeric
1
0
numeric
numeric
1
0.84999999999999998
character
character
1
90
name
name
1
labelex
numeric
numeric
1
0.5
character
character
6




cex
adj
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0
call
call
3
name
name
1
c
numeric
numeric
1
-0.90000000000000002
numeric
numeric
1
-1
.named_I
structure
2
call
call
6
name
name
1
text
numeric
numeric
1
0
numeric
numeric
1
-0.82499999999999996
character
character
1
270
name
name
1
labelex
numeric
numeric
1
0.5
character
character
6




cex
adj
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
-1
numeric
numeric
1
-0.90000000000000002
call
call
3
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0
.named_I
structure
2
call
call
5
name
name
1
text
numeric
numeric
1
-0.80000000000000004
numeric
numeric
1
0.0074999999999999997
character
character
1
180
name
name
1
labelex
character
character
5




cex
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0.90000000000000002
numeric
numeric
1
1
call
call
3
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0
.named_I
structure
2
call
call
5
name
name
1
text
numeric
numeric
1
0.84999999999999998
numeric
numeric
1
0.0074999999999999997
character
character
1
0
name
name
1
labelex
character
character
5




cex
if
if
2
call
call
3
name
name
1
>
name
name
1
minor
integer
integer
1
0
{
{
2
<-
<-
2
name
name
1
increment
call
call
3
name
name
1
/
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
minor
integer
integer
1
1
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
minor
{
{
4
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
name
name
1
pi
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
name
name
1
pi
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
name
name
1
pi
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
name
name
1
pi
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
/
call
call
3
name
name
1
*
integer
integer
1
3
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
/
call
call
3
name
name
1
*
integer
integer
1
3
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
/
call
call
3
name
name
1
*
integer
integer
1
3
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
/
call
call
3
name
name
1
*
integer
integer
1
3
name
name
1
pi
integer
integer
1
2
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
cos
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
cos
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
3
name
name
1
c
call
call
3
name
name
1
*
numeric
numeric
1
0.94999999999999996
call
call
2
name
name
1
sin
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
call
call
2
name
name
1
sin
call
call
3
name
name
1
*
name
name
1
tick
name
name
1
increment
character
character
3
tick


<-
<-
2
name
name
1
count.data
call
call
3
name
name
1
-
name
name
1
data
call
call
3
name
name
1
$
name
name
1
.Machine
character
character
1
double.eps
<-
<-
2
name
name
1
bins.count
call
call
2
name
name
1
c
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
bins
<-
<-
2
name
name
1
mids
call
call
2
name
name
1
c
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
bins
<-
<-
2
name
name
1
sectors
call
call
2
name
name
1
c
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
bins
<-
<-
2
name
name
1
arc
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
name
name
1
bins
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
bins
{
{
1
if
if
3
call
call
3
name
name
1
==
name
name
1
i
integer
integer
1
1
{
{
3
<-
<-
2
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
call
call
2
name
name
1
sum
call
call
3
name
name
1
|
call
call
3
name
name
1
<=
name
name
1
count.data
call
call
3
name
name
1
/
name
name
1
arc
integer
integer
1
2
call
call
3
name
name
1
>
name
name
1
count.data
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
call
call
3
name
name
1
/
name
name
1
arc
integer
integer
1
2
<-
<-
2
call
call
3
name
name
1
[
name
name
1
mids
name
name
1
i
integer
integer
1
0
<-
<-
2
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
call
call
3
name
name
1
/
name
name
1
arc
integer
integer
1
2
{
{
3
<-
<-
2
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
call
call
2
name
name
1
sum
call
call
3
name
name
1
&
call
call
3
name
name
1
<=
name
name
1
count.data
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
call
call
3
name
name
1
/
integer
integer
1
1
integer
integer
1
2
name
name
1
arc
call
call
3
name
name
1
>
name
name
1
count.data
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
call
call
3
name
name
1
/
integer
integer
1
1
integer
integer
1
2
name
name
1
arc
<-
<-
2
call
call
3
name
name
1
[
name
name
1
mids
name
name
1
i
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
mids
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
name
name
1
arc
<-
<-
2
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
call
call
3
name
name
1
-
name
name
1
i
integer
integer
1
1
name
name
1
arc
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
character
character
3
i


if
if
2
call
call
3
name
name
1
==
name
name
1
type
character
character
1
direction
{
{
1
if
if
3
call
call
3
name
name
1
==
name
name
1
orientation
character
character
1
clock
{
{
2
<-
<-
2
name
name
1
mids
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
name
name
1
mids
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
<-
<-
2
name
name
1
sectors
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
name
name
1
sectors
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
{
{
2
<-
<-
2
name
name
1
mids
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
mids
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
<-
<-
2
name
name
1
sectors
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
sectors
call
call
3
name
name
1
/
name
name
1
pi
integer
integer
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
integer
integer
1
2
name
name
1
pi
if
if
2
call
call
3
name
name
1
==
name
name
1
stack
name
name
1
T
{
{
2
<-
<-
2
name
name
1
index
call
call
3
name
name
1
/
name
name
1
cex
name
name
1
dotsep
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
bins
{
{
1
if
if
2
call
call
3
name
name
1
!=
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
numeric
numeric
1
0
{
{
1
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
0
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
bins.count
name
name
1
i
numeric
numeric
1
1
{
{
4
<-
<-
2
name
name
1
r
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
*
name
name
1
j
name
name
1
index
<-
<-
2
name
name
1
x
call
call
3
name
name
1
*
name
name
1
r
call
call
2
name
name
1
cos
call
call
3
name
name
1
[
name
name
1
mids
name
name
1
i
<-
<-
2
name
name
1
y
call
call
3
name
name
1
*
name
name
1
r
call
call
2
name
name
1
sin
call
call
3
name
name
1
[
name
name
1
mids
name
name
1
i
.named_I
structure
2
call
call
5
name
name
1
points
name
name
1
x
name
name
1
y
name
name
1
cex
name
name
1
pch
character
character
5



cex
pch
character
character
3
j


character
character
3
i


<-
<-
2
name
name
1
radius
call
call
3
name
name
1
*
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
/
name
name
1
bins.count
name
name
1
n
name
name
1
prop
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
integer
integer
1
1
name
name
1
bins
{
{
1
if
if
2
call
call
3
name
name
1
!=
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
numeric
numeric
1
0
{
{
3
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
call
call
3
name
name
1
c
numeric
numeric
1
0
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
call
call
3
name
name
1
lines
call
call
3
name
name
1
c
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
cos
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
call
call
3
name
name
1
c
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
call
call
3
name
name
1
*
call
call
3
name
name
1
[
name
name
1
radius
name
name
1
i
call
call
2
name
name
1
sin
call
call
3
name
name
1
+
call
call
3
name
name
1
[
name
name
1
sectors
name
name
1
i
name
name
1
arc
character
character
3
i


character
character
15
data
bins
prop
type
orientation
label
labelex
minor
main
pch
stack
cex
dotsep
shrink

rtri
.named_I
structure
2
function
function
3
missing
missing
0
missing
missing
0
{
{
5
<-
<-
2
name
name
1
u
.named_I
structure
2
call
call
3
name
name
1
matrix
call
call
2
name
name
1
c
call
call
2
name
name
1
runif
name
name
1
n
numeric
numeric
1
1
character
character
3


ncol
<-
<-
2
name
name
1
get.theta
.named_I
structure
2
function
function
3
missing
missing
0
missing
missing
0
{
{
1
if
if
3
call
call
3
name
name
1
<
name
name
1
u
numeric
numeric
1
0.5
{
{
6
<-
<-
2
name
name
1
a
call
call
3
name
name
1
*
name
name
1
pi
name
name
1
r
<-
<-
2
name
name
1
b
call
call
2
name
name
1
-
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
4
call
call
3
name
name
1
*
call
call
3
name
name
1
^
name
name
1
pi
numeric
numeric
1
2
name
name
1
r
<-
<-
2
name
name
1
c
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
8
name
name
1
pi
name
name
1
u
<-
<-
2
name
name
1
theta1
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
2
name
name
1
-
name
name
1
b
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
-
call
call
3
name
name
1
^
name
name
1
b
numeric
numeric
1
2
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
4
name
name
1
a
name
name
1
c
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
a
<-
<-
2
name
name
1
theta2
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
2
name
name
1
-
name
name
1
b
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
-
call
call
3
name
name
1
^
name
name
1
b
numeric
numeric
1
2
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
4
name
name
1
a
name
name
1
c
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
a
call
call
3
name
name
1
min
name
name
1
theta1
name
name
1
theta2
{
{
6
<-
<-
2
name
name
1
a
call
call
3
name
name
1
*
name
name
1
pi
name
name
1
r
<-
<-
2
name
name
1
b
call
call
3
name
name
1
-
numeric
numeric
1
4
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
3
call
call
3
name
name
1
^
name
name
1
pi
numeric
numeric
1
2
name
name
1
r
<-
<-
2
name
name
1
c
call
call
3
name
name
1
-
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
call
call
3
name
name
1
^
name
name
1
pi
numeric
numeric
1
3
name
name
1
r
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
8
name
name
1
pi
name
name
1
u
<-
<-
2
name
name
1
theta1
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
2
name
name
1
-
name
name
1
b
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
-
call
call
3
name
name
1
^
name
name
1
b
numeric
numeric
1
2
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
4
name
name
1
a
name
name
1
c
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
a
<-
<-
2
name
name
1
theta2
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
2
name
name
1
-
name
name
1
b
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
-
call
call
3
name
name
1
^
name
name
1
b
numeric
numeric
1
2
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
4
name
name
1
a
name
name
1
c
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
a
call
call
3
name
name
1
max
name
name
1
theta1
name
name
1
theta2
character
character
3
u
r

<-
<-
2
name
name
1
theta
call
call
5
name
name
1
apply
name
name
1
u
numeric
numeric
1
1
name
name
1
get.theta
name
name
1
r
<-
<-
2
call
call
3
name
name
1
[
name
name
1
theta
call
call
3
name
name
1
>
name
name
1
theta
name
name
1
pi
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
theta
call
call
3
name
name
1
>
name
name
1
theta
name
name
1
pi
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
name
name
1
theta
character
character
3
n
r

rvm
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
7
<-
<-
2
name
name
1
vm
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
<-
<-
2
name
name
1
a
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
*
numeric
numeric
1
4
(
(
2
name
name
1
(
call
call
3
name
name
1
^
name
name
1
k
numeric
numeric
1
2
numeric
numeric
1
0.5
<-
<-
2
name
name
1
b
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
a
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
a
numeric
numeric
1
0.5
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
k
<-
<-
2
name
name
1
r
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
^
name
name
1
b
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
b
<-
<-
2
name
name
1
obs
numeric
numeric
1
1
while
while
2
call
call
3
name
name
1
<=
name
name
1
obs
name
name
1
n
{
{
6
<-
<-
2
name
name
1
U1
call
call
4
name
name
1
runif
numeric
numeric
1
1
numeric
numeric
1
0
numeric
numeric
1
1
<-
<-
2
name
name
1
z
call
call
2
name
name
1
cos
call
call
3
name
name
1
*
name
name
1
pi
name
name
1
U1
<-
<-
2
name
name
1
f
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
*
name
name
1
r
name
name
1
z
(
(
2
name
name
1
(
call
call
3
name
name
1
+
name
name
1
r
name
name
1
z
<-
<-
2
name
name
1
c
call
call
3
name
name
1
*
name
name
1
k
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
r
name
name
1
f
<-
<-
2
name
name
1
U2
call
call
4
name
name
1
runif
numeric
numeric
1
1
numeric
numeric
1
0
numeric
numeric
1
1
if
if
3
call
call
3
name
name
1
>
call
call
3
name
name
1
-
call
call
3
name
name
1
*
name
name
1
c
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
2
name
name
1
c
name
name
1
U2
numeric
numeric
1
0
{
{
4
<-
<-
2
name
name
1
U3
call
call
4
name
name
1
runif
numeric
numeric
1
1
numeric
numeric
1
0
numeric
numeric
1
1
<-
<-
2
call
call
3
name
name
1
[
name
name
1
vm
name
name
1
obs
call
call
3
name
name
1
+
call
call
3
name
name
1
*
call
call
2
name
name
1
sign
call
call
3
name
name
1
-
name
name
1
U3
numeric
numeric
1
0.5
call
call
2
name
name
1
acos
name
name
1
f
name
name
1
mean
<-
<-
2
call
call
3
name
name
1
[
name
name
1
vm
name
name
1
obs
call
call
3
name
name
1
%%
call
call
3
name
name
1
[
name
name
1
vm
name
name
1
obs
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
obs
call
call
3
name
name
1
+
name
name
1
obs
numeric
numeric
1
1
{
{
1
if
if
2
call
call
3
name
name
1
>=
call
call
3
name
name
1
-
call
call
3
name
name
1
+
call
call
2
name
name
1
log
call
call
3
name
name
1
/
name
name
1
c
name
name
1
U2
numeric
numeric
1
1
name
name
1
c
numeric
numeric
1
0
{
{
4
<-
<-
2
name
name
1
U3
call
call
4
name
name
1
runif
numeric
numeric
1
1
numeric
numeric
1
0
numeric
numeric
1
1
<-
<-
2
call
call
3
name
name
1
[
name
name
1
vm
name
name
1
obs
call
call
3
name
name
1
+
call
call
3
name
name
1
*
call
call
2
name
name
1
sign
call
call
3
name
name
1
-
name
name
1
U3
numeric
numeric
1
0.5
call
call
2
name
name
1
acos
name
name
1
f
name
name
1
mean
<-
<-
2
call
call
3
name
name
1
[
name
name
1
vm
name
name
1
obs
call
call
3
name
name
1
%%
call
call
3
name
name
1
[
name
name
1
vm
name
name
1
obs
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
obs
call
call
3
name
name
1
+
name
name
1
obs
numeric
numeric
1
1
name
name
1
vm
character
character
4
n
mean
k

rwrpcauchy
.named_I
structure
2
function
function
4
missing
missing
0
numeric
numeric
1
0
call
call
2
name
name
1
exp
numeric
numeric
1
-1
{
{
2
if
if
3
call
call
3
name
name
1
==
name
name
1
rho
numeric
numeric
1
0
<-
<-
2
name
name
1
result
call
call
4
name
name
1
runif
name
name
1
n
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
if
if
3
call
call
3
name
name
1
==
name
name
1
rho
numeric
numeric
1
1
<-
<-
2
name
name
1
result
call
call
3
name
name
1
rep
name
name
1
location
name
name
1
n
{
{
2
<-
<-
2
name
name
1
scale
call
call
2
name
name
1
-
call
call
2
name
name
1
log
name
name
1
rho
<-
<-
2
name
name
1
result
call
call
3
name
name
1
%%
call
call
4
name
name
1
rcauchy
name
name
1
n
name
name
1
location
name
name
1
scale
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
name
name
1
result
character
character
4
n
location
rho

rwrpnorm
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
2
if
if
3
call
call
3
name
name
1
==
name
name
1
rho
numeric
numeric
1
0
<-
<-
2
name
name
1
result
call
call
4
name
name
1
runif
name
name
1
n
numeric
numeric
1
0
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
if
if
3
call
call
3
name
name
1
==
name
name
1
rho
numeric
numeric
1
1
<-
<-
2
name
name
1
result
call
call
3
name
name
1
rep
name
name
1
mu
name
name
1
n
{
{
2
<-
<-
2
name
name
1
sd
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
*
numeric
numeric
1
-2
call
call
2
name
name
1
log
name
name
1
rho
<-
<-
2
name
name
1
result
call
call
3
name
name
1
%%
call
call
4
name
name
1
rnorm
name
name
1
n
name
name
1
mu
name
name
1
sd
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
name
name
1
result
character
character
4
n
mu
rho

rwrpstab
.named_I
structure
2
function
function
4
missing
missing
0
missing
missing
0
missing
missing
0
{
{
1
call
call
3
name
name
1
%%
call
call
4
name
name
1
rstab
name
name
1
n
name
name
1
index
name
name
1
skewness
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
character
character
4
n
index
skewness

trig.moment
.named_I
structure
2
function
function
4
missing
missing
0
numeric
numeric
1
1
name
name
1
F
{
{
9
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
sinr
call
call
2
name
name
1
sum
call
call
2
name
name
1
sin
name
name
1
data
<-
<-
2
name
name
1
cosr
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
name
name
1
data
<-
<-
2
name
name
1
circmean
call
call
3
name
name
1
atan
name
name
1
sinr
name
name
1
cosr
<-
<-
2
name
name
1
sin.p
call
call
3
name
name
1
/
call
call
2
name
name
1
sum
call
call
2
name
name
1
sin
call
call
3
name
name
1
*
name
name
1
p
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
data
call
call
3
name
name
1
*
name
name
1
circmean
name
name
1
center
name
name
1
n
<-
<-
2
name
name
1
cos.p
call
call
3
name
name
1
/
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
call
call
3
name
name
1
*
name
name
1
p
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
data
call
call
3
name
name
1
*
name
name
1
circmean
name
name
1
center
name
name
1
n
<-
<-
2
name
name
1
mu.p
call
call
3
name
name
1
atan
name
name
1
sin.p
name
name
1
cos.p
<-
<-
2
name
name
1
rho.p
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
sin.p
numeric
numeric
1
2
call
call
3
name
name
1
^
name
name
1
cos.p
numeric
numeric
1
2
call
call
5
name
name
1
data.frame
name
name
1
mu.p
name
name
1
rho.p
name
name
1
cos.p
name
name
1
sin.p
character
character
4
data
p
center

v0.test
.named_I
structure
2
function
function
3
missing
missing
0
numeric
numeric
1
0
{
{
8
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
r0.bar
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
2
name
name
1
sum
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
data
name
name
1
mu0
name
name
1
n
<-
<-
2
name
name
1
z0
call
call
3
name
name
1
*
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
n
name
name
1
r0.bar
<-
<-
2
name
name
1
pz
call
call
2
name
name
1
pnorm
name
name
1
z0
<-
<-
2
name
name
1
fz
call
call
2
name
name
1
dnorm
name
name
1
z0
<-
<-
2
name
name
1
p.value
call
call
3
name
name
1
+
call
call
3
name
name
1
-
numeric
numeric
1
1
name
name
1
pz
call
call
3
name
name
1
*
name
name
1
fz
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
*
numeric
numeric
1
3
name
name
1
z0
call
call
3
name
name
1
^
name
name
1
z0
numeric
numeric
1
3
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
16
name
name
1
n
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
+
call
call
3
name
name
1
*
numeric
numeric
1
15
name
name
1
z0
call
call
3
name
name
1
*
numeric
numeric
1
305
call
call
3
name
name
1
^
name
name
1
z0
numeric
numeric
1
3
call
call
3
name
name
1
*
numeric
numeric
1
125
call
call
3
name
name
1
^
name
name
1
z0
numeric
numeric
1
5
call
call
3
name
name
1
*
numeric
numeric
1
9
call
call
3
name
name
1
^
name
name
1
z0
numeric
numeric
1
7
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
4608
call
call
3
name
name
1
^
name
name
1
n
numeric
numeric
1
2
<-
<-
2
name
name
1
result
.named_I
structure
2
call
call
3
name
name
1
list
name
name
1
r0.bar
name
name
1
p.value
character
character
3

r0.bar
p.value
name
name
1
result
character
character
3
data
mu0

vm.bootstrap.ci
.named_I
structure
2
function
function
6
missing
missing
0
name
name
1
F
numeric
numeric
1
0.050000000000000003
numeric
numeric
1
1000
name
name
1
T
{
{
22
<-
<-
2
name
name
1
old.seed
name
name
1
.Random.seed
<-
<-
2
name
name
1
mean.bs
.named_I
structure
2
call
call
5
name
name
1
bootstrap
name
name
1
data
name
name
1
circ.mean
name
name
1
reps
name
name
1
print
character
character
5



B
trace
<-
<-
2
name
name
1
mean.reps
call
call
3
name
name
1
$
name
name
1
mean.bs
character
character
1
replicates
<-
<-
2
name
name
1
mean.reps
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
mean.reps
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
B
call
call
3
name
name
1
$
name
name
1
mean.bs
character
character
1
B
<-
<-
2
name
name
1
spacings
call
call
3
name
name
1
c
call
call
2
name
name
1
diff
name
name
1
mean.reps
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
[
name
name
1
mean.reps
numeric
numeric
1
1
call
call
3
name
name
1
[
name
name
1
mean.reps
name
name
1
B
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
max.spacing
call
call
3
name
name
1
[
(
(
2
name
name
1
(
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
B
call
call
3
name
name
1
==
name
name
1
spacings
call
call
2
name
name
1
max
name
name
1
spacings
<-
<-
2
name
name
1
off.set
call
call
3
name
name
1
-
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
3
name
name
1
[
name
name
1
mean.reps
call
call
3
name
name
1
+
name
name
1
max.spacing
numeric
numeric
1
1
if
if
3
call
call
3
name
name
1
!=
name
name
1
max.spacing
name
name
1
B
<-
<-
2
name
name
1
mean.reps2
call
call
3
name
name
1
+
name
name
1
mean.reps
name
name
1
off.set
<-
<-
2
name
name
1
mean.reps2
name
name
1
mean.reps
<-
<-
2
name
name
1
mean.reps2
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
mean.reps2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
mean.ci
call
call
3
name
name
1
quantile
name
name
1
mean.reps2
call
call
3
name
name
1
c
call
call
3
name
name
1
/
name
name
1
alpha
numeric
numeric
1
2
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
/
name
name
1
alpha
numeric
numeric
1
2
if
if
2
call
call
3
name
name
1
!=
name
name
1
max.spacing
name
name
1
B
<-
<-
2
name
name
1
mean.ci
call
call
3
name
name
1
-
name
name
1
mean.ci
name
name
1
off.set
<-
<-
2
name
name
1
kappa.bs
.named_I
structure
2
call
call
7
name
name
1
bootstrap
name
name
1
data
name
name
1
est.kappa
name
name
1
reps
call
call
2
name
name
1
list
name
name
1
bias
name
name
1
old.seed
name
name
1
print
character
character
7



B
args.stat
seed
trace
<-
<-
2
name
name
1
kappa.reps
call
call
3
name
name
1
$
name
name
1
kappa.bs
character
character
1
replicates
<-
<-
2
name
name
1
kappa.ci
call
call
3
name
name
1
quantile
name
name
1
kappa.reps
call
call
3
name
name
1
c
call
call
3
name
name
1
/
name
name
1
alpha
numeric
numeric
1
2
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
/
name
name
1
alpha
numeric
numeric
1
2
<-
<-
2
name
name
1
result
call
call
5
name
name
1
list
name
name
1
mean.ci
name
name
1
kappa.ci
call
call
2
name
name
1
c
name
name
1
mean.reps
call
call
2
name
name
1
c
name
name
1
kappa.reps
<-
<-
2
call
call
2
name
name
1
names
name
name
1
result
call
call
5
name
name
1
c
character
character
1
mu.ci
character
character
1
kappa.ci
character
character
1
mu.reps
character
character
1
kappa.reps
call
call
3
name
name
1
cat
character
character
1
Bootstrap Confidence Intervals for Mean Direction and Concentration
character
character
1
\n
call
call
5
name
name
1
cat
character
character
1
Confidence Level:  
call
call
3
name
name
1
*
numeric
numeric
1
100
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
1
name
name
1
alpha
character
character
1
%
character
character
1
\n
call
call
7
name
name
1
cat
character
character
1
Mean Direction:           
character
character
1
Low =
call
call
3
name
name
1
round
call
call
3
name
name
1
[
name
name
1
mean.ci
numeric
numeric
1
1
numeric
numeric
1
2
character
character
1
  High =
call
call
3
name
name
1
round
call
call
3
name
name
1
[
name
name
1
mean.ci
numeric
numeric
1
2
numeric
numeric
1
2
character
character
1
\n
call
call
7
name
name
1
cat
character
character
1
Concentration Parameter:  
character
character
1
Low =
call
call
3
name
name
1
round
call
call
3
name
name
1
[
name
name
1
kappa.ci
numeric
numeric
1
1
numeric
numeric
1
2
character
character
1
  High =
call
call
3
name
name
1
round
call
call
3
name
name
1
[
name
name
1
kappa.ci
numeric
numeric
1
2
numeric
numeric
1
2
character
character
1
\n
name
name
1
result
character
character
6
data
bias
alpha
reps
print

vm.ml
.named_I
structure
2
function
function
3
missing
missing
0
name
name
1
F
{
{
4
<-
<-
2
name
name
1
mu
call
call
2
name
name
1
circ.mean
name
name
1
data
<-
<-
2
name
name
1
kappa
call
call
2
name
name
1
A1inv
call
call
2
name
name
1
mean
call
call
2
name
name
1
cos
call
call
3
name
name
1
-
name
name
1
data
name
name
1
mu
if
if
2
call
call
3
name
name
1
==
name
name
1
bias
name
name
1
T
{
{
4
<-
<-
2
name
name
1
kappa.ml
name
name
1
kappa
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
if
if
2
call
call
3
name
name
1
<
name
name
1
kappa.ml
numeric
numeric
1
2
<-
<-
2
name
name
1
kappa
call
call
3
name
name
1
max
call
call
3
name
name
1
-
name
name
1
kappa.ml
call
call
3
name
name
1
*
numeric
numeric
1
2
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
*
name
name
1
n
name
name
1
kappa.ml
numeric
numeric
1
-1
numeric
numeric
1
0
if
if
2
call
call
3
name
name
1
>=
name
name
1
kappa.ml
numeric
numeric
1
2
<-
<-
2
name
name
1
kappa
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
n
numeric
numeric
1
1
numeric
numeric
1
3
name
name
1
kappa.ml
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
n
numeric
numeric
1
3
name
name
1
n
call
call
3
name
name
1
data.frame
name
name
1
mu
name
name
1
kappa
character
character
3
data
bias

watson
.named_I
structure
2
function
function
4
missing
missing
0
numeric
numeric
1
0
character
character
1
uniform
{
{
1
if
if
3
call
call
3
name
name
1
==
name
name
1
dist
character
character
1
uniform
{
{
13
call
call
5
name
name
1
cat
character
character
1
\n
character
character
1
      Watson's Test for Circular Uniformity
character
character
1
\n
character
character
1
\n
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
u
call
call
3
name
name
1
/
call
call
2
name
name
1
sort
name
name
1
data
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
u.bar
call
call
2
name
name
1
mean
name
name
1
u
<-
<-
2
name
name
1
i
call
call
2
name
name
1
seq
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
<-
<-
2
name
name
1
sum.terms
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
3
name
name
1
-
name
name
1
u
name
name
1
u.bar
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
i
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
n
numeric
numeric
1
0.5
numeric
numeric
1
2
<-
<-
2
name
name
1
u2
call
call
3
name
name
1
+
call
call
2
name
name
1
sum
name
name
1
sum.terms
call
call
3
name
name
1
/
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
12
name
name
1
n
<-
<-
2
name
name
1
u2
call
call
3
name
name
1
*
(
(
2
name
name
1
(
call
call
3
name
name
1
+
call
call
3
name
name
1
-
name
name
1
u2
call
call
3
name
name
1
/
numeric
numeric
1
0.10000000000000001
name
name
1
n
call
call
3
name
name
1
/
numeric
numeric
1
0.10000000000000001
(
(
2
name
name
1
(
call
call
3
name
name
1
^
name
name
1
n
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
/
numeric
numeric
1
0.80000000000000004
name
name
1
n
<-
<-
2
name
name
1
crits
call
call
6
name
name
1
c
numeric
numeric
1
99
numeric
numeric
1
0.26700000000000002
numeric
numeric
1
0.221
numeric
numeric
1
0.187
numeric
numeric
1
0.152
if
if
2
call
call
3
name
name
1
<
name
name
1
n
numeric
numeric
1
8
{
{
1
call
call
4
name
name
1
cat
character
character
1
Total Sample Size < 8:  Results are not valid
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
Test Statistic:
call
call
3
name
name
1
round
name
name
1
u2
numeric
numeric
1
4
character
character
1
\n
if
if
2
call
call
3
name
name
1
==
call
call
2
name
name
1
sum
call
call
3
name
name
1
==
name
name
1
alpha
call
call
6
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0.01
numeric
numeric
1
0.025000000000000001
numeric
numeric
1
0.050000000000000003
numeric
numeric
1
0.10000000000000001
numeric
numeric
1
0
call
call
2
name
name
1
stop
character
character
1
Invalid input for alpha
if
if
3
call
call
3
name
name
1
==
name
name
1
alpha
numeric
numeric
1
0
{
{
1
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
numeric
numeric
1
0.26700000000000002
call
call
4
name
name
1
cat
character
character
1
P-value < 0.01
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
numeric
numeric
1
0.221
call
call
4
name
name
1
cat
character
character
1
0.01 < P-value < 0.025
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
numeric
numeric
1
0.187
call
call
4
name
name
1
cat
character
character
1
0.025 < P-value < 0.05
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
numeric
numeric
1
0.152
call
call
4
name
name
1
cat
character
character
1
0.05 < P-value < 0.10
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
P-value > 0.10
character
character
1
\n
character
character
1
\n
{
{
5
<-
<-
2
name
name
1
index
call
call
3
name
name
1
[
(
(
2
name
name
1
(
call
call
3
name
name
1
:
numeric
numeric
1
1
numeric
numeric
1
5
call
call
3
name
name
1
==
name
name
1
alpha
call
call
6
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0.01
numeric
numeric
1
0.025000000000000001
numeric
numeric
1
0.050000000000000003
numeric
numeric
1
0.10000000000000001
<-
<-
2
name
name
1
Critical
call
call
3
name
name
1
[
name
name
1
crits
name
name
1
index
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
name
name
1
Critical
<-
<-
2
name
name
1
Reject
character
character
1
Reject Null Hypothesis
<-
<-
2
name
name
1
Reject
character
character
1
Do Not Reject Null Hypothesis
call
call
6
name
name
1
cat
character
character
1
Level
name
name
1
alpha
character
character
1
Critical Value:
call
call
3
name
name
1
round
name
name
1
Critical
numeric
numeric
1
4
character
character
1
\n
call
call
4
name
name
1
cat
name
name
1
Reject
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
==
name
name
1
dist
character
character
1
vm
{
{
15
call
call
5
name
name
1
cat
character
character
1
\n
character
character
1
      Watson's Test for the von Mises Distribution
character
character
1
\n
character
character
1
\n
<-
<-
2
name
name
1
u2.crits
call
call
5
name
name
1
cbind
call
call
8
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0.5
numeric
numeric
1
1
numeric
numeric
1
1.5
numeric
numeric
1
2
numeric
numeric
1
4
numeric
numeric
1
100
call
call
8
name
name
1
c
numeric
numeric
1
0.051999999999999998
numeric
numeric
1
0.056000000000000001
numeric
numeric
1
0.066000000000000003
numeric
numeric
1
0.076999999999999999
numeric
numeric
1
0.084000000000000005
numeric
numeric
1
0.092999999999999999
numeric
numeric
1
0.096000000000000002
call
call
8
name
name
1
c
numeric
numeric
1
0.060999999999999999
numeric
numeric
1
0.066000000000000003
numeric
numeric
1
0.079000000000000001
numeric
numeric
1
0.091999999999999998
numeric
numeric
1
0.10100000000000001
numeric
numeric
1
0.113
numeric
numeric
1
0.11700000000000001
call
call
8
name
name
1
c
numeric
numeric
1
0.081000000000000003
numeric
numeric
1
0.089999999999999997
numeric
numeric
1
0.11
numeric
numeric
1
0.128
numeric
numeric
1
0.14199999999999999
numeric
numeric
1
0.158
numeric
numeric
1
0.16400000000000001
<-
<-
2
name
name
1
n
call
call
2
name
name
1
length
name
name
1
data
<-
<-
2
name
name
1
mu.hat
call
call
2
name
name
1
circ.mean
name
name
1
data
<-
<-
2
name
name
1
kappa.hat
call
call
2
name
name
1
est.kappa
name
name
1
data
<-
<-
2
name
name
1
data
call
call
3
name
name
1
%%
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
data
name
name
1
mu.hat
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
data
.named_I
structure
2
call
call
3
name
name
1
matrix
name
name
1
data
numeric
numeric
1
1
character
character
3


ncol
<-
<-
2
name
name
1
z
call
call
6
name
name
1
apply
name
name
1
data
numeric
numeric
1
1
name
name
1
pvm
numeric
numeric
1
0
name
name
1
kappa.hat
<-
<-
2
name
name
1
z
call
call
2
name
name
1
sort
name
name
1
z
<-
<-
2
name
name
1
z.bar
call
call
2
name
name
1
mean
name
name
1
z
<-
<-
2
name
name
1
i
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
<-
<-
2
name
name
1
sum.terms
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
z
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
i
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
n
numeric
numeric
1
2
<-
<-
2
name
name
1
Value
call
call
3
name
name
1
+
call
call
3
name
name
1
-
call
call
2
name
name
1
sum
name
name
1
sum.terms
call
call
3
name
name
1
*
name
name
1
n
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
z.bar
numeric
numeric
1
0.5
numeric
numeric
1
2
call
call
3
name
name
1
/
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
12
name
name
1
n
if
if
3
call
call
3
name
name
1
<
name
name
1
kappa.hat
numeric
numeric
1
0.25
<-
<-
2
name
name
1
row
numeric
numeric
1
1
if
if
3
call
call
3
name
name
1
<
name
name
1
kappa.hat
numeric
numeric
1
0.75
<-
<-
2
name
name
1
row
numeric
numeric
1
2
if
if
3
call
call
3
name
name
1
<
name
name
1
kappa.hat
numeric
numeric
1
1.25
<-
<-
2
name
name
1
row
numeric
numeric
1
3
if
if
3
call
call
3
name
name
1
<
name
name
1
kappa.hat
numeric
numeric
1
1.75
<-
<-
2
name
name
1
row
numeric
numeric
1
4
if
if
3
call
call
3
name
name
1
<
name
name
1
kappa.hat
numeric
numeric
1
3
<-
<-
2
name
name
1
row
numeric
numeric
1
5
if
if
3
call
call
3
name
name
1
<
name
name
1
kappa.hat
numeric
numeric
1
5
<-
<-
2
name
name
1
row
numeric
numeric
1
6
<-
<-
2
name
name
1
row
numeric
numeric
1
7
if
if
3
call
call
3
name
name
1
!=
name
name
1
alpha
numeric
numeric
1
0
{
{
6
if
if
3
call
call
3
name
name
1
==
name
name
1
alpha
numeric
numeric
1
0.10000000000000001
<-
<-
2
name
name
1
col
numeric
numeric
1
2
if
if
3
call
call
3
name
name
1
==
name
name
1
alpha
numeric
numeric
1
0.050000000000000003
<-
<-
2
name
name
1
col
numeric
numeric
1
3
if
if
3
call
call
3
name
name
1
==
name
name
1
alpha
numeric
numeric
1
0.01
<-
<-
2
name
name
1
col
numeric
numeric
1
4
{
{
2
call
call
4
name
name
1
cat
character
character
1
Invalid input for alpha
character
character
1
\n
character
character
1
\n
break
break
0
<-
<-
2
name
name
1
Critical
call
call
4
name
name
1
[
name
name
1
u2.crits
name
name
1
row
name
name
1
col
if
if
3
call
call
3
name
name
1
>
name
name
1
Value
name
name
1
Critical
<-
<-
2
name
name
1
Reject
character
character
1
Reject Null Hypothesis
<-
<-
2
name
name
1
Reject
character
character
1
Do Not Reject Null Hypothesis
call
call
4
name
name
1
cat
character
character
1
Test Statistic:
call
call
3
name
name
1
round
name
name
1
Value
numeric
numeric
1
4
character
character
1
\n
call
call
6
name
name
1
cat
character
character
1
Level
name
name
1
alpha
character
character
1
Critical Value:
call
call
3
name
name
1
round
name
name
1
Critical
numeric
numeric
1
4
character
character
1
\n
call
call
4
name
name
1
cat
name
name
1
Reject
character
character
1
\n
character
character
1
\n
{
{
2
call
call
4
name
name
1
cat
character
character
1
Test Statistic:
call
call
3
name
name
1
round
name
name
1
Value
numeric
numeric
1
4
character
character
1
\n
if
if
3
call
call
3
name
name
1
<
name
name
1
Value
call
call
4
name
name
1
[
name
name
1
u2.crits
name
name
1
row
numeric
numeric
1
2
call
call
4
name
name
1
cat
character
character
1
P-value > 0.10
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
&&
(
(
2
name
name
1
(
call
call
3
name
name
1
>=
name
name
1
Value
call
call
4
name
name
1
[
name
name
1
u2.crits
name
name
1
row
numeric
numeric
1
2
(
(
2
name
name
1
(
call
call
3
name
name
1
<
name
name
1
Value
call
call
4
name
name
1
[
name
name
1
u2.crits
name
name
1
row
numeric
numeric
1
3
call
call
4
name
name
1
cat
character
character
1
0.05 < P-value > 0.10
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
&&
(
(
2
name
name
1
(
call
call
3
name
name
1
>=
name
name
1
Value
call
call
4
name
name
1
[
name
name
1
u2.crits
name
name
1
row
numeric
numeric
1
3
(
(
2
name
name
1
(
call
call
3
name
name
1
<
name
name
1
Value
call
call
4
name
name
1
[
name
name
1
u2.crits
name
name
1
row
numeric
numeric
1
4
call
call
4
name
name
1
cat
character
character
1
0.01 < P-value > 0.05
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
P-value < 0.01
character
character
1
\n
character
character
1
\n
call
call
2
name
name
1
stop
character
character
1
Distribution must be either uniform or von Mises
character
character
4
data
alpha
dist

watson.two
.named_I
structure
2
function
function
5
missing
missing
0
missing
missing
0
numeric
numeric
1
0
name
name
1
F
{
{
20
<-
<-
2
name
name
1
n1
call
call
2
name
name
1
length
name
name
1
data1
<-
<-
2
name
name
1
n2
call
call
2
name
name
1
length
name
name
1
data2
<-
<-
2
name
name
1
n
call
call
3
name
name
1
+
name
name
1
n1
name
name
1
n2
if
if
2
call
call
3
name
name
1
<
name
name
1
n
numeric
numeric
1
18
call
call
4
name
name
1
cat
character
character
1
Total Sample Size < 18:  Consult tabulated critical values
character
character
1
\n
character
character
1
\n
if
if
2
call
call
3
name
name
1
==
name
name
1
plot
name
name
1
T
{
{
5
<-
<-
2
name
name
1
data1
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
data1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
<-
<-
2
name
name
1
data2
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
data2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
.named_I
structure
2
call
call
5
name
name
1
plot.edf
name
name
1
data1
character
character
1
Comparison of Empirical CDFs
character
character
1

character
character
1

character
character
5


main
xlab
ylab
.named_I
structure
2
call
call
2
name
name
1
par
name
name
1
T
character
character
2

new
.named_I
structure
2
call
call
6
name
name
1
plot.edf
name
name
1
data2
character
character
1

character
character
1

name
name
1
F
numeric
numeric
1
2
character
character
6


xlab
ylab
axes
lty
call
call
5
name
name
1
cat
character
character
1
\n
character
character
1
      Watson's Two-Sample Test of Homogeneity
character
character
1
\n
character
character
1
\n
<-
<-
2
name
name
1
data1
call
call
3
name
name
1
cbind
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
data1
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
3
name
name
1
rep
numeric
numeric
1
1
name
name
1
n1
<-
<-
2
name
name
1
data2
call
call
3
name
name
1
cbind
call
call
2
name
name
1
sort
call
call
3
name
name
1
%%
name
name
1
data2
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
3
name
name
1
rep
numeric
numeric
1
2
name
name
1
n2
<-
<-
2
name
name
1
data
call
call
3
name
name
1
rbind
name
name
1
data1
name
name
1
data2
<-
<-
2
name
name
1
rank
call
call
2
name
name
1
order
call
call
4
name
name
1
[
name
name
1
data
missing
missing
0
numeric
numeric
1
1
<-
<-
2
name
name
1
data
call
call
3
name
name
1
cbind
call
call
4
name
name
1
[
name
name
1
data
name
name
1
rank
missing
missing
0
call
call
2
name
name
1
seq
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
<-
<-
2
name
name
1
a
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
<-
<-
2
name
name
1
b
call
call
2
name
name
1
c
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
.named_I
structure
2
for
for
3
NULL
NULL
0
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
n
{
{
2
<-
<-
2
call
call
3
name
name
1
[
name
name
1
a
name
name
1
i
call
call
2
name
name
1
sum
call
call
3
name
name
1
==
call
call
4
name
name
1
[
name
name
1
data
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
i
numeric
numeric
1
2
numeric
numeric
1
1
<-
<-
2
call
call
3
name
name
1
[
name
name
1
b
name
name
1
i
call
call
2
name
name
1
sum
call
call
3
name
name
1
==
call
call
4
name
name
1
[
name
name
1
data
call
call
3
name
name
1
:
numeric
numeric
1
1
name
name
1
i
numeric
numeric
1
2
numeric
numeric
1
2
character
character
3
i


<-
<-
2
name
name
1
d
call
call
3
name
name
1
-
call
call
3
name
name
1
/
name
name
1
b
name
name
1
n2
call
call
3
name
name
1
/
name
name
1
a
name
name
1
n1
<-
<-
2
name
name
1
dbar
call
call
2
name
name
1
mean
name
name
1
d
<-
<-
2
name
name
1
u2
call
call
3
name
name
1
*
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
name
name
1
n1
name
name
1
n2
call
call
3
name
name
1
^
name
name
1
n
numeric
numeric
1
2
call
call
2
name
name
1
sum
call
call
3
name
name
1
^
(
(
2
name
name
1
(
call
call
3
name
name
1
-
name
name
1
d
name
name
1
dbar
numeric
numeric
1
2
<-
<-
2
name
name
1
crits
call
call
6
name
name
1
c
numeric
numeric
1
99
numeric
numeric
1
0.38500000000000001
numeric
numeric
1
0.26800000000000002
numeric
numeric
1
0.187
numeric
numeric
1
0.152
call
call
4
name
name
1
cat
character
character
1
Test Statistic:
call
call
3
name
name
1
round
name
name
1
u2
numeric
numeric
1
4
character
character
1
\n
if
if
3
call
call
3
name
name
1
==
call
call
2
name
name
1
sum
call
call
3
name
name
1
==
name
name
1
alpha
call
call
6
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0.001
numeric
numeric
1
0.01
numeric
numeric
1
0.050000000000000003
numeric
numeric
1
0.10000000000000001
numeric
numeric
1
0
call
call
2
name
name
1
stop
character
character
1
Invalid input for alpha
if
if
3
call
call
3
name
name
1
==
name
name
1
alpha
numeric
numeric
1
0
{
{
1
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
numeric
numeric
1
0.38500000000000001
call
call
4
name
name
1
cat
character
character
1
P-value < 0.001
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
numeric
numeric
1
0.26800000000000002
call
call
4
name
name
1
cat
character
character
1
0.001 < P-value < 0.01
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
numeric
numeric
1
0.187
call
call
4
name
name
1
cat
character
character
1
0.01 < P-value < 0.05
character
character
1
\n
character
character
1
\n
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
numeric
numeric
1
0.152
call
call
4
name
name
1
cat
character
character
1
0.05 < P-value < 0.10
character
character
1
\n
character
character
1
\n
call
call
4
name
name
1
cat
character
character
1
P-value > 0.10
character
character
1
\n
character
character
1
\n
{
{
5
<-
<-
2
name
name
1
index
call
call
3
name
name
1
[
(
(
2
name
name
1
(
call
call
3
name
name
1
:
numeric
numeric
1
1
numeric
numeric
1
5
call
call
3
name
name
1
==
name
name
1
alpha
call
call
6
name
name
1
c
numeric
numeric
1
0
numeric
numeric
1
0.001
numeric
numeric
1
0.01
numeric
numeric
1
0.050000000000000003
numeric
numeric
1
0.10000000000000001
<-
<-
2
name
name
1
Critical
call
call
3
name
name
1
[
name
name
1
crits
name
name
1
index
if
if
3
call
call
3
name
name
1
>
name
name
1
u2
name
name
1
Critical
<-
<-
2
name
name
1
Reject
character
character
1
Reject Null Hypothesis
<-
<-
2
name
name
1
Reject
character
character
1
Do Not Reject Null Hypothesis
call
call
6
name
name
1
cat
character
character
1
Level
name
name
1
alpha
character
character
1
Critical Value:
call
call
3
name
name
1
round
name
name
1
Critical
numeric
numeric
1
4
character
character
1
\n
call
call
4
name
name
1
cat
name
name
1
Reject
character
character
1
\n
character
character
1
\n
character
character
5
data1
data2
alpha
plot

wrpcauchy.ml
.named_I
structure
2
function
function
5
missing
missing
0
missing
missing
0
missing
missing
0
numeric
numeric
1
1.0000000000000001e-015
{
{
9
<-
<-
2
name
name
1
mu1.old
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
rho0
call
call
2
name
name
1
cos
name
name
1
mu0
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
^
name
name
1
rho0
numeric
numeric
1
2
<-
<-
2
name
name
1
mu2.old
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
*
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
rho0
call
call
2
name
name
1
sin
name
name
1
mu0
(
(
2
name
name
1
(
call
call
3
name
name
1
+
numeric
numeric
1
1
call
call
3
name
name
1
^
name
name
1
rho0
numeric
numeric
1
2
<-
<-
2
name
name
1
w.old
call
call
3
name
name
1
/
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
*
name
name
1
mu1.old
call
call
2
name
name
1
cos
name
name
1
data
call
call
3
name
name
1
*
name
name
1
mu2.old
call
call
2
name
name
1
sin
name
name
1
data
<-
<-
2
name
name
1
flag
numeric
numeric
1
0
while
while
2
call
call
3
name
name
1
==
name
name
1
flag
numeric
numeric
1
0
{
{
5
<-
<-
2
name
name
1
mu1.new
call
call
3
name
name
1
/
call
call
2
name
name
1
sum
call
call
3
name
name
1
*
name
name
1
w.old
call
call
2
name
name
1
cos
name
name
1
data
call
call
2
name
name
1
sum
name
name
1
w.old
<-
<-
2
name
name
1
mu2.new
call
call
3
name
name
1
/
call
call
2
name
name
1
sum
call
call
3
name
name
1
*
name
name
1
w.old
call
call
2
name
name
1
sin
name
name
1
data
call
call
2
name
name
1
sum
name
name
1
w.old
<-
<-
2
name
name
1
diff1
call
call
2
name
name
1
abs
call
call
3
name
name
1
-
name
name
1
mu1.new
name
name
1
mu1.old
<-
<-
2
name
name
1
diff2
call
call
2
name
name
1
abs
call
call
3
name
name
1
-
name
name
1
mu2.new
name
name
1
mu2.old
if
if
3
call
call
3
name
name
1
&&
(
(
2
name
name
1
(
call
call
3
name
name
1
<
name
name
1
diff1
name
name
1
acc
(
(
2
name
name
1
(
call
call
3
name
name
1
<
name
name
1
diff2
name
name
1
acc
<-
<-
2
name
name
1
flag
numeric
numeric
1
1
{
{
3
<-
<-
2
name
name
1
mu1.old
name
name
1
mu1.new
<-
<-
2
name
name
1
mu2.old
name
name
1
mu2.new
<-
<-
2
name
name
1
w.old
call
call
3
name
name
1
/
numeric
numeric
1
1
(
(
2
name
name
1
(
call
call
3
name
name
1
-
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
*
name
name
1
mu1.old
call
call
2
name
name
1
cos
name
name
1
data
call
call
3
name
name
1
*
name
name
1
mu2.old
call
call
2
name
name
1
sin
name
name
1
data
<-
<-
2
name
name
1
mu.const
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
+
call
call
3
name
name
1
^
name
name
1
mu1.new
numeric
numeric
1
2
call
call
3
name
name
1
^
name
name
1
mu2.new
numeric
numeric
1
2
<-
<-
2
name
name
1
rho
call
call
3
name
name
1
/
(
(
2
name
name
1
(
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
2
name
name
1
sqrt
call
call
3
name
name
1
-
numeric
numeric
1
1
call
call
3
name
name
1
^
name
name
1
mu.const
numeric
numeric
1
2
name
name
1
mu.const
<-
<-
2
name
name
1
mu
call
call
3
name
name
1
%%
call
call
3
name
name
1
atan
name
name
1
mu2.new
name
name
1
mu1.new
(
(
2
name
name
1
(
call
call
3
name
name
1
*
numeric
numeric
1
2
name
name
1
pi
call
call
3
name
name
1
data.frame
name
name
1
mu
name
name
1
rho
character
character
5
data
mu0
rho0
acc