Estimados,

estoy tratando de ajustar un modelo lineal simple de dos segmentos con
pendientes, ordenadas al origen y punto de quiebre desconocidos, es decir
que pueden variar libremente.
Intente usar la función "segmented" del paquete del mismo nombre, pero me da
siempre el siguiente error:

Error en model.matrix(seg.Z, data = eval(obj$call$data)) :
        el argumento "seg.Z" está ausente, sin valor por omisión

Abajo pego el script que estoy usando y el objeto con los datos "lcru"


> lcru

     logAC logLC
1    0.71  0.50
2    0.81  0.59
3    0.82  0.61
4    0.83  0.61
5    0.92  0.73
6    0.97  0.77
7    0.99  0.68
8    1.04  0.83
9    1.10  0.90
10   1.10  0.87
11   1.11  0.89
12   1.18  0.97
13   1.19  0.94
14   1.19  0.99
15   1.19  0.97
16   1.20  1.01
17   1.27  1.07
18   1.30  1.11
19   1.33  1.13
20   1.33  1.15
21   1.35  1.16
22   1.35  1.15
23   1.36  1.17
24   1.37  1.18
25   1.38  1.19
26   1.38  1.18
27   1.38  1.21
28   1.38  1.19
29   1.40  1.17
30   1.40  1.23
31   1.41  1.22
32   1.41  1.23
33   1.41  1.22
34   1.42  1.16
35   1.43  1.25
36   1.45  1.27
37   1.46  1.28
38   1.47  1.28
39   1.48  1.33
40   1.49  1.32
41   1.49  1.33
42   1.49  1.29
43   1.49  1.33
44   1.50  1.32
45   1.50  1.34
46   1.52  1.34
47   1.52  1.35
48   1.53  1.34
49   1.53  1.35
50   1.54  1.38
51   1.55  1.39
52   1.56  1.40
53   1.57  1.39
54   1.57  1.41
55   1.59  1.42
56   1.60  1.43
57   1.60  1.41
58   1.61  1.46
59   1.61  1.44
60   1.62  1.44
61   1.62  1.45
62   1.62  1.48
63   1.63  1.47
64   1.63  1.46
65   1.64  1.53
66   1.64  1.51
67   1.65  1.49
68   1.65  1.49
69   1.65  1.46
70   1.66  1.51
71   1.66  1.52
72   1.66  1.51
73   1.66  1.48
74   1.66  1.51
75   1.67  1.51
76   1.67  1.52
77   1.68  1.51
78   1.69  1.54
79   1.69  1.54
80   1.70  1.56
81   1.70  1.57
82   1.70  1.55
83   1.71  1.56
84   1.71  1.57
85   1.71  1.59
86   1.73  1.59
87   1.73  1.61
88   1.74  1.61
89   1.74  1.61
90   1.76  1.65
91   1.76  1.59
92   1.76  1.58
93   1.77  1.66
94   1.78  1.63
95   1.78  1.66
96   1.78  1.66
97   1.78  1.66
98   1.79  1.68
99   1.79  1.64
100  1.79  1.68
101  1.80  1.70
102  1.81  1.73
103  1.81  1.71
104  1.81  1.73
105  1.83  1.73
106  1.83  1.67
107  1.83  1.74
108  1.83  1.72
109  1.83  1.74
110  1.83  1.72
111  1.84  1.76
112  1.84  1.77
113  1.84  1.77
114  1.85  1.80
115  1.86  1.81
116  1.86  1.87
117  1.87  1.84
118  1.88  1.83
119  1.89  1.89
120  1.90  1.83
121  1.92  1.88
122  1.94  1.91
123  1.95  1.97


> attach(lcru)
> Lmodel<-lm(logLC~logAC)
> Lmodel

Call:
lm(formula = y ~ x)

Coefficients:
(Intercept)            x
    -0.4026       1.1607


> segmented(Lmodel, Z=logAC, psi=1.8)

Error en model.matrix(seg.Z, data = eval(obj$call$data)) :
        el argumento "seg.Z" está ausente, sin valor por omisión

En que me equivoco?


Desde ya gracias!

-- 
Lic. Nahuel E. Farias
Laboratorio de Invertebrados
Departamento de Biología
Facultad de Ciencias Exactas y Naturales,
Universidad Nacional de Mar Del Plata, CC 1245, CP 7600 Mar del Plata,
Argentina
Tel. + 54 - 0223 - 4752426 (Int. 466)

	[[alternative HTML version deleted]]

