[R-br] Derivar equações com R

Quinta Dezembro 6 20:31:52 BRST 2012

Fernando, você pode usar a função D() do R como no exemplo abaixo:
D(expression(x^2),"x")
No seu exemplo o R não vai reconhecer a função log10(x) para derivar. Mas como já sabemos da propriedade matemática que log10(x)=ln(x)/ln(10) e que a derivada de ln(x)=1/x então a derivada da expressão ln(x)/ln(10) = 1/x*ln(10)
No R você pode fazer da seguinte maneira:
D(expression(a+b*(log(PCVZ)/log(10))),"PCVZ")
Espero ter contribuído para sanar a sua dúvida qualquer duvida estou a disposição.
Att.
Tiago. 

Date: Thu, 6 Dec 2012 19:41:35 -0200
From: nandodesouza em gmail.com
To: r-br em listas.c3sl.ufpr.br
Subject: [R-br] Derivar equações com R

Caros colegas

Gerei uma equação da forma = log10(PB)~a+b*log10(PCVZ) o qual relaciona o conteúdo de proteína corporal animal em função de seu peso corporal vazio. Segundo a metodologia que estou seguindo esta equação deve ser derivada para se estimar o conteúdo de PB por Kg de PCVZ. Há uma forma de fazer esta derivada? Alguém poderia me ensinar como fazer?

grato pela atenção.

modelo5<-lm(log10(PB)~log10(PCVZ),data=PB)
summary(modelo5)

dput(PB)
structure(list(Animal = c(25L, 38L, 42L, 50L, 53L, 18L, 22L, 
34L, 41L, 44L, 48L, 2L, 12L, 14L, 39L, 46L, 4L, 5L, 21L, 26L, 
49L, 3L, 32L, 37L, 43L, 47L, 8L, 15L, 30L, 40L, 1L, 17L, 23L, 
28L, 7L, 9L, 19L, 35L), Gest = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 
0L, 0L, 0L, 0L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 100L, 
100L, 100L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 
140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L), Manej = c(1L, 
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 
2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 
1L, 2L, 2L, 2L, 2L), rep = c(4L, 5L, 3L, 1L, 2L, 1L, 2L, 6L, 
4L, 3L, 5L, 3L, 2L, 5L, 1L, 4L, 2L, 4L, 1L, 3L, 5L, 1L, 5L, 2L, 
3L, 4L, 1L, 4L, 3L, 2L, 2L, 3L, 4L, 1L, 3L, 2L, 4L, 1L), PV = c(40.5, 
40.9, 35.5, 34.5, 36.7, 26.9, 35.5, 31.2, 31.5, 30, 35.5, 40.8, 
36.6, 36.6, 36.5, 37.9, 41.9, 35, 35.5, 32.5, 42, 51.8, 50, 46.5, 
38.7, 42.7, 35.4, 38.4, 36.5, 31.2, 49, 46.4, 48.4, 42.8, 40.5, 
42.3, 50.2, 33.5), PCVZ = c(34.7811, 33.6882, 22.8156, 26.4014, 
26.7112, 19.1745, 26.6094, 22.4639, 25.6466, 22.1737, 25.773, 
35.4571, 30.8157, 28.1626, 31.8387, 30.3712, 32.9897, 26.5665, 
27.9632, 23.4118, 38.1716, 42.3432, 38.2408, 36.4839, 31.3696, 
36.7408, 25.9655, 28.1035, 28.9001, 26.0481, 39.8412, 37.9058, 
40.8529, 43.4908, 30.3399, 28.4752, 34.3194, 24.6319), PB = c(9.18910573, 
9.11864608, 5.653976965, 6.870682232, 6.838822605, 5.119182736, 
6.932858931, 6.047701201, 5.593529478, 5.379896679, 6.611844309, 
9.230423109, 8.184234256, 7.20898144, 8.703687287, 8.645435859, 
7.017519902, 5.976945686, 7.366941986, 4.644961891, 9.085415538, 
11.06095042, 9.461960295, 8.85340357, 8.008399724, 7.058114299, 
5.333367984, 5.670736687, 7.363494857, 6.014652251, 7.077788459, 
7.890982758, 8.935434043, 9.484189031, 6.804250745, 7.034024007, 
7.628728687, 5.124003503)), .Names = c("Animal", "Gest", "Manej", 
"rep", "PV", "PCVZ", "PB"), class = "data.frame", row.names = c(NA, 
-38L))

_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
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. 		 	   		  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121207/e863516a/attachment.html>