- R-br - DINF UFPR

Modelo em aov semelhante ao do ExpDes.pt
by Marcelo Laia 10 Jan '19

10 Jan '19

Colegas, Minha aluna está rodando no ExpDes o seguinte modelo: attach(clorofila) fat3.dbc(Isolado,Genótipo,Tempo,Bloco,Clorofila, quali = c(TRUE, TRUE, TRUE), fac.names = c("Isolado","Genótipo","Tempo"), sigF = "0.05", sigT = "0.05" Qual seria o modelo semelhante para ela rodar diretamente no aov? Por exemplo: aov(Clorofila ~ Genótipo*Tempo*Isolado ... etc Muito obrigado -- Marcelo

2 2

Particionar soma de quadrado
by Andre Oliveira 10 Jan '19

10 Jan '19

Boa noite, dados que houve a internação entre fatores, preciso particionar a soma de quadrado. É apenas um exercício que estou fazendo co dados que está na web. Evidentemente que posso construir os contraste e testar o que me interessa e também poderia usar a ExpDes para isto. Alguém poderia ajudar com um script para tal desdobramento? Grato Pequei este exemplo! da <- read.table("http://www.leg.ufpr.br/~walmes/cursoR/anovareg.txt", header=TRUE, sep="\t") str(da) attach(da) # Este exemplo testa regressão! fit=aov(indice ~ cultivar*ordered(dose)) summary(fit) summary.lm(fit) # Gostaria de ajuda para desdobrar e testar doses para este modelo abaixo! fit=aov(indice ~ cultivar*factor(dose)+Error(bloco:(cultivar+dose)) summary(fit) Grato

3 5

Distribuição gama bivariada no R
by Ana Beatriz Souza Souza 20 Dez '18

20 Dez '18

Boa tarde pessoal! Estou trabalhando com distribuições gama bivariada e outras distribuições relacionadas. Alguma indicação de pacotes no R que abordam esses modelos? Desde já agradeço Att., Ana

2 1

Mojema (1977) - ponto de corte
by Ze Henrique 20 Dez '18

20 Dez '18

Boa tarde Srs. gostaria de saber se alguém sabe onde posso conseguir como fazer a análise de Mojema (1977) no r. Teste este que consegue dar a significância de ponto de corte nos nós dos agrupamentos. Desde já agradece a pela ajuda. -- *José Henrique Soler Guilhen* Doutorando em Genética e Melhoramento Universidade Federal do Espírito Santo - CCAUFES

1 0

capturar a saída do r
by Cinara Santos 19 Dez '18

19 Dez '18

Olá! Gostaria de saber da possibilidade de capturar o [console] do R de forma dinâmica, como no software STATA - tudo o que sai na tela é copiado para um arquivo de log. Se houver interrupção, o log registra até o momento da interrupção. E se fosse possível habilitar isso no R, me ajudaria muito!!! Digo isso porque estou usando o comando [sink()] sink(file = log_file, append = TRUE, type = "output") Ocorre que o código lida com dados muito volumosos, dentro de um loop (com sub-loops também). Assim, na minha inocência tentei o seguinte - são 3 loops, um dentro do outro. Ao final do loop mais interno pus: - desliga o [sink()], - fecha o arquivo de saida (log_file - no intuito de salvar o que já tem) - abre o arquivo de saida (log_file) - "liga" o [sink] -------------------------- saida <- paste(....) #monta o nome do arquivo de saida log_file <- file(saida, open = "wt") sink(file = log_file, append = TRUE, type = "output") # começa a "captura" LOOP1- inicio pegar referencia no [arquivo_a] LOOP2- inicio montar nome [arquivo_b] LOOP3(varrer [arquivo_b] - particionando) - inicio tarefa1 tarefa2 : tarefa_n sink() close(log_file) log_file <- file(saida, open = "at") sink(file = log_file, append = TRUE, type = "output") LOOP3( varrer [arquivo_b] - particionando ) - fim LOOP2 - fim LOOP1- fim -------------------------- Mas parece que o arquivo continua fechando só ao final do código, e com isso, se acaba a luz no meio do caminho, ou ocorre algum outro evento.... fico na saudade... Tem remédio????? ------------------------ Cinara Santos whatsApp (32) 9.8416-0441

3 2

Identificar cpf's
by Felipe Barletta 18 Dez '18

18 Dez '18

Prezados, Boa tarde, Tenho uma conjunto de dados que são mensagens de textos e preciso identificar os cpf's digitados nestas mensagens. Quando os cpf's são digitados sem pontos ou traços, consigo extrair esta informação. mas preciso selecionar também os casos como por exemplo: 000.000.000-00 Uma amostra de 69 mensagens dos dados estão nest link https://www.dropbox.com/s/cphsc02rz242z57/amostra.csv?dl=0 ## CMR d<-unlist(regmatches(dados[1:69,], gregexpr("[0-9,]+", dados[1:69,]))) d<- gsub(",", "", d) d <- d[d != ""] d[str_length(d)==11]

2 1

Calculo Distancia entre Coordenadas [loop]
by Bruno Rhuan 17 Dez '18

17 Dez '18

Boa Tarde, Estou escrevendo um algorítimo que calcule a distância entre uma serie de coordenadas e retorne qual o ponto mais próximo. Para isso, o código utiliza dois datasets de coordenadas geográficas (lat/long) e o calculo da distancia é feito utilizando a biblioteca "geosphere". Porém, basicamente todos os resultados do algorítimo sempre apontam para a ultimo ponto do segundo dataset como sendo o mais próximo, assim creio ter algum erro no loop. Alguma sugestão? Abaixo o código ______ #biblioteca necessária #install.packages("geosphere") library("geosphere") #carregar arquivos x <- read.csv(file.choose(),sep=";",dec=",") # carregar lista de pontos nominais TIM y <- read.csv(file.choose(),sep=";",dec=",") # carregar base de sites #converter site id em characther x$SiteID <- as.character(x$SiteID) y$SITE.ID <- as.character(y$SITE.ID) #Mostrar tabelas View(y) View(x) dista1 <- c() #distancia do ponto nominal ao site # Calculo Distancia for (i in 1:216) { coord1 <- c(x$LONG[i],x$LAT[i]) for (j in 1:11195 ) { coord2 <- c(y$Longitude[j],y$Latitude[j]) dista1[j] <- as.numeric(distm(coord1,coord2, fun = distGeo)) #Calculo da menor distancia if (j>1) { k <- j-1 if (dista1[j]<dista1[k]) { menor_dist <- dista1[j] x$Distancia[i] <- menor_dist x$SiteID[i] <- y$SITE.ID[j] } else { menor_dist <- dista1[k] x$Distancia[i] <- menor_dist x$SiteID[i] <- y$SITE.ID[k] } } } } write.csv(x, "Resultado.csv", row.names = FALSE) _______________ Resultado retornado: ï*..END_ID* *LAT* *LONG* *Distancia (m)* *SiteID* BASDR_1 -1.296.434 -3.847.421 144.669.250.579.953 YBARWL BALFS_2 -1.289.767 -3.832.365 146.401.362.045.206 YBARWL CEFLA_3 -377.709 -3.855.802 222.607.630.292.628 YBARWL Distancia (m) = distancia ao ponto mais próximo SiteID = nome do ponto mais próximo. At.te Bruno

2 1

Fwd: Criação de um dataframe a partir das diferenças
by SILVIA PATRICIA BARRETO SANTANA 14 Dez '18

14 Dez '18

----- Mensagem encaminhada ----- De: "SILVIA PATRICIA BARRETO SANTANA" <silpat(a)ufrb.edu.br> Para: "r-br" <r-br(a)listas.c3sl.ufpr.br> Enviadas: Terça-feira, 11 de dezembro de 2018 18:05:24 Assunto: Criação de um dataframe a partir das diferenças Bom dia pessoal. Estou com dificuldade em obter os dados que necessito e qualque r ajuda se r á bem vinda, po r isso já começo ag r adecendo. Com este banco de dados (abaixo) gosta r ia de obte r a seguinte info r mação (mais abaixo): Ho r a r io Nivel Status 01/01/2011 00:00 29 0 01/01/2011 00:01 29 0 01/01/2011 00:02 30 1 01/01/2011 00:03 31 1 01/01/2011 00:04 29 0 01/01/2011 00:05 29 0 01/01/2011 00:06 29 0 01/01/2011 00:07 29 0 01/01/2011 00:08 30 1 01/01/2011 00:09 32 1 01/01/2011 00:10 34 1 01/01/2011 00:11 35 1 01/01/2011 00:12 35 1 01/01/2011 00:13 25 0 O banco abaixo consiste em guardar a informação do horário em que o status muda de 0 para 1 com as diferenças dos níveis. Ho r a r io Nivel Status 01/01/2011 00:02 01 1 01/01/2011 00:03 01 1 01/01/2011 00:08 01 1 01/01/2011 00:09 02 1 01/01/2011 00:10 02 1 01/01/2011 00:11 01 1 01/01/2011 00:12 00 1 Muito grata! Silvia

2 1

IC para previsao
by Edmar Caldas 10 Dez '18

10 Dez '18

Preciso de uma ajuda para criar o IC para ajuste do modelo de regressão linear simples. não dá nenhuma mensagem de erro, cria somente a previsao do modelo em novos dados. lixo$previsao <-predict(simples, newdata=lixo,interval="confidence",level=0.95)lixo Edmar

4 3

Re: [R-br] Digest R-br, volume 96, assunto 8
by Edmar Caldas 10 Dez '18

10 Dez '18

Não funcionou a função prediction. O engraçado é quem nem mensagem de erro mostra. lixo$previsao <-predict(simples, newdata=lixo,interval = "prediction")lixo$previsao <-predict(simples, newdata=lixo,interval = "prediction", level=0.95) segue meu codigo, attach(waste)names(waste)simples <-lm(waste$wastee ~ restrnts)simplesajuste<-fitted(simples)waste$ajuste <-ajusteerro <- (waste$ajuste - waste$wastee)waste$erro<-erroprova <- waste$ajuste-waste$errowaste$prova <-provahead(waste) y<- (0.14685 + 0.01014 * waste$restrnts)waste$y<-y plot(wastee ~ ajuste) head(waste) #Prevendo em novos dados. library(readxl)lixo <- read_excel("Dieretorio/lixo.xlsx")View(lixo) attach(lixo) lixo$previsao <-predict(simples, newdata=lixo,interval="prediction",level=0.95)lixo eu consigo fazer a previsão, mas preciso do intervalo de confiança para previsão. Edmar nsegue Em segunda-feira, 10 de dezembro de 2018 12:00:18 BRST, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Re: IC para previsao (Elias T. Krainski) 2. Problema para plotar IC para modelo glm Gamma usando stat_smooth() no ggplot2 (ASANTOS) ---------------------------------------------------------------------- Message: 1 Date: Sun, 9 Dec 2018 16:54:45 +0000 (UTC) From: "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] IC para previsao Message-ID: <815218723.1276327.1544374485235(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Coloca "prediction" On Fri, Dec 7, 2018 at 22:39, Edmar Caldas por (R-br)<r-br(a)listas.c3sl.ufpr.br> wrote: Preciso de uma ajuda para criar o IC para ajuste do modelo de regressão linear simples. não dá nenhuma mensagem de erro, cria somente a previsao do modelo em novos dados. lixo$previsao <-predict(simples, newdata=lixo,interval="confidence",level=0.95)lixo Edmar_______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181209/9931f2c9/atta…> -------------- Próxima Parte ---------- Um texto embutido e sem conjunto de caracteres especificado foi limpo... Nome: Untitled Url: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181209/9931f2c9/atta…> ------------------------------ Message: 2 Date: Mon, 10 Dec 2018 01:47:12 -0300 From: ASANTOS <alexandre.santos(a)cas.ifmt.edu.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Problema para plotar IC para modelo glm Gamma usando stat_smooth() no ggplot2 Message-ID: <8d326872-10b9-06f7-cc1b-564caf54f3e7(a)cas.ifmt.edu.br> Content-Type: text/plain; charset=utf-8; format=flowed Prezados Membros do R-br, Estou tentando plotar um modelo glm Gamma com stat_smooth() no ggplot2 sem sucesso. Minha ideia é representar o intervalo de confiança para um modelo onde realizei a junção de dois níveis para o fator Feature. Segundo exemplo fictício do meu problema: #Simulação de um banco de dados, onde a variável resposta Production tem distribuição Gamma set.seed(123) d<-NULL N<-50 d$Production <- rgamma(N,10) d$Feature <- ifelse(d$Production >7 & d$Production<10, c("green"),ifelse(d$Production>11, c("red"), c("blue"))) d$Temp<-rnorm(N,20,5) d<-as.data.frame(d) # # Ajusto um modelo de Gamma completo mG<- glm(Production~ Feature + Temp, family= Gamma, data = d) summary(mG) anova(mG,test="Chi") # Comparo os níveis do fator Feature library(multcomp) PW<-summary(glht(mG, linfct = mcp(Feature= "Tukey"))) PW cld(PW) # Faço de conta que os níveis green = blue e realizo um novo ajuste Feature2<-d$Feature levels(Feature2) levels(Feature2)[1]<-"blue&green" levels(Feature2)[2]<-"blue&green" levels(Feature2) d$Feature2<-Feature2 mG2<- glm(Production~ Feature2 + Temp, family= Gamma, data = d) # Crio os valores de predição pred.data = data.frame( Feature2<-d$Feature2, Temp<-d$Temp ) pred.data$Production = predict(mG2, newdata=pred.data, type="response") #Realizo o plot library("ggplot2") ggplot(d, aes(Temp, Production, colour = Feature2)) + geom_point() + geom_line(data=pred.data) + stat_smooth(method = "glm", formula = y ~ x, family = Gamma) # E tenho uma representação errada do intervalo de confiança do meu modelo mG2 porque não consigo fazer a função stat_smooth() considerar o ajuste mG2. Como posso resolver isso? Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ====================================================================== ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------ Fim da Digest R-br, volume 96, assunto 8 ****************************************

3 2