Pessoal;
Novamente recorro a ajuda de vcs...Estou tendo dificuldades no ajuste de um
modelo não linear aos meus dados.
Vejam mensagem de erro que está aparecendo.
> a=6> b=0.16> c=-0.125
> plot(HT~DAP,data=dados)> curve(1.30+(x^2/a+b*x+c*x^2),add=T)
> n <- nls(HT~(DAP^2/a+b*DAP+c*DAP^2)+1.3,data=dados,+
> start=list(a=a, b=b, c=c))Erro em nlsModel(formula, mf, start, wts) :
> matriz gradiente singular nas estimativas iniciais dos parâmetros
Em anexo segue a base de dados!
Muito obrigado pela ajuda!
Thiago de Paula Protásio
Acadêmico de Engenharia Florestal
Universidade Federal de Lavras
Laboratório de Energia da Biomassa Florestal
Laboratório de Biomateriais
(035) 9183-2246
Arvore DAP HT
1 8,30 9,00
2 13,60 8,65
3 13,40 10,75
4 20,40 11,87
5 10,10 9,14
6 14,40 9,95
7 13,10 9,70
8 15,60 10,50
9 13,40 9,00
10 13,60 8,55
11 14,40 9,90
12 4,40 5,15
13 10,40 9,25
14 10,40 8,78
15 11,60 8,30
16 12,10 8,80
17 12,40 9,65
18 13,40 10,93
19 10,60 9,25
20 14,60 9,48
21 12,20 10,25
22 14,10 10,00
23 10,40 9,26
24 7,20 7,47
25 4,40 5,83
26 10,10 8,40
27 14,60 10,45
28 8,10 7,67
29 15,10 11,10
30 14,10 10,95
31 11,60 11,10
32 16,30 11,70
33 11,40 10,36
34 13,60 11,15
35 11,40 8,80
36 13,10 10,05
37 13,40 9,60
38 13,10 11,50
39 11,40 10,40
40 10,10 10,65
41 12,40 10,77
42 17,40 10,62
43 10,30 11,68
44 10,40 8,45
45 7,40 5,64
46 14,60 11,37
47 15,30 11,10
48 13,60 12,47
49 10,40 9,65
50 12,10 9,85
51 10,60 10,30
52 12,20 9,86
53 9,40 5,90
54 10,60 9,92
55 15,30 9,36
56 7,10 8,24
57 16,80 11,68
58 16,10 12,13
59 19,60 12,21
60 8,40 8,20
61 12,40 10,00
62 15,40 11,16
63 12,60 10,35
64 12,60 8,10
65 12,60 10,30
66 9,40 8,87
67 10,50 8,75
68 10,10 8,80
69 13,40 9,45
70 8,10 8,75
71 11,60 10,45
72 14,40 11,82
73 11,80 10,55
74 13,60 10,15
75 15,60 11,10
76 10,30 10,23
77 8,60 8,25
78 9,40 10,30
79 9,10 9,70
80 12,40 9,72
81 11,40 10,27
82 5,20 5,90
83 11,60 9,60
84 12,50 10,10
85 7,50 8,12
86 8,40 8,55
87 9,30 9,17
88 5,10 5,47
89 18,40 12,68
90 16,10 10,50
91 6,00 6,88
92 9,40 9,72
93 12,50 11,50
94 9,50 10,18
95 11,20 10,50
96 11,50 10,26
97 6,50 7,54
98 11,60 8,70
99 9,50 9,70
100 11,50 8,43
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
