Hola Nahuel,
Leyendo la ayuda, veo que debes especificar el modelo y luego seg.Z (tienes
Z) como una formula. En tu caso seria lo siguiente:
segmented(Lmodel, seg.Z=~logAC, psi=1.8)
# Call: segmented.lm(obj = Lmodel, seg.Z = ~logAC, psi = 1.8)
#
# Meaningful coefficients of the linear terms:
# (Intercept) logAC U1.logAC
# -0.3209 1.0989 0.6926
#
#Estimated Break-Point(s) psi1.logAC : 1.753
Espero sea de utilidad,
Jorge Ivan Velez
2010/8/27 Nahuel Farias <>
> 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]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
[[alternative HTML version deleted]]