Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

Link para esta página de comparações

roteiros:matriz [2012/05/25 18:38] – [Entrando na matriz...] mortararoteiros:matriz [2024/01/09 18:18] (atual) – edição externa 127.0.0.1
Linha 70: Linha 70:
  
 ==== Uma ajuda ==== ==== Uma ajuda ====
-Abaixo uma função para projetar populações a partir da matriz de transição e do estado inicial (tmax é o tempo máximo de projeção). É basicamente o que fizemos anteriormente, mas agora com a função. +Abaixo uma tem uma função para projetar populações a partir da matriz de transição e do estado inicial (tmax é o tempo máximo de projeção). É basicamente o que fizemos anteriormente, mas agora com a ajuda da função. 
  
 <box blue 90% | Função proj.mat> <box blue 90% | Função proj.mat>
Linha 83: Linha 83:
  for(i in 2:(tmax+1))  for(i in 2:(tmax+1))
  {  {
- res.mat[i,]=res.mat[(i-1),%*% matproj+ res.mat[i,]=matproj %*% res.mat[(i-1),]
  }  }
 return(res.mat) return(res.mat)
Linha 100: Linha 100:
 </code> </code>
 </box> </box>
- 
 ===== Taxa de Crescimento  ===== ===== Taxa de Crescimento  =====
  
Linha 119: Linha 118:
   * projete a população a tempos mais longos!   * projete a população a tempos mais longos!
   * veja como se comporta a taxa de crescimento da população $\lambda= \frac{N_t}{N_{t-1}}$   * veja como se comporta a taxa de crescimento da população $\lambda= \frac{N_t}{N_{t-1}}$
-  * faça o mesmo variando a alguma parâmetro da matriz de transição+  * faça o mesmo variando algum parâmetro da matriz de transição
   * como se comporta essa taxa ao longo do tempo? Ele muda (qualitativamente) quando muda algum parâmetro da população (transições, estado inicial)?   * como se comporta essa taxa ao longo do tempo? Ele muda (qualitativamente) quando muda algum parâmetro da população (transições, estado inicial)?
  
Linha 231: Linha 230:
 return(mat) return(mat)
 } }
-#estado inical  +cory 
-n0=c(10,5,2)+# estado inical  
 +n0=matrix(c(10,5,2), ncol=1)
 ## tempo 1 ## tempo 1
-n1 = n0 %*% cory +n1 =  cory %*% n0  
-n1d= n0 %*% ddf(n0,cory,h=10, st=3)+n1 
 +n1d= ddf(n0,cory,h=10, st=3)%*%  n0 
 n1 n1
 n1d n1d
 ## tempo 2 ## tempo 2
-n2 = n1 %*% cory +n2 =  cory %*% n1  
-n2d=n1d %*% ddf(n1d,cory,h=10, st=3)+n2 
 +n2d= ddf(n1,cory,h=10, st=3)%*%  n1 
 n2 n2
 n2d n2d
 ## tempo 3 ## tempo 3
-n3 = n2 %*% cory +n3 =  cory %*% n2 
-n3d=n2d %*% ddf(n2d,cory,h=10, st=3)+
 n3 n3
-n3d+n3d= ddf(n2,cory,h=10, st=3)%*%  n1  
 +n2 
 +n2d
 </code> </code>
  
Linha 262: Linha 265:
  for(i in 2:(tmax+1))  for(i in 2:(tmax+1))
  {  {
- res.mat[i,]=res.mat[(i-1),] %*% ddf(res.mat[(i-1),], matproj,h, st)+ res.mat[i,]=ddf(res.mat[(i-1),], matproj,h, st)  %*%  res.mat[(i-1),]
  }  }
 return(res.mat) return(res.mat)
Linha 272: Linha 275:
 prop.estdd<-res.corydd/apply(res.corydd,1,sum) prop.estdd<-res.corydd/apply(res.corydd,1,sum)
 matplot(0:20,prop.estdd, type="l", lty=2:4, col=2:4) matplot(0:20,prop.estdd, type="l", lty=2:4, col=2:4)
 +## tmax = 100 
 +res.corydd<-proj.dd(n0,matproj=cory, h=100, st=3, tmax=100) 
 +res.corydd 
 +matplot(0:100,res.corydd, type="l", ylab="N"
 +prop.estdd<-res.corydd/apply(res.corydd,1,sum) 
 +matplot(0:100,prop.estdd, type="l", lty=2:4, col=2:4, ylab="Proporção das classes")
 </code> </code>
  
Linha 280: Linha 288:
   *3. o estado inicial da população influencia a projeção quanto ao estado final da população? Há alguma similaridade com relação ao modelo sem densidade-dependência?   *3. o estado inicial da população influencia a projeção quanto ao estado final da população? Há alguma similaridade com relação ao modelo sem densidade-dependência?
 </box> </box>
- 
  
 ===== Autovalores e autovetores ===== ===== Autovalores e autovetores =====
Linha 450: Linha 457:
  
 </code> </code>
 +
 +===== Extração e Manejo =====
 +
 +[[exercicios:manejo|]]
 +
roteiros/matriz.1337971090.txt.gz · Última modificação: 2024/01/09 18:17 (edição externa)
www.chimeric.de Creative Commons License Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0