[R-br] erro de objecto
alanarocha em sapo.pt
alanarocha em sapo.pt
Sexta Julho 12 13:23:08 BRT 2013
boa tarde
porque me dá este erro?
all.equal(integrate(f=fdp_expweibull,par=c(1,1,1),0,Inf)$value, 1)
Error in match.fun(f) : object 'fdp_expweibull' not found
Ana
Quoting r-br-request em listas.c3sl.ufpr.br:
> Enviar submissões para a lista de discussão R-br para
> r-br em listas.c3sl.ufpr.br
>
> Para se cadastrar ou descadastrar via WWW, visite o endereço
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> ou, via email, envie uma mensagem com a palavra 'help' no assunto ou
> corpo da mensagem para
> r-br-request em listas.c3sl.ufpr.br
>
> Você poderá entrar em contato com a pessoa que gerencia a lista pelo
> endereço
> r-br-owner em listas.c3sl.ufpr.br
>
> Quando responder, por favor edite sua linha Assunto assim ela será
> mais específica que "Re: Contents of R-br digest..."
>
>
> Tópicos de Hoje:
>
> 1. Re: Simulação estocástica (Walmes Zeviani)
> 2. Re: Coordenadas e distâncias (Tito Conte)
> 3. Re: Soma de Quadrado Tipo III no R para analisar um fatorial
> desbalaceado (Ivan Bezerra Allaman)
> 4. Re: Soma de Quadrado Tipo III no R para analisar um fatorial
> desbalaceado (Cesar Rabak)
> 5. Gráfico de intervalo de confiança de modelo não linear
> (Fernando Antonio de souza)
> 6. Ordenação de boxplots pela mediana no ggplot2 (Nicolay Cunha)
> 7. Re: Gráfico de intervalo de confiança de modelo não linear
> (walmes .)
> 8. Re: Ordenação de boxplots pela mediana no ggplot2 (walmes .)
> 9. Re: Ordenação de boxplots pela mediana no ggplot2 (Nicolay Cunha)
> 10. Re: Gráfico de intervalo de confiança de modelo não linear
> (Fernando Antonio de souza)
> 11. Dados semanais (ARMIN KOENIG)
> 12. Re: Dados semanais (Edson Lira)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 11 Jul 2013 12:14:10 -0300
> From: Walmes Zeviani <walmes em ufpr.br>
> To: Hildemar Calixto <hildemarcalixto em yahoo.com>,
> r-br em listas.c3sl.ufpr.br
> Subject: Re: [R-br] Simulação estocástica
> Message-ID:
> <CAFU=Eka_ka3W7uF_i3rzRMmdkgNDK+SDDMv_mxGLEqTpGdTSng em mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hildemar,
>
> Envie sua mensagem para o endereço da lista (r-br em listas.c3sl.ufpr.br) e
> não para os administradores dela. Envite enviar arquivos anexos. A lista
> recomenda que você hospede em algum lugar (na public do DropBox por
> exemplo) e envie o link para acesso.
>
> À disposição.
> Walmes.
>
> ==========================================================================
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: walmes em ufpr.br
> skype: walmeszeviani
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
>
> 2013/7/11 Hildemar Calixto <hildemarcalixto em yahoo.com>
>
>> Prezados, gostaria de sua ajuda para resolver um problema de simulação, no
>> qual utilizo o método da transformação inversa para gerar valores
>> de uma v.a. com distribuição hipergeométrica. Conforme material anexo
>> (pág. 80) fiz uma função para tal fim, porém está fornece alguns valores
>> "NA" e não sei como solucionar.
>>
>> Agradeço desde já a ajuda.
>>
>> --------------------------------
>> Hildemar Calixto
>> Graduando em Estatística
>> Universidade Federal do Ceará
>> --------------------------------
>>
> -------------- Próxima Parte ----------
> Um anexo em HTML foi limpo...
> URL:
> <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/0e9d581e/attachment-0001.html>
>
> ------------------------------
>
> Message: 2
> Date: Thu, 11 Jul 2013 12:17:13 -0300
> From: Tito Conte <tito.conte em gmail.com>
> To: "r-br em listas.c3sl.ufpr.br" <r-br em listas.c3sl.ufpr.br>
> Subject: Re: [R-br] Coordenadas e distâncias
> Message-ID:
> <CACqq46xQ0Hm8PnPE0jHiChvqu8vhRtiRpocdX9G3HqKEAAH=9g em mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> utilize um laço,
> o comando for ou while resolvem, dependendo do que você deseje fazer
>
> Tito Conte
>
>
>
> 2013/7/8 Roberto Guevara <robertoguevara.bio em gmail.com>
>
>> Olá pessoal do R-br.
>>
>> Meus dados consistem basicamente em coordenadas e o horário em que
>> cada ponto foi obtido, como exemplificado a seguir:
>>
>> obs dia min x y
>> 1 1 400 425 475
>> 2 1 410 425 425
>> 3 1 420 475 325
>> 4 1 430 525 325
>> 5 1 435 525 275
>> 6 1 443 525 275
>>
>> obs se refere a cada unidade de amostra, dia a cada dia de coleta (um
>> total de 73), min é o horário em unidades de minutos e x e y cada
>> elemento do ponto cartesiano.
>>
>> O que preciso obter desses dados é a distribuição de distâncias
>> euclidianas entre um ponto escolhido aleatoriamente e todos os pontos
>> seguintes de um mesmo dia. Ou seja, preciso escolher aleatoriamente um
>> ponto por dia, calcular, para todos os dias, a euclidiana entre este
>> ponto e todos os pontos posteriores do mesmo dia e relacionar a
>> distância espacial com a diferença de tempo entre os pontos e o dia em
>> que cada distância está associada.
>>
>> Eu já elaborei um código pra isso. Caso tenham sugestões para melhorar
>> o código, ficarei agradecido. Porém, não consigo passar para uma
>> próxima etapa. Como poderia fazer para repetir esse procedimento o
>> número de vezes que fosse necessário? Por exemplo, existem milhares de
>> combinações de pontos escolhidos aleatoriamente ao longo dos 73 dias.
>> Eu gostaria de repetir o procedimento 1000 vezes para ter uma ideia da
>> variação de distribuição de distâncias que posso obter em decorrência
>> da escolha aleatória de um ponto de origem por dia.
>>
>> Código já elaborado:
>>
>> arquivo<-read.table("arquivo.txt",head=T)
>>
>> fixo<-numeric() #vetor com origens de cada dia
>> ult<-numeric() #vetor com últimos pontos de cada dia
>> xy<-list() #lista para armazenar a "parte dos pontos de cada dia"
>> delimitadas por fixo e ult
>> hor<-list() #lista para armazenar a "parte dos horarios de cada dia"
>> delimitadas por fixo e ult
>> distancias<-list()#distancias euclidianas
>> tempo<-list()#diferença de tempo entre os pontos
>> dia<-list()#dia de coleta
>> for (i in 1:max(unique(arquivo$dia))){ #repetindo até o total de dias
>> fixo[i]<-sample(arquivo$obs[arquivo$dia==i],1)#sorteio ponto de
>> origem
>> ult[i]<-max(arquivo$obs[arquivo$dia==i])#armazenando último ponto
>> de cada dia
>> xy[[i]]<-arquivo[fixo[i]:ult[i],4:5]#armazenando parte dos pontos
>> de cada dia
>> hor[[i]]<-arquivo[fixo[i]:ult[i],6]#armazenando parte dos horários
>> de cada dia
>> distancias[[i]]<-dist(xy[[i]])[1:nrow(xy[[i]])-1]#somente as
>> distâncias em relação ao primeiro ponto interessam
>> tempo[[i]]<-dist(horarios[[i]])[1:length(horarios[[i]])-1]#somente
>> diferenças ao primeiro ponto
>> dia[[i]]<-rep(i,length(arquivo$dia[fixo[i]:ult[i]])-1)$cada dia de
>> coleta
>> }
>>
>> distancias<-unlist(distancias)
>> tempo<-unlist(tempo)
>> dia<-unlist(dia)
>>
>> dtd<-data.frame(distancias,tempo,dia)
>>
>> Este data.frame seria como uma unidade de amostra das milhares de
>> distribuições que são possíveis de obter. Como otimizar para que
>> qualquer número desejado não seja dispendioso de obter?
>>
>> Obrigado pela atenção
>>
>> --
>> Roberto Guevara Ferreira Lima
>>
>> Biólogo - CRBio 6ª Região: 73146/06-D
>> Mestre em Zoologia
>> Laboratório de Ecologia e Zoologia de Vertebrados
>> Instituto de Ciências Biológicas
>> Universidade Federal do Pará
>> Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110
>> Currículo Lattes: http://lattes.cnpq.br/7363382159007600
>> _______________________________________________
>> 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/20130711/69054e0a/attachment-0001.html>
>
> ------------------------------
>
> Message: 3
> Date: Thu, 11 Jul 2013 08:23:33 -0700 (PDT)
> From: Ivan Bezerra Allaman <ivanalaman em yahoo.com.br>
> To: R Brasil <r-br em listas.c3sl.ufpr.br>
> Subject: Re: [R-br] Soma de Quadrado Tipo III no R para analisar um
> fatorial desbalaceado
> Message-ID:
> <1373556213.31227.YahooMailNeo em web162301.mail.bf1.yahoo.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Esta dúvida já foi sanada a "n" tópicos atrás. Para obter a mesma
> soma de quadrados do tipo III do SAS, você precisa usar a mesma
> reparametrização do "DEMO", ou seja, SOMA ZERO.
>
> (s,f,p)
> Allaman
>
>
> \begin{signature}
> <<>>=
> Prof. Dr. Ivan Bezerra Allaman
> Universidade Estadual de Santa Cruz
> Departamento de Ciências Exatas e Tecnológicas
> Ilhéus/BA - Brasil
> Fone: +55 73 3680-5596
> E-mail: ivanalaman em yahoo.com.br/ivanalaman em gmail.com
> @
> \end{signature}
> -------------- Próxima Parte ----------
> Um anexo em HTML foi limpo...
> URL:
> <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/31c1bdcd/attachment-0001.html>
>
> ------------------------------
>
> Message: 4
> Date: Thu, 11 Jul 2013 15:11:17 -0200
> From: Cesar Rabak <cesar.rabak em gmail.com>
> To: r-br em listas.c3sl.ufpr.br, Ivan Bezerra Allaman
> <ivanalaman em yahoo.com.br>
> Subject: Re: [R-br] Soma de Quadrado Tipo III no R para analisar um
> fatorial desbalaceado
> Message-ID:
> <CAKrF98kGO-9eeDX00x9Q59MrHfTqj+LXqGQwqCE-CVcq-+7Tww em mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Ivan,
>
> Estou entendendo que você acha que já vale a pena termos um FAQ desta lista?
>
> []s
> --
> Cesar Rabak
>
>
>
> 2013/7/11 Ivan Bezerra Allaman <ivanalaman em yahoo.com.br>
>
>> Esta dúvida já foi sanada a "n" tópicos atrás. Para obter a mesma soma de
>> quadrados do tipo III do SAS, você precisa usar a mesma reparametrização do
>> "DEMO", ou seja, SOMA ZERO.
>>
>> (s,f,p)
>> Allaman
>> *
>> *
>> \begin{signature}
>> <<>>=
>> Prof. Dr. Ivan Bezerra Allaman
>> Universidade Estadual de Santa Cruz
>> Departamento de Ciências Exatas e Tecnológicas
>> Ilhéus/BA - Brasil
>> Fone: +55 73 3680-5596
>> E-mail: ivanalaman em yahoo.com.br/ivanalaman em gmail.com
>> @
>> \end{signature}
>>
>> _______________________________________________
>> 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/20130711/03299d22/attachment-0001.html>
>
> ------------------------------
>
> Message: 5
> Date: Thu, 11 Jul 2013 17:17:52 -0300
> From: Fernando Antonio de souza <nandodesouza em gmail.com>
> To: r-br em listas.c3sl.ufpr.br
> Subject: [R-br] Gráfico de intervalo de confiança de modelo não linear
> Message-ID:
> <CAFrWFs=1ozgQT9Y7vfJaczLPwkX8LNsj+taEiFFuj+1unxR0Kg em mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Caros amigos
>
> Estou tentando criar um gráfico de um modelo não linear utilizando a função
> gnls (pacote nlme). Utilizei essa função por que meus dados não possuem
> variâncias homogêneas e a função gnls me permite lidar melhor com isso.
>
> Pois bem. Dei uma lida no post do Walmes em seu blogo (Ridículas) sobre
> como fazer (
> http://ridiculas.wordpress.com/2011/05/19/bandas-de-confianca-para-modelo-de-regressao-nao-linear/)
> no entanto não consigo obter a matriz gradiente e Hessiana neste
> modelo
> e por isso não consigo dar continuidade.
>
> Abaixo segue o CMR, bem como os comando fornecido pelo Walmes até onde
> consegui evoluir. Este é meu gargalo!
>
> Desde já agradeço a atenção
>
> GM<-structure(list(Animal = c(1L, 60L, 62L, 68L, 72L, 7L, 73L, 98L,
> 99L, 100L, 91L, 92L, 93L, 51L, 83L, 89L, 86L, 87L, 94L, 95L,
> 17L, 82L, 84L, 97L, 3L, 76L, 78L, 48L, 52L, 55L, 69L, 13L, 22L,
> 33L, 63L, 12L, 67L, 77L, 80L, 41L, 42L, 46L, 50L, 25L, 53L, 57L,
> 64L, 28L, 35L, 58L, 61L, 21L, 24L, 54L, 59L, 20L, 30L, 38L, 45L,
> 10L, 11L, 19L, 9L, 14L, 18L, 32L, 56L, 2L, 4L, 5L, 23L, 27L,
> 29L, 40L, 47L, 6L, 8L, 31L, 34L, 36L, 39L, 44L, 65L), Manejo = c(1L,
> 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L,
> 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L,
> 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L,
> 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L,
> 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 2L, 2L), Gest = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 90L,
> 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L, 90L,
> 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L, 110L,
> 110L, 110L, 110L, 110L, 130L, 130L, 130L, 130L, 130L, 130L, 130L,
> 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 130L, 140L, 140L,
> 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L,
> 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L, 140L,
> 140L, 140L, 140L, 140L), Fetos = c(0L, 0L, 0L, 0L, 0L, 0L, 0L,
> 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), NaGLAMA = c(0.515523054,
> 0.274319898, 0.260406616, 0.345393579, 0.941416811, 3.242298603,
> 0.543581754, NA, NA, 0.310913983, 0.103119122, 0.242179757, 1.511982572,
> 0.248561812, 1.377624119, 1.239927094, NA, 0.347712662, NA, 0.907352204,
> NA, 1.629736193, 0.707560235, 1.260423046, 0.6423459, 1.314054397,
> NA, 0.728530524, 0.398189115, 0.749234643, 0.544891837, NA, 5.434432736,
> 1.273166074, NA, 1.055026755, 0.937103858, 1.297767845, 0.691943258,
> 1.336532109, NA, 1.69233929, 2.762457777, 1.639361362, 1.220678864,
> 0.696944641, 0.555373873, 2.694735252, 1.593673698, 4.8664716,
> 1.14279477, 0.952757874, 1.538877833, 3.417410969, NA, 3.827188165,
> 4.90094219, NA, 1.254578538, 3.997852624, 1.562503764, 2.689082873,
> 8.459237956, 8.787103701, 2.535628503, 4.715973548, 2.187817632,
> 1.860875669, 2.656679208, 4.929889812, 2.548334487, 4.177756252,
> 4.293566946, 5.717494612, 1.55775636, 4.922141642, 4.650077264,
> 6.112568346, 5.235183032, 1.907048841, 2.814221997, 5.477504341,
> 5.465705692)), .Names = c("Animal", "Manejo", "Gest", "Fetos",
> "NaGLAMA"), row.names = c(NA, -83L), class = "data.frame")
>
> library(nlme)
> NAGLM1<-gnls(NaGLAMA~a*exp(c*Gest),start=list(a=0.04,c=0.03),weight=varIdent(form=~1|Gest),na.action=na.omit,data=GM[-c(59,57,61),],subset=c(Gest>0
> & Fetos==1))
>
>
> #Códigos fornecidos pelo Walmes Adaptado
>
> criando função que representa o modelo, retorna gradidente e hessiano
>
> expo.der <- deriv3(~a*exp(c*Gest),
> c("a", "c"),
> function(Gest, a,c) NULL)
> str(expo.der)
>
> #-----------------------------------------------------------------------------
> # diagnose gráfica e primeiro chute
>
> plot(NaGLAMA~Gest, data=GM, xlab="idade Gestacional (dias)",
> ylab="Conteudo de sodio(g)")
> a <- 0.04; c <-0.03;
> curve(expo.der(x, a, c), add=TRUE, col=2)
> start <- list(a=a, c=c)
>
> #-----------------------------------------------------------------------------
> # ajustando o modelo aos dados a partir dos valores iniciais via gráfico
>
> n0 <- gnls(NaGLAMA~expo.der(Gest,
> a,c),weight=varIdent(form=~1|Gest),na.action=na.omit,
> data=GM[-c(59,57,61),], start=start)
> summary(n0)
> confint(n0)
> #----------------------------------------------------------------------------------------------------
> #Não consigo acessar estas matrizes. Não consegui progredi.
> str(n0)
> str(n0$m$fitted())
> c(n0$m$fitted())
> attr(n0$m$fitted(), "gradient")
> attr(n0$m$fitted(), "hessian")
> -------------- Próxima Parte ----------
> Um anexo em HTML foi limpo...
> URL:
> <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/daa99be6/attachment-0001.html>
>
> ------------------------------
>
> Message: 6
> Date: Thu, 11 Jul 2013 18:20:10 -0400
> From: Nicolay Cunha <nicolaycunha em gmail.com>
> To: r-br em listas.c3sl.ufpr.br
> Subject: [R-br] Ordenação de boxplots pela mediana no ggplot2
> Message-ID:
> <CAAPK02Akk8Yzv6bprcRXRf7ZyMig0AjeXC2TgZ9odhDr=LjFTg em mail.gmail.com>
> Content-Type: text/plain; charset=ISO-8859-1
>
> # Prezados,
> # tenho uma dúvida em relação a ordenação dos boxplots em relação à
> mediana em um gráfico de três fatores.
> # quero que os boxplots fiquem ordenados da maior para a menor mediana
> (escada), dentro de cada nível do fator que confere as cores no
> gráfico
> # abaixo o CMR
>
>
> ################################################################################################################################
> #install.packages("ggplot2")
> library(ggplot2)
> cmr <- read.table ('http://pastebin.com/raw.php?i=61J9W4Gt', header = T)
> str(cmr)
>
> #ordenando fatores
> cmr$height <- factor(cmr$height, levels = c('h','m','l'))
> cmr$Pulse <- factor (cmr$Pulse, levels = c("1LF_Low_Short",
> "2PF_Low_Short", "2PF_Medium_Medium",
> "3GF_Medium_Short","3GF_Medium_Long"))
>
> #plot
> ggplot(cmr, aes( height , variable,
> colour = factor(pop, levels = pop[order(factor(Pulse))]),
> fill = factor(Pulse)))+
> xlab ('Anther levels')+
> ylab ('Length (mm)')+
> scale_colour_manual( values = rep("black", 24), legend = FALSE) +
> scale_fill_manual(name = 'Flood Pulse intensity', values =
> c("#E5E5E5", "#ACACAC", "#747474", "#3B3B3B", "#030303"))+
> geom_boxplot(mapping = NULL, data = NULL, stat = "boxplot", position =
> "dodge", outlier.colour = "black",
> outlier.shape = 8, outlier.size = 1.5, notch = FALSE, notchwidth =
> 0.5, show_guide = TRUE)+
> theme_bw()
> ################################################################################################################################
>
>
> # a dúvida é, como ordenar os boxplots dentro de cada nível de Flood
> em ordem descrescente da mediana (varias escadinhas)?
> # []
>
>
> --
> Nicolay Leme da Cunha
>
> Biólogo, Mestre, Doutorando em Ecologia e Conservação
> Universidade Federal de Mato Grosso do Sul, 79070-900
> Campo Grande, MS, Brasil
> E-mail: nicolaycunha em gmail.com
> lattes.cnpq.br/5916316648872099
>
>
> ------------------------------
>
> Message: 7
> Date: Thu, 11 Jul 2013 19:32:00 -0300
> From: "walmes ." <walmeszeviani em gmail.com>
> To: r-br em listas.c3sl.ufpr.br
> Subject: Re: [R-br] Gráfico de intervalo de confiança de modelo não
> linear
> Message-ID:
> <CAFU=EkZZ7TPMxfvYfgh8d2GohxdFeVBSZibMRqDs1zgnU1-SvA em mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Para você extrair n0 tem que ser de classe 'nls' e você tá usando a 'gnls'.
> Então as matrizes por fora assim
>
> x <- seq(0,140,20)
> args(expo.der)
> m <- expo.der(Gest=x, a=coef(n0)["a"], c=coef(n0)["c"])
>
> c(m)
> attr(m, "gradient")
> attr(m, "hessian")
>
> À disposição.
> Walmes.
>
> ==========================================================================
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: walmes em ufpr.br
> skype: walmeszeviani
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
> -------------- Próxima Parte ----------
> Um anexo em HTML foi limpo...
> URL:
> <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/acbdb490/attachment-0001.html>
>
> ------------------------------
>
> Message: 8
> Date: Thu, 11 Jul 2013 20:31:54 -0300
> From: "walmes ." <walmeszeviani em gmail.com>
> To: r-br em listas.c3sl.ufpr.br
> Subject: Re: [R-br] Ordenação de boxplots pela mediana no ggplot2
> Message-ID:
> <CAFU=Eka7H0-6v2=zivwm0B_3G-r-R9msY+P++WsPA+JQczhbzQ em mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Eu sei fazer com a lattice. Está nesse script
>
> http://www.leg.ufpr.br/~walmes/cursoR/cnj/script2.R
>
> Vá executando desde o começo por causa das manipuções com a base de dados
> que é feita em sequência e tals. A parte que interessa é em torno da linha
> 405.
>
> À disposição.
> Walmes.
>
> ==========================================================================
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: walmes em ufpr.br
> skype: walmeszeviani
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
> -------------- Próxima Parte ----------
> Um anexo em HTML foi limpo...
> URL:
> <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/fcfcbeba/attachment-0001.html>
>
> ------------------------------
>
> Message: 9
> Date: Thu, 11 Jul 2013 20:20:04 -0400
> From: Nicolay Cunha <nicolaycunha em gmail.com>
> To: r-br em listas.c3sl.ufpr.br
> Subject: Re: [R-br] Ordenação de boxplots pela mediana no ggplot2
> Message-ID:
> <CAAPK02Cd18FgjkUwaFPCz7dOomB13=GtTVBGfQKDNK_kYE_aTw em mail.gmail.com>
> Content-Type: text/plain; charset=ISO-8859-1
>
> Obrigado Walmes, se não houver como fazer isso no ggplot2 partirei pro
> lattice como você sugere.
> Nicolay
>
>
>
> 2013/7/11 walmes . <walmeszeviani em gmail.com>:
>> Eu sei fazer com a lattice. Está nesse script
>>
>> http://www.leg.ufpr.br/~walmes/cursoR/cnj/script2.R
>>
>> Vá executando desde o começo por causa das manipuções com a base de dados
>> que é feita em sequência e tals. A parte que interessa é em torno da linha
>> 405.
>>
>> À disposição.
>> Walmes.
>>
>> ==========================================================================
>> Walmes Marques Zeviani
>> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
>> Departamento de Estatística - Universidade Federal do Paraná
>> fone: (+55) 41 3361 3573
>> VoIP: (3361 3600) 1053 1173
>> e-mail: walmes em ufpr.br
>> skype: walmeszeviani
>> twitter: @walmeszeviani
>> homepage: http://www.leg.ufpr.br/~walmes
>> linux user number: 531218
>> ==========================================================================
>>
>>
>> _______________________________________________
>> 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.
>
>
>
> --
> Nicolay Leme da Cunha
>
> Biólogo, Mestre, Doutorando em Ecologia e Conservação
> Universidade Federal de Mato Grosso do Sul, 79070-900
> Campo Grande, MS, Brasil
> E-mail: nicolaycunha em gmail.com
> lattes.cnpq.br/5916316648872099
>
>
> ------------------------------
>
> Message: 10
> Date: Thu, 11 Jul 2013 21:56:57 -0300
> From: Fernando Antonio de souza <nandodesouza em gmail.com>
> To: r-br em listas.c3sl.ufpr.br
> Subject: Re: [R-br] Gráfico de intervalo de confiança de modelo não
> linear
> Message-ID:
> <CAFrWFsmgR4r0YkFBbHw-7jvL3NoavXXnohOBJT60peUR5m4u-A em mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Mais uma vez, muito obrigado Walmes! Acabei de rodar aqui e pelo jeito é
> muito mais simples que pelo nls.
>
> Muito obrigado
>
>
>
>
> Em 11 de julho de 2013 19:32, walmes . <walmeszeviani em gmail.com> escreveu:
>
>> Para você extrair n0 tem que ser de classe 'nls' e você tá usando a
>> 'gnls'. Então as matrizes por fora assim
>>
>> x <- seq(0,140,20)
>> args(expo.der)
>> m <- expo.der(Gest=x, a=coef(n0)["a"], c=coef(n0)["c"])
>>
>> c(m)
>> attr(m, "gradient")
>> attr(m, "hessian")
>>
>> À disposição.
>> Walmes.
>>
>> ==========================================================================
>> Walmes Marques Zeviani
>> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
>> Departamento de Estatística - Universidade Federal do Paraná
>> fone: (+55) 41 3361 3573
>> VoIP: (3361 3600) 1053 1173
>> e-mail: walmes em ufpr.br
>> skype: walmeszeviani
>> twitter: @walmeszeviani
>> homepage: http://www.leg.ufpr.br/~walmes
>> linux user number: 531218
>> ==========================================================================
>>
>>
>> _______________________________________________
>> 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/20130711/c9902af6/attachment-0001.html>
>
> ------------------------------
>
> Message: 11
> Date: Thu, 11 Jul 2013 18:16:53 -0700 (PDT)
> From: ARMIN KOENIG <arminkoenig19 em yahoo.com.br>
> To: "r-br em listas.c3sl.ufpr.br" <r-br em listas.c3sl.ufpr.br>
> Subject: [R-br] Dados semanais
> Message-ID:
> <1373591813.54894.YahooMailNeo em web162004.mail.bf1.yahoo.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Prezados,
> Estou lidando com uma série temporal, que tem 55 semanas.
> Para realizar os cálculos de previsão, utilizei:
> frequency=7, start=c(....), end=c(...) e funcionou.
> Mas, quando pedi plot (gerar o gráfico), o R devolveu-me um eixo X
> muito estranho.
> Minha dúvida é: como faço para que o "começo" (semana) e o "final"
> (também semana), seja plotado adequadamente?
> Obrigado pela ajuda,
> Armin
> -------------- Próxima Parte ----------
> Um anexo em HTML foi limpo...
> URL:
> <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130711/61ca3c78/attachment-0001.html>
>
> ------------------------------
>
> Message: 12
> Date: Fri, 12 Jul 2013 05:51:41 -0700 (PDT)
> From: Edson Lira <edinhoestat em yahoo.com.br>
> To: "r-br em listas.c3sl.ufpr.br" <r-br em listas.c3sl.ufpr.br>, ARMIN
> KOENIG <arminkoenig19 em yahoo.com.br>
> Subject: Re: [R-br] Dados semanais
> Message-ID:
> <1373633501.42969.YahooMailNeo em web161206.mail.bf1.yahoo.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Poste um CMR. Umas linhas de seu banco de dados, por exemplo no
> DatafileHost, DropBox e poste o link para nossos amigos possam
> ajudá-lo, inclusive o seu script também pode ser postado.
>
>
> [ ]'s.
>
>
> Edson Lira
> Estatístico
> Manaus-Amazonas
>
>
> ________________________________
> De: ARMIN KOENIG <arminkoenig19 em yahoo.com.br>
> Para: "r-br em listas.c3sl.ufpr.br" <r-br em listas.c3sl.ufpr.br>
> Enviadas: Quinta-feira, 11 de Julho de 2013 21:16
> Assunto: [R-br] Dados semanais
>
>
>
> Prezados,
> Estou lidando com uma série temporal, que tem 55 semanas.
> Para realizar os cálculos de previsão, utilizei:
> frequency=7, start=c(....), end=c(...) e funcionou.
> Mas, quando pedi plot (gerar o gráfico), o R devolveu-me um eixo X
> muito estranho.
> Minha dúvida é: como faço para que o "começo" (semana) e o "final"
> (também semana), seja plotado adequadamente?
> Obrigado pela ajuda,
> Armin
> _______________________________________________
> 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/20130712/41966ff1/attachment-0001.html>
>
> ------------------------------
>
> Subject: Legenda do Digest
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
> ------------------------------
>
> Fim da Digest R-br, volume 31, assunto 15
> *****************************************
Mais detalhes sobre a lista de discussão R-br