Outubro 2012 - R-br

Dúvida pacote URCA
by Marcel Quintana 02 Out '12

02 Out '12

Caros, alguém saberia me informar se na função ur. za, para detecção de raiz unitária, as possível quebra estrutural acusada pelo teste é do tipo I, II, ou III (additive outlier model, innovational outlier model ou ambos)? Grato.

1 0

Somar mais de uma variável
by Carlos Mendonça 01 Out '12

01 Out '12

Caros, existe uma forma mais prática de somar várias variáveis por um determinado by, ao invés da maneira que faço atualmente? Arquivo de dados (arquivo) codigo valor1 valor2 valor3 a 3 25 30 a 5 10 50 a 2 15 20 b 15 85 20 b 5 15 10 c 2 5 25 c 4 7 5 c 2 5 8 c 2 3 12 1) Utilizo o comando arq1 = aggregate(arq1$valor1, list(arquivo$codigo), sum) para somar o valor1 por cada codigo. 2) Utilizo o comando arq2 = aggregate(arq1$valor2, list(arquivo$codigo), sum) para somar o valor2 por cada codigo. 3) Utilizo o comando arq3 = aggregate(arq1$valor3, list(arquivo$codigo), sum) para somar o valor3 por cada codigo. 4) Depois junto o arq1, arq2 e arq3 e fico com um arquivo igual ao abaixo: codigo valor1 valor2 valor3 a 10 50 100 b 20 100 30 c 10 20 50

2 2

Re: [R-br] Função Gnls
by Walmes Zeviani 01 Out '12

01 Out '12

Rafael, Talvez você esteja passando chutes iníciais ruins. Faça gráficos de dispersão da resposta conta as covariáveis para obter chutes melhores ou use outro meio de pré-estimação. Envia sua resposta para lista de discussão do R, R-br, pois mais pessoas podem aproveitar-se da mesma dúvida, mais pessoas podem colaborar com você, você pode ser respondido em menos tempo, você pode contemplar diferentes soluções para o seu problema e sua mensagem fica no histórico da lista para consulta futura. Detalhes de inscrição na lista disponíveis em http://www.leg.ufpr.br/doku.php/software:rbr Não esqueça de ler o guia de postagem e submeter um CMR na sua mensagem disponíveis em http://www.leg.ufpr.br/doku.php/software:rbr#guia_de_postagem_da_r-br_em_co… https://gist.github.com/1088208 À 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(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2012/10/1 <rafaelbrco(a)hotmail.com> > Olá, então to com o seguinte problema no R: > > Tenho este modelo abaixo: > > Gompertz=nls(Ht~SSgompertz(Dap,b0,b1,b2)) > summary(Gompertz) > > Com os seguintes resultados para os seus parâmetros > > Formula: Ht ~ SSgompertz(Dap, b0, b1, b2) > > Parameters: > Estimate Std. Error t value Pr(>|t|) > b0 16.35440 0.98953 16.53 <2e-16 *** > b1 1.67231 0.67153 2.49 0.0133 * > b2 0.81996 0.05814 14.10 <2e-16 *** > --- > Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > Residual standard error: 4.704 on 310 degrees of freedom > > Number of iterations to convergence: 7 > Achieved convergence tolerance: 9.538e-07 > > Porem quando incluo a covariavel Idade (I) aparece o seguinte erro: > > > Gompertz.covariavel=gnls(Ht~SSgompertz(Dap,b0,b1,b2),params=list(b0~I,b1~I,b2~I),start=c(b0=16,0,b1=2,0,b2=1,0)) > > Erro em gnls(Ht ~ SSgompertz(Dap, b0, b1, b2), params = list(b0 ~ I, : > Step halving factor reduced below minimum in NLS step > > Alguem poderia me dar algumas dicas > >

1 0

Colocar frequências num gráfico de barras do ggplot
by Eduardo Carvalho da Silva 01 Out '12

01 Out '12

Olá a todos! Estou precisando gerar um gráfico de barras de uma variável X, particionando cada barra pela distribuição de uma outra variável Z. Até aí, tudo bem. O problema é que preciso colocar as frequências relativas de cada categoria da variável Z nas barras das categorias da variável X. Só consegui colocar as frequências considerando todas as observações. O que preciso é colocar as frequências considerando o número de observações em cada barra. Alguém saberia como resolver isso? Desde já agradeço, __ * Eduardo Carvalho da Silva *

3 2