- R-br - DINF UFPR

Desdobrar modelos de regresão em fatorial duplo (um fator quantitativo)
by Marcelo Cardoso mello 12 Jun '11

12 Jun '11

Estou com dificuldades em desdobrar a interação em uma regressão em fatorial. Segue um exemplo abaixo. Como obter um modelo de regressão de r em função de b para a1 e a2?????? a<-as.factor(c(rep(1,9),rep(2,9))) b<-c(rep(0,3),rep(5,3),rep(10,3), rep(0,3),rep(5,3),rep(10,3)) r<-c(1,1.2,1.3,2,2.3,2.25,4,3.9,3.8,5,4.8,4.7,3.5,3,3.4,1.7,1.5,1.2) b<-as.factor(b) m<-lm(r~a*b) anova(m) b<-as.numeric(b) regressão<-lm(r~b) regressão anova(regressão) Como desdobrar b dentro de a1 e a2????? Obrigado. Marcelo.

4 9

Problema com o Excell.
by Edson Lira 11 Jun '11

11 Jun '11

Boa noite a todos, estou lendo um banco de dados do excell que transformei em csv em um Mac OS 10.6.6. Estou com o problema abaixo: id Nome Amostra idade cor est_civ escol alcool dro_ili 1 1 D.P.S Amostra 2 e 3 26 Faiodermo Solteiro 1 G N\xe3o N\xe3o 2 2 C.S.S.F Amostra 4 41 Faiodermo Casado 1 G Sim N\xe3o 3 3 A.M.S.A Amostra 5 e 6 37 Faiodermo Solteiro 2 G N\xe3o N\xe3o 4 4 F.M.M Amostra 7e 8 34 Faiodermo Amasiado 2 G N\xe3o N\xe3o 5 5 S.B.F Amostra 9 e 10 45 Faiodermo Solteiro 3 G Sim N\xe3o 6 6 W.S.A Amostra 11 e 12 32 Faiodermo Solteiro 2 G N\xe3o N\xe3o 7 7 A.S.M Amostra 13 e 14 50 Faiodermo Solteiro 2 G N\xe3o sim 8 8 T.F.D Amostra 15 e 16 21 Leucodermo Solteiro 2 G N\xe3o sim 9 9 E.S.C Amostra 17 e 18 32 Leucodermo Solteiro 2 G N\xe3o N\xe3o 10 10 D.A.V Amostra 19 e 20 44 Faiodermo Solteiro 1 G Sim N\xe3o O não está aparecendo N\xe3o. O excell que foi gerado este banco de dados em em espanhol, tentei corrigir e não consegui, alguém já teve esse problema? O que tentar fazer? Edson Lira Estatístico Manaus-Amazonas

2 1

Como adicionar ponto(s) de sela(s) em superfície resposta?
by Ivan Bezerra Allaman 11 Jun '11

11 Jun '11

Boa noite senhores!! Após ler quase todos os tópicos do nabble internacional de como adicionar pontos com "wireframe" e manipular de tudo qualquer jeito as funções panel.3d.INFERNO e outras milhões de funções da lattice para tentar fazer uma coisa simples, venho recorrer novamente a ajuda dos senhores. Tenho o seguinte CMR. require(plotrix) jet.colors <- colorRampPalette( c("blue", "green") ) x <- seq(-1.95, 1.95, length=30) y <- seq(-1.95, 1.95, length=35) da <- expand.grid(x=x, y=y) da$z <- with(da, x*y^2) require(lattice) panel.3d.contour <- function(x, y, z, rot.mat, distance, nlevels = 20, zlim.scaled, ...) { add.line <- trellis.par.get("add.line") panel.3dwire(x, y, z, rot.mat, distance, zlim.scaled = zlim.scaled, ...) clines <- contourLines(x, y, matrix(z, nrow = length(x), byrow = TRUE), nlevels = nlevels) for (ll in clines) { m <- ltransform3dto3d(rbind(ll$x, ll$y, zlim.scaled[1]), rot.mat, distance) panel.lines(m[1,], m[2,], col = add.line$col, lty = add.line$lty, lwd = add.line$lwd) } } wireframe(z~x+y, da, drape=TRUE, scales=list(arrows=FALSE),col.regions=jet.colors(100),panel.3d.wireframe="panel.3d.contour") Gostaria de indicar no gráfico o(s) ponto(s) (coordenadas x,y e z) de sela, para tornar o gráfico autoexplicativo. Desde já agradeço mais uma vez a atenção de todos. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 2

