- R-br - DINF UFPR

Valores lógicos.
by Juliana Freitas de Mello e Silva 18 Nov '11

18 Nov '11

Bom dia! Existe um valor lógico que signifique "até"? Eu queria selecionar todos os valores de "dd_stn" quando "frange" varia de 0 a 24, depois de 27 a 36, 39 a 48, ... prev = c(0,24,27,36,39,48,51,72,75,96,99,120) subset(dd_stn,dd_stn$frange==prev[j] ?? dd_stn$frange==prev[j+1]) Obrigada!

5 4

Pacotes não compatíveis com R 2.14
by Carlos Andrade 18 Nov '11

18 Nov '11

Meus cumprimentos a todos da lista. Atualizei a minha versão do R para 2.14. Utilizo Debian 6.0. A atualização via synaptic foi um pouco posterior ao aviso no CRAN de nova versão. Alguns dias depois acusou a atualização. Tentei instalar dois pacotes para análise envoltória de dados, Como root,tentei instalar dea, e veio a mensagem In getDependencies(pkgs, dependencies, available, lib) : package ‘dea’ is not available (for R version 2.14.0) Tentei instalar nonparaeff e veio a mensagem > install.packages("nonparaeff") Installing package(s) into ‘/usr/local/lib/R/site-library’ (as ‘lib’ is unspecified) tentando a URL ' http://cran-r.c3sl.ufpr.br/src/contrib/nonparaeff_0.5-2.tar.gz' Content type 'application/x-gzip' length 12470 bytes (12 Kb) URL aberta ================================================== downloaded 12 Kb * installing *source* package ‘nonparaeff’ ... ** Creating default NAMESPACE file ** R ** preparing package for lazy loading Error in library(pkg, character.only = TRUE, logical.return = TRUE, lib.loc = lib.loc) : package ‘Design’ does not have a NAMESPACE and should be re-installed ERROR: lazy loading failed for package ‘nonparaeff’ * removing ‘/usr/local/lib/R/site-library/nonparaeff’ The downloaded packages are in ‘/tmp/RtmpjwfNMx/downloaded_packages’ Mensagens de aviso perdidas: In install.packages("nonparaeff") : installation of package ‘nonparaeff’ had non-zero exit status > Tentei instalar o Design e veio a mensagem > install.packages("Design") Installing package(s) into ‘/usr/local/lib/R/site-library’ (as ‘lib’ is unspecified) Mensagens de aviso perdidas: In getDependencies(pkgs, dependencies, available, lib) : package ‘Design’ is not available (for R version 2.14.0) > Acho que o problema é com a versão 2.14 do R. Os pacotes não acompanharam ainda a atualização do R? Ou é problema da versão 2.14 do R? Agradeço qualquer ajuda. -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

2 5

Re: [R-br] Regressão segmentada!
by Ivan Bezerra Allaman 18 Nov '11

18 Nov '11

Valeu Walmes e PJ! Realmente fornecendo os pontos de corte a convergência ficou mais fácil. Mauro, Tem que ler a discussão desde o início! A função requerida está no início da discussão. Abraços! (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-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Re: [R-br] Regressão segmentada!
by Ivan Bezerra Allaman 17 Nov '11

17 Nov '11

Descobri o erro na função. O correto seria: segment_lpl <- function(x,b0,b1,b2,X0,X1){ (b0+b1*x)*(x<=X0)+ (b0+b1*X0)*(X0<x & x<X1)+ (b0+b1*X0-b2*(X1-x))*(x>=X1) } Só não entendo o porque que não estou conseguindo convergência, uma vez que estou fornecendo exatamente os valores que o gráfico interativo está fornecendo. Se alguém descobrir, fico grato. (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-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

4 3

Dúvida chisq.test
by Rafael Simões 17 Nov '11

17 Nov '11

Boa tarde! Estou tentando fazer um trabalho e está realmente dificíl. Será que alguém poderia me ajudar? Dois vetores, cada um com três elementos: *Curitiba <-c(24673, 505, 27) Lapa <- c(682, 2, 0)* Quando obtenho o teste obtenho os seguintes resultados: *chisq.test(Curitiba, Lapa)* "> chisq.test(Curitiba, Lapa) Pearson's Chi-squared test data: Curitiba and Lapa X-squared = 6, df = 4, p-value = 0.1991 Warning message: In chisq.test(Curitiba, Lapa) : Chi-squared approximation may be incorrect" Vi nos argumentos da chisq.test que para se fazer o teste, deve se usar Matriz ou vetores. Acontece que quando coloco os dois dados em formato de matriz, o resultado dá bem diferente: *m <- matrix(c(Curitiba, Lapa),2,3, byrow=TRUE) chisq.test(m)* "> chisq.test(m) Pearson's Chi-squared test data: m X-squared = 10.9144, df = 2, p-value = 0.004266 Warning message: In chisq.test(m) : Chi-squared approximation may be incorrect" A idéia é observar se existe algum tipo de associação entre Tipo de Gravidez e Munícipio. Os dados representam Um, dois, três ou mais filhos respectivamente. Algum dos meus testes está certo? Caso não estejam, alguém pode me informar onde eu estou errando? Essa é a primeira vez que trabalho com testes envolvendo a distribuição X² -- att, Rafael Simões (41) 9836 1010 http://alienado.com.br http://twitter.com/cwbrafael

5 6

D de Dunnet
by Peterson Gardin 17 Nov '11

17 Nov '11

Pessoal, alguém sabe me dizer como encontrar o D de Dunnet, tabelado no R. Obrigado. João Peterson Pereira Gardin Pesquisador Epagri - Videira - SC Doutor em Fisiologia Vegetal Engenheiro Agrônomo Fones: (49) 3566-0054

2 1

[Dúvida] Pintar polygono de Mapa
by Pedro Rafael 17 Nov '11

17 Nov '11

Pessoal plotei um arquivo shapfile de um mapa do Estado da Paraíba. Esse mapa apresente 16 regiões e cada região eu quero pintar com uma cor diferente (verde, azul, amarelo), corres não sequenciais. Estou conseguindo pintar os poligonos de cada região em separado. Criei uma variável no banco de dados (DBF) em que classifico o município enter 1 a 16. O comando abaixo eu consigo pintar de vermelho aqueles municípios que tem a variável dados_regionais igual à 1 e de azul os que tem a variável igual à 2. plot(mapa, col=ifelse(dados_regionais == 1, "red", "white"), axes = "TRUE") plot(mapa, col=ifelse(dados_regionais == 2, "blue", "white"), axes = "TRUE") O problema é que a cada plot tenho um mapa diferente. Como destacar esses municípios em um mesmo mapa? -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

3 5

Regressão segmentada!
by Ivan Bezerra Allaman 17 Nov '11

17 Nov '11

Bom dia senhores!! Estou encontrando um problema no ajuste de uma regressão segmentada, e não estou encontrando onde está este problema. Já fiz vários ajustes deste tipo, mais com este modelo particularmente tem alguma coisa errada que não sei onde é. No estudo dos pontos de chutes, o erro é claro, mais não estou sabendo como contornar. Deveria haver três segmentos, mais pelo estudo com janelas interativas, parece haver 5 segmentos. Se alguém puder ajudar, agradeço desde já. Segue um CMR. x <- c(0,0,0,10,10,10,20,20,20,35,35,35) y <- c(7.77, 7.81, 7.81, 7.60, 7.67, 7.65, 7.60, NA, 7.66, 7.60, 7.42, 7.50) segment_lpl <- function(x,b0,b1,b2,K,X0,X1){ (b0+b1*x)*(x<=X0)+ (K)*(X0<x & x<X1)+ (b0+b1*X0-b2*(X1-x))*(x>=X1) } teste <- nls(y ~ segment_lpl(x,b0,b1,b2,K,X0,X1), start=list(b0=7.8,b1=-0.02,b2=-0.01,K=7.5,X0=10,X1=20)) library(gWidgetsRGtk2) limits <- list(b0=c(7,8), b1=c(-0.05,0), b2=c(-0.04,0), K=c(7,8), X0=c(8,10), X1=c(18,25)) start <- list() # lista com os valores para chute #----------------------------------------------------------------------------- #função que será atualizada a cada movimento do deslizador #parâmetros dentro de svalue() são controlados, nomes igual aos da lista plot.chute <- function(...){ #faz o gráfico de dispersão plot(y ~ x) #sobrepõe a curva com os valores dos deslizadores curve(segment_lpl(x, svalue(b0), svalue(b1), svalue(b2), svalue(K),svalue(X0), svalue(X1)), add=TRUE, col=2) #reescreve o start com os valores dos delizadores, para usar na nls() start <- list(b0=svalue(b0), b1=svalue(b1), b2=svalue(b2), K=svalue(K), X0=svalue(X0), X1=svalue(X1)) } #----------------------------------------------------------------------------- #criação da janela com deslizadores #na primeira chamada escolher uma das opções (sempre escolho a 1) #Select a GUI toolkit #essa função pode estar num arquivo fn.R e carregada com source("fn.R") w <- gwindow("Caixa com deslizadores para controlar parâmetros") tbl <- glayout(cont=w) for(i in 1:length(limits)){ tbl[i,1] <- paste("Controle", names(limits)[i]) tbl[i,2, expand=TRUE] <- (assign(names(limits)[i], gslider(from=limits[[i]][1], to=limits[[i]][2], by=diff(limits[[i]])/20, value=mean(limits[[i]]), container=tbl, handler=plot.chute))) } # #----------------------------------------------------------------------------- #agora com a caixa criada, basta chamar a função e mover os deslizadores plot.chute() (S,P,f) 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-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

R rodando no terminal??
by Marcelo Cardoso mello 16 Nov '11

16 Nov '11

Por este dias iniciei a utilização do linux e muitos de vocês me deram excelentes dicas. No terminal (do ubunto) usei: $ sudo apt-get install r-base Depois: $ R E o R abre e roda no terminal. No linix o R roda no terminal ou faltou algo? Obrigado mais uma vez por alguma luz.

2 1

ks.test ajuda em resultados
by julio cesar araujo 16 Nov '11

16 Nov '11

Kolmogorov - Smirnov test no R, Pessoal, quero verificar se um conjunto de dados tem distribuição uniforme (0,1). Daí através da função ks.test, uso, por exemplo, o comando ks.test(x, "punif"). Para um conjutno de dados que possuo esse comando apresenta valor abaixo de 0.5 o que significa que meus dados não possuem a distribuição uniforme. Só que essa mesma função pode ser utilizada para comparar se duas amostras possuem a mesmadistribuição, por exemplo ks.test(x, y), daí criei uma variável y, através de um runif (uma uniforme randomica) e testei com a X e resultado foi aceito (acima de 0.20, no meu caso). Como isso é possível???? seria x então uniforme ou não, como é possível quando utilizado "punif" dar que não oe quando comparado a uma variável uniforme criada aleatorimente por mim acusar ter a mesma distribuição?? Desde já agradeço e peço desculpas pela "trivialidade" da pergunta.

6 7