Agosto 2011 - R-br

Modelos não lineares....
by Thiago De paula protásio 09 Ago '11

09 Ago '11

Pessoal; Novamente recorro a ajuda de vcs...Estou tendo dificuldades no ajuste de um modelo não linear aos meus dados. Vejam mensagem de erro que está aparecendo. > a=6> b=0.16> c=-0.125 > plot(HT~DAP,data=dados)> curve(1.30+(x^2/a+b*x+c*x^2),add=T) > n <- nls(HT~(DAP^2/a+b*DAP+c*DAP^2)+1.3,data=dados,+ start=list(a=a, b=b, c=c))Erro em nlsModel(formula, mf, start, wts) : matriz gradiente singular nas estimativas iniciais dos parâmetros Em anexo segue a base de dados! Muito obrigado pela ajuda! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

3 2

mudança de nomes de colunas (variáveis) em uma lista
by Cassiano 09 Ago '11

09 Ago '11

Olá pessoal, sou iniciante no R e estou quebrando a cabeça em um problema que, acho, deve ter uma solução simples. Tenho uma lista de data.frames com os mesmos nomes nas colunas (nomes das variáveis). No total são 8 itens (data.frames). Quero mudar o nome de algumas variáveis e estou usando para isso, gsub(). Com um laço for consegui executar a tarefa. O código que usei foi for(i in 1:length(Indicador)) { names(Indicador[[i]]) <- gsub("Pontos", "Total.pontos.atend", names(Indicador[[i]])) names(Indicador[[i]]) <- gsub("Comercio", "Comercio", names(Indicador[[i]])) names(Indicador[[i]]) <- gsub("(.)*Folha", "Liberada.folha", names(Indicador[[i]])) names(Indicador[[i]]) <- gsub("(.)*Atendimento", "Ponto.atend", names(Indicador[[i]])) } onde, Indicador = list(), e as strings dentro do gsub são os nomes da variáveis originais e os nomes novos, respectivamente. Porém acredito que exista uma forma mais eficaz, provavelmente sem usar o laço for. Tentei de outra maneira (para o primeiro gsub do laço anterior) com > gsub("Pontos", "Total.pontos.atend", lapply(Indicador, names)) Esse último gsub fez as mudanças, porém gera uma saída como [1] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [2] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [3] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [4] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [5] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [6] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [7] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [8] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" [9] "c(\"Nome\", \"Codigo\", \"Lotericos\", \"SAA\", \"SNC\", \"PAE\", \"Total.pontos.atend\")" Aqui, "Pontos" foi alterado para "Total.pontos.atend". Se isso estiver certo como faço para atribuir o novo nome para os nomes das variáveis dos data.frames em cada item da lista "Indicador". Abraços! Cassiano Novais Consultor economista

2 1

Barras de erro em barplot
by Marcelo Claro de Souza 09 Ago '11

09 Ago '11

Possuo uma planilha distribuída da seguinte forma: cidade pH desvpad..... Brasília 4,5 0,5 Itirapina 4,3 0,2 Como faço p plotar acrescentar o desvio padrão superior e inferios a minha figura partindo de barplot(pH~cidade,data=dados, .....) Muito obrigado.

3 2

Re: [R-br] Leitura arquivo rtf
by Roney Fraga Souza 08 Ago '11

08 Ago '11

Eder, Tentei abrir esse arquivo (ST-26-E.rtf) no Mac OS pelo TextEdit, Word, LibreOffice, Pages e Xcode, alguns apresentaram que o arquivo pode estar corrompido e outros não deram mensagem de erro, mas também não abrem. Só visualizei o conteúdo do documento pelo Komodo Edit, ou quando modifiquei a extensão para .txt, será que esse arquivo não esta corrompido? Att Roney

3 2

Re: [R-br] Barras de erro em barplot
by Ivan Bezerra Allaman 08 Ago '11

08 Ago '11

Ou, PH <- c(4.5, 4.3) error <- c(0.5,0.2) cidade <- c('Brasília','Itirapina') da <- data.frame(cidade,PH,error) library(ggplot2) limits <- aes(ymax=PH + error, ymin=PH - error) p <- ggplot(da, aes(cidade,PH)) p + geom_bar(fill='blue') + geom_errorbar(limits) Ultimamente estou ficando muito satisfeito com esta library. Para maiores detalhes, http://had.co.nz/ggplot2/ Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

função diff com variavel data tempo
by Vinicius Brito Rocha 08 Ago '11

08 Ago '11

Estou trabalhando em um código onde basicamente preciso calcular as diferenças de ordem 1 para uma variavel com o formato "2011-08-08 12:30:00" "2011-08-09 13:30:00" . . . "2011-08-15 09:30:00" . . . O problema surgiu quando quando variava o parâmetor lag da função e armazenar o primeiro valor das diferenças. guarda=0 for ( i in 1:10){ guarda[i]=diff(vetor_data_tempo,lag=i)[1] } pois dependendo da distância/tempo entra linhas a notação da resposta muda de minuto, para hora e por aí vai.. perciso fixar esta notação, alguém saberia como fazer? Abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

1 0

Leitura arquivo rtf
by Eder David Borges da Silva 06 Ago '11

06 Ago '11

Pessoal, tenho arquivos neste formato rtf (primo distante do .doc e CIA).... Arquivo em: http://www.leg.ufpr.br/~eder/ST-26-E.rtf gostaria de saber se alguém já conseguiu carregar um documento com estas características, pois aparentemente é uma tabela, porem parece que cada célula é uma caixa de texto, muito estranho. Quanto utilizo o readLines vem uma serie de comandos de formatação: 1] "{\\rtf1\\ansi{\\fonttbl{\\f0\\fcharset0\\fprq0 Arial;}}{\\colortbl;\\red0\\green0\\blue0;\\red255\\green255\\blue255}" Att Éder

3 5

Re: [R-br] Teste de Normalidade Multivariado! Erro!
by Thiago De paula protásio 06 Ago '11

06 Ago '11

Ok Ivan... Obrigado pela dica....Eu encontrei o erro aqui......Estava usando a função matrix() e na verdade no meu caso a função é as.matrix(). C=as.matrix(carvão[,(9:14)]) mvShapiro.Test(C) #Teste de normalidade multivariado Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Re: [R-br] Teste de Normalidade Multivariado! Erro!
by Ivan Bezerra Allaman 06 Ago '11

06 Ago '11

A mensagem de erro já diz tudo!! X <- matrix(rnorm(22),ncol=2) mvShapiro.Test(X) X1<-matrix(rnorm(24),ncol=2) mvShapiro.Test(X1) Sem um exemplo reproduzível fica difícil verificar qual o real problema dos seus dados. Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

1 0

Teste de Normalidade Multivariado! Erro!
by Thiago De paula protásio 06 Ago '11

06 Ago '11

Prezados membros da R-Br; Estou tendo um problema com o teste de normalidade multivariado da função mvShapiro.Test() do pacote mvShapiroTest Estou trabalhando com 6 variáveis e 19 observações por variável....Talvez, vcs possam me ajudar a resolver esse erro. C=matrix(carvão[,(2:7)])> is.matrix(D)[1] TRUE > C=matrix(carvão[,(2:7)])> mvShapiro.Test(C)Erro em mvShapiro.Test(C) : Sample size must be between 12 and 5000. Muito obrigado! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0