Re: [R-br] Atualização do R no Ubuntu 11.04
by Junior Beleti 11 Jun '11

11 Jun '11

Bom dia.. Estou tentando instalar o R 2-13 em outra máquina, compilando por meio do script que foi descrito neste tópico, porém estou recebendo as seguintes mensagens: ============================================= ============================================= ARQUIVOS FONTE R - VERSAO 2-13 ============================================= ============================================= ============================================= ============================================= ARQUIVOS FONTES R - VERSAO 2-13 ============================================= ============================================= Rpatch.sh: 46: aptitude: not found ============================================= ============================================= ARQUIVOS FONTE R-VERSAO 2-13 ============================================= ============================================= Rpatch.sh: 69: svn: not found Rpatch.sh: 69: ./tools/rsync-recommended: not found Rpatch.sh: 69: ./tools/link-recommended: not found Rpatch.sh: 88: ../R/configure: not found make: *** Nenhum alvo indicado e nenhum arquivo make encontrado. Pare. make: *** Sem regra para processar o alvo `vignettes'. Pare. make: *** Sem regra para processar o alvo `dvi'. Pare. make: *** Sem regra para processar o alvo `pdf'. Pare. make: *** Sem regra para processar o alvo `info'. Pare. ============================================= ============================================= ARQUIVOS FONTE R - VERSAO 2-13 ============================================= ============================================= make: *** Sem regra para processar o alvo `install'. Pare. make: *** Sem regra para processar o alvo `install-dvi'. Pare. make: *** Sem regra para processar o alvo `install-pdf'. Pare. make: *** Sem regra para processar o alvo `install-info'. Pare. ============================================= ============================================= ARQUIVOS FONTE R - VERSAO 2-13 ============================================= ============================================= Gostaria de saber o que pode estar errado, pois não consigo instalá-lo. Versão Ubuntu: 11.04 Obrigado. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

4 3

Eliminar NA e merge com colunas de dif tamanho
by Alexandre Heinemann 11 Jun '11

11 Jun '11

Caros Colegas, Estou com um problema de programação no R. No seguinte arquivo: http://sharex.xpg.com.br/files/6386884446/teste.txt.html, a coluna 1 representa a data. Essa data não considera os dias 31 dos respectivos meses. Nessa coluna Data ha NAs que representam datas inxistentes (30/2 e 29/2 para os anos que esse não existe). Necessito eliminar as linhas que possuem NA (somente para a coluna Data). Após eliminar os NAs, necessito fazer um "merge" com um vetor criado abaixo Data_1. Esse vetor Data_1 representa as datas "reais", ou seja, o ano com 365 e ou 366 dias. Esse merge deve ser feito em função da coluna teste$Data, com a finalidade de criar as datas referentes aos dias 31, preenchendo com NA as restantes. A saída que eu espero é: Data prec_mm tmax_C tmin_C rad_W.m2 file_name Lat Long Cenario rad_MJ.m2 tmax_memo tmin_memo rad_memo 1977-08-30 3.40726e-06 36.9722 23.0109 279.2510 lat_0.5_lon_46.2_P.dat -0.5 -46.2 P 24.127286 35.89341 22.94294 23.02423 1977-08-31 NA NA NA NA NA NA NA NA NA NA NA NA 1977-09-01 3.87843e-04 36.8744 23.0071 263.9960 lat_0.5_lon_46.2_P.dat -0.5 -46.2 P 22.809254 36.51460 22.99073 23.90967 O arquivo no qual estou trabalhando tem 1517400 linhas. http://sharex.xpg.com.br/files/6386884446/teste.txt.html (download) teste <- read.table("teste.txt",sep=";", header=T, colClasses=c("Date", "numeric","numeric","numeric","numeric","character","numeric","numeric","character", "numeric","numeric","numeric","numeric")) str(teste) head(teste) teste$Data = format(as.Date(teste$Data,'%d-%m-%y')) dim(teste) current = as.Date("1975-01-01") final = as.Date("1977-10-10") Data_1 <- seq(from=current,to = final, 1) length(Data_1) Obrigado, Alexandre

2 1

Re: [R-br] Pacote doMC ou doSMP
by Junior Beleti 11 Jun '11

11 Jun '11

Não conheço tais pacotes, mas vou tentar utilizá-los. Obrigado -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

1 0

modificar lavels de uma variável
by Edson Lira 10 Jun '11

10 Jun '11

Boa tarde a todos, Tenho um banco de dados, com estrutura abaixo, gostaria de mudar os níveis da variável Lab, por exemplo, aonde tem "ACIDO URICO" ou "ALBUMINA" ou "AMILASE" gostaria de atribuir o nome "Bioquímica", aonde é "17 ALFA HIDROXI PROGESTERONA" atribuir "Citoquímica" e etc... Como fazer isto? Descrição Convenio Quantidade Lab 1 17 ALFA HIDROXI PROGESTERONA <NA> NA 17 ALFA HIDROXI PROGESTERONA 2 17 ALFA HIDROXI PROGESTERONA GOLDEN CROSS 1 17 ALFA HIDROXI PROGESTERONA 3 17 ALFA HIDROXI PROGESTERONA PARTICULAR 1 17 ALFA HIDROXI PROGESTERONA 4 17 ALFA HIDROXI PROGESTERONA SUS 1 17 ALFA HIDROXI PROGESTERONA 5 17 ALFA HIDROXI PROGESTERONA UNIMED 2 17 ALFA HIDROXI PROGESTERONA 6 <NA> <NA> 5 <NA> 7 ACIDO URICO <NA> NA ACIDO URICO 8 ACIDO URICO AMIL 1 ACIDO URICO 9 ACIDO URICO CASSIS 1 ACIDO URICO 10 ACIDO URICO GEAP 5 ACIDO URICO []'s Edson Lira Estatístico Manaus-Amazonas

4 3

Plot em R studio
by Marcelo Claro de Souza 10 Jun '11

10 Jun '11

Estou rodando umas PCAs utilizando o R studio no Ubuntu 11, utilizando o pacote vegan pelos seguintes comandos: > dados=read.csv2(file.choose(),h=T,row.names=1) > dados=t(dados) > dados=as.data.frame(dados) > dados.pca=rda(dados,scale=TRUE) > summary(dados.pca) > simbolos=c(rep(19,7),rep(2,7),rep(3,7)) > plot(dados.pca,type="none",xlab="PC1 = 36%",ylab="PC2 = 14%") > points(dados.pca,display="sites",pch=simbolos) > legend("topleft",c("D","Sd","Eg"),pch=c(19,2,3))A análise está rodando certinha, mas na hora de exportar a figura não aprece a legenda na figura. Nem exportando como pdf e nem como png. Tem algum comando errado? Muito obrigado.

2 1

Pacote doMC ou doSMP
by Junior Beleti 09 Jun '11

09 Jun '11

Boa noite a todos. Estou utilizando os pacotes doMC e doSMP para paralelizar loops "for". No entanto, estou com problema de consumo de memória. Ambos os pacotes estão causando falta de memória. Tenho 4GB de memória real, e cheguei a deixar 20GB para swap. Gostaria de saber se alguém já utilizou tais pacotes ou já enfrentou problema semelhante. Obrigado.. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

Transformar tabelas em data.frames
by Edson Lira 09 Jun '11

09 Jun '11

Bom dia a todos, Tenho a estrutura abaixo: Descrição Convenio Quantidade 1 PROCEDIMENTO : 17 ALFA HIDROXI PROGESTERONA 2 GOLDEN CROSS 03 3 PARTICULAR 10 4 SUS 100 5 UNIMED 200 Como transformar a tabela em dados, sendo cada linha uma observação. Edson Lira Estatístico Manaus-Amazonas

2 1