- R-br - DINF UFPR

Re: [R-br] correlação de Pearson
by Ivan Bezerra Allaman 02 Out '12

02 Out '12

cor.prob <- function(X, dfr=nrow(X)-2){ R <- cor(X,use='na.or.complete') above <- row(R) < col(R) r2 <- R[above]^2 Fstat <- r2*dfr/(1-r2) R[above] <- 1-pf(Fstat, 1, dfr) R } \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@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

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

Chave de licença pública para R: Ubuntu 12.04
by Vinícius Lionel Mateus 30 Set '12

30 Set '12

Prezados, Não consigo fazer atualizações através dos repositórios por conta da ausência da chave de licença pública. Segui as dicas que constam em http://cran.fiocruz.br/, mas na hora de atualizar o sistema "diz": "Necessita instalação de pacotes não confiáveis" "A ação deve requerer a instalação de pacotes de fontes não autenticadas." ...e por fim não consigo fazer o download. Alguém encontrou o mesmo problema no Ubuntu? Vinícius -- -- Atenciosamente, VINÍCIUS LIONEL MATEUS, M.Sc (http://lattes.cnpq.br/6501001637020665) Bacharel em Química - Doutorando em Química Analítica Laboratório de Química Atmosférica - Departamento de Química Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900 Telefone: (+55) (21) 3527-1327 (+55) (21) 9358-8051 www.puc-rio.br

3 5

Não consigo ler shapefile
by Emerson Cotta Bodevan 28 Set '12

28 Set '12

Prezados, tenho mais uma dúvida. Primeiramente, uso o UBUNTU e o EMACS. Estou tentando ler um shapefile que baixei do site do IBGE. São dois arquivos: 31mu2500gsd.shp e 31mu2500gsr.shp Instalei o pacote maptools e não consegui instalar o rgeos (dúvida enviada para a lista em outro e-mail) ao digitar o comando: mg<-readShapePoly("31mu2500gsr.shp") ou mg<-readShapePoly("31mu2500gsd.shp") Aparece a mensagem mg<-readShapePoly("31mu2500gsd.shp") Erro em make.names(onames, unique = TRUE) : string multibyte inválida 5 mg<-readShapePoly("31mu2500gsr.shp") Erro em make.names(onames, unique = TRUE) : string multibyte inválida 5 Qualquer ajuda será bem vinda. Att., Emerson.

4 5

Re: [R-br] (sem assunto)
by Ivan Bezerra Allaman 28 Set '12

28 Set '12

?descr \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@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

4 3

Modo mais inteligente?(STRING)
by Ivan Bezerra Allaman 28 Set '12

28 Set '12

Bom dia senhores!! Suponha que eu tenha o seguinte vetor: vet <- paste(seq(1:12),1:5,sep='.') Eu terei sempre um 'ponto' e um número a direita de cada elemento do vetor. O meu objetivo, é pegar exatamente qualquer coisa que esteja à esquerda deste ponto. Então fiz: unlist(lapply(strsplit(vet,split='[[:punct:]][[:alnum:]]'),function(x)sample(x[1]))) Tem algum modo mais inteligente? Desde já, grato! \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@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

4 4

Remover coluna de um data.frame dentro de uma lista!
by Ivan Bezerra Allaman 28 Set '12

28 Set '12

Bom dia senhores! Novamente venho recorrer a ajuda dos senhores! Vamos supor o seguinte: DADOSS <- expand.grid(sitio=1,ponto=as.character(c('SS')),serie=2,camara=as.character(LETTERS[1:5]),tempo=c(0,15,30,45),armazen=c('VIDRO','HEAD','PP')) which <- 'armazen' x <- DADOSS[,c(1,3,5)] y <- DADOSS vars <- list() for(i in levels(y[[which]])){ vars[[i]] <- subset(cbind(x,groups=y[[which]]), groups==i) vars[[i]][,!(colnames(dd[[i]]) %in% c('groups'))] } vars Vejam que o meu objetivo é remover a coluna 'groups' de cada lista. No entanto, a linha " vars[[i]][,!(colnames(dd[[i]]) %in% c('groups'))]" não funciona! Alguém tem alguma sugestão? OBS: 'Éder, peguei o seu CMR emprestado ok?' Desde já, grato novamente! \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@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 3