- R-br - DINF UFPR

Re: [R-br] tratamento de pontos discrepantes (pesos de bovinos)
by sznelwar＠uol.com.br 20 Jul '17

20 Jul '17

2 1

CFP - IV WORKSHOP DE CIÊNCIA DE DADOS EDUCACIONAIS (WCDE’2017)
by Jorge Cavalcanti 20 Jul '17

20 Jul '17

------------------------------------------------------------ ---------------------------------------------- *IV WORKSHOP DE CIÊNCIA DE DADOS EDUCACIONAIS (WCDE’2017)* (antigo WORKSHOP DE MINERAÇÃO DE DADOS EDUCACIONAIS) VI Congresso Brasileiro de Informática em Educação (CBIE’2017) Recife, 30 de Outubro a 02 de Novembro de 2017 https://sites.google.com/view/wcde/ ------------------------------------------------------------ ---------------------------------------------- *Sobre o evento* A análise de dados educacionais, de uma maneira geral, representa uma área de pesquisa emergente em Informática em Educação para o desenvolvimento de métodos que exploram dados oriundos de ambientes educacionais e também administrativos com a finalidade de entender melhor os estudantes e os cenários em que eles aprendem. Dentro deste contexto, nas últimas três edições do Congresso Brasileiro de Informática em Educação (CBIE) esta temática estava sendo parcialmente tratada no Workshop de Mineração de Dados Educacionais (WMDE), por meio de tópicos de interesse que se restringem a dados gerados em ambientes de ensino-aprendizagem e com a aplicação de tarefas de mineração de dados como predição (regressão, séries temporais e classificação), agrupamento ou associação de dados, resultando em descobertas de conhecimento intrínseco nos dados. Contudo, ainda há outras áreas de pesquisa nesta temática como é o caso de Learning Analytics (LA) que faz uso de técnicas de análises de dados,incluindo tarefas de mineração de dados, para confirmar hipóteses colocadas em atividades de ensino que envolvem a aprendizagem do aluno, fomentando assim recursos analíticos ao tutor para entendimento e aprimoramento do ensino-aprendizagem. Por fim, existe também a temática de Academic Analytics (AA) que tem na essência as mesmas abordagens de análise utilizadas em MDE e LA, porém com alteração na fonte de dados, nesse caso advindo de sistemas educacionais administrativos. Como proposta de englobar as três áreas de pesquisa, uma vez que o foco principal de todas estas é analisar dados educacionais para o aprimoramento do ensino e aprendizagem, neste ano se fez a alteração do nome do workshop para Ciência de Dados Educacionais, ampliando o escopo do evento para uma abordagem mais ampla da temática análise de dados e fomentando a Comunidade de Informática em Educação um fórum para atualização, discussão e aplicação de pesquisas que tem como premissa o uso de dados oriundos da educação. O WCDE é parte do Congresso Brasileiro de Informática da Educação (CBIE), um evento anual da SBC, que ocorrerá este ano em Recife entre os dias 30 de Outubro e 02 de Novembro. Informações adicionais sobre o CBIE estão em https://sabertecnologias.com.br/cbie/. *Tópicos de Interesse* A lista abaixo apresenta alguns tópicos iniciais de interesse deste workshop. Porém, ela não é estritamente limitante quanto aos assuntos que podem ser tratados em Ciência de Dados Educacionais. • Modelagem analítica de Bancos de Dados Educacionais e integração com Dados Abertos (Open Data); • Mineração de Dados Educacionais e Personalização da Aprendizagem; • Infra-estrutura para análise de dados educacionais; • Aplicação de análise de dados educacionais no ensino (Superior ou Educação Básica); • Integrando dados educacionais e Internet das Coisas (IoT); • Análise de dados em ambientes emergentes pedagógicos como jogos educacionais, MOOCs e Sistemas Tutores Inteligentes; • Análise de dados na aprendizagem social e colaborativa; • Melhores práticas para adaptar o estado da arte das técnicas analíticas de recuperação de informação, sistemas de recomendação, análise de redes sociais, mineração de opinião para o domínio educativo; • Abordagens, métodos, algoritmos e ferramentas de analítica; • Análise de dados em apoio a tarefas educacionais administrativas; • Uso de dados combinados com teorias para a compreensão da aprendizagem; • Métricas de aprendizagem e avaliação; • Aprendizagem adaptativa e personalizada; • Valores, éticas e leis sobre os dados educacionais; • Visualização de Dados Educacionais; • Processos e metodologias para análises de dados educacionais. *Datas importantes* • Data limite para submissão de artigos: *20/08/2017* • Notificação de artigos aceitos: *08/09/2017* • Entrega da versão final dos artigos (com copyright): *15/09/2017* • Data limite para inscrição de pelo menos um autor por trabalho: * 15/09/2017*

1 0

Maximizar plot window no RStudio
by diogojose21＠yahoo.com.br 13 Jul '17

13 Jul '17

Salve Maria professor!!! Apenas complementando os colegas: na função "par", existe uma série de parâmetros que, associados a essa função, permitem delimitar as margens do papel (parâmetro “oma”) e da área em que o gráfico está inserido (parâmetro “mar”). Irás definir as margens para baixo (“bottom”), à esquerda (“left”), acima (“top”) e à direita (“right”). Uma observação de quem já passou por isso: dependendo do número de painéis, pode não ser útil os gráficos da função “graphics”. Nesse caso, pode ter de apelar para packages alternativas, como “lattice”, “ggplot2” e “plotly”. Esperando ter ajudado... In corde Jesu semper!!! _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

1 0

GLM multivariado
by ana paula coelho madeira 12 Jul '17

12 Jul '17

Bom dia! Alguém tem alguma referência de estudos com GLM multivariado para me passar? Desde já agradeço Ana Paula

5 4

Opções do lattice (histograma).
by Diogo Jerônimo 12 Jul '17

12 Jul '17

Prezados, boa tarde!!! Estou empregando a library "lattice" para construir alguns histogramas, como no exemplo abaixo: ###dados res=c(rnorm(1000,20,2),rnorm(1000,15,1.5),rnorm(1000,25,2.5),rnorm(1000,25,2.5)) ###fatores fac=c(rep("F1",1000),rep("F2",1000),rep("F3",1000),rep("F4",1000)) fac=factor(fac,levels=c("F1","F2","F3","F4"),labels=c("F1","F2","F3","F4")) ###gráfico em lattice require(lattice) histogram(~res|fac,xlab="Resultado",ylab="Freqüência",type=c("count"),layout=c(1,4),endpoints=c(0,40),nint=30,scales=list(x=list(at=c(0,10,20,30,40))),panel = function(x, ...) { panel.histogram(x,col=rgb(62,88,172,max=255),...) panel.abline(v=20,col=rgb(13,141,70,max=255),lwd=2)} ) Alguém sabe como fazer para inserir a curva normal no gráfico do histograma, quando se tem no eixo y a frequência ("count") em vez da densidade ("density")? Obrigado, Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQICONRE: 8514 - SÉRIE ACel:(21)99680-8535Email: diogojose21@yahoo.com.brhttp://lattes.cnpq.br/8996149312896520

2 1

Ajuda com BoxPlot
by Yury Duarte 11 Jul '17

11 Jul '17

Bom dia colegas listeiros, como vão? Nos últimos dias me vi com necessidade de elaborar gráficos do tipo BoxPlot pra poder representar a variabilidade da produção agrícola de cordo com a época de plantio. Não tive problemas em fazer o gráfico e visualizar ele no output do R-studio. Minha dúvida é em relação ao formato com que devo salvar a figura. Gostaria de salvar os gráficos em formato de imagem e com um tamanho pré definido, pois o default de salvamento das imagens fez com que eles ocupem toda uma pagina (quando salvei em PDF). Ao fim do email, segue parte do código que usei para gerar o BoxPlot. Desde já, agradeço pela ajuda de todos! a = read.table('sdesiderio.txt', header = T, sep = '\t') media = matrix(NA, nrow = 36, ncol = 4) colnames(media) = c('Desv.Pad', 'Media', 'M.mais', 'M.menos') for(i in 1:36){ d = sd(a[i,1:(ncol(a)-1)]) media[i,1] = d m = rowMeans(a, na.rm = T, dims = 1) media[i,2] = m[i] media[i,3] = media[i,2]+media[i,1] media[i,4] = media[i,2]-media[i,1] ifelse(media[i,]<0, 0, media[i,]) } a.box = t(a) colnames(a.box) = c('J1','J2','J3','F1','F2','F3','M1','M2','M3','A1','A2','A3','M1','M2','M3', 'J1','J2','J3','J1','J2','J3','A1','A2','A3','S1','S2','S3','O1','O2','O3', 'N1','N2','N3','D1','D2','D3') boxplot(a.box, las = 1, xlab = 'Semeadura', ylab = 'kg/ha', cex=0.5, pch=16) title('Variação de Ya de acordo com a data de semeadura') Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 4

Re: [R-br] incompatibilidade entre Geoda e R
by Edimeire Alexandra Pinto 10 Jul '17

10 Jul '17

Oi Elias. Pois é, olha só o que aparece: este.nb= read.gal("municipios.gal" , override.id=TRUE) print(este.nb) > print(este.nb) Neighbour list object: Number of regions: 5564 Number of nonzero links: 32228 Percentage nonzero weights: 0.104102 Average number of links: 5.792236 3 regions with no links: ######### w.nb<-nb2listw(este.nb,style="W", zero.policy=TRUE) | > print(w.nb, zero.policy=TRUE) Characteristics of weights list object: Neighbour list object: Number of regions: 5564 Number of nonzero links: 32228 Percentage nonzero weights: 0.104102 Average number of links: 5.792236 3 regions with no links: 3520400 5300108 2605459 Weights style: W Weights constants summary: n nn S0 S1 S2 W 5561 30924721 5561 2044.729 23174.73 | | | | | >w.nb ###se peço para aparecer o executado da erro como se vê abaixo Error in print.listw(x) : regions with no neighbours found, use zero.policy=TRUE Preciso desses pesos para calcular o indice de Moran da variavel Y, que nada mais do que o Renda municipal, usando matriz de pesos w.nb dados<-read.table("res.csv",sep=";",dec=",",header = TRUE) head(dados) | MUNICIPIO Y PIBa PIB cod 1 MANCIO LIMA 591239.9 33426.61 150136.08 1200336 2 MANOEL URBANO 58934.4 18044.54 80018.27 1200344 3 ASSIS BRASIL 1978559.5 17371.21 64601.58 1200054 4 FEIJO 285000.0 45966.75 267202.96 1200302 5 XAPURI 355711.4 55298.02 168708.32 1200708 6 TARAUACA 2852586.8 65297.63 335519.90 1200609 | | | | | > | | summary(dados) moran.test(dados$Y, listw=w.nb, randomisation=FALSE, alternative="greater",zero.policy=TRUE) Error in na.fail.default(x) : missing values in object | | Em Segunda-feira, 10 de Julho de 2017 15:22, Elias T. Krainski via R-br <r-br(a)listas.c3sl.ufpr.br> escreveu: Complementando: On 10/07/2017 15:20, Elias T. Krainski wrote: > nbw <- nb2listw(este.nb,style="W", zero.policy=TRUE) print(nbw, zero.policy=TRUE) Elias _______________________________________________ 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 fornea cdigo mnimo reproduzvel.

2 1

incompatibilidade entre Geoda e R
by Edimeire Alexandra Pinto 10 Jul '17

10 Jul '17

Boa tarde. Quem já tiver conhecimento sobre isso, ajudará muito no que puder colaborar. Alguém, por favor, que já mexeu com o pacote spdep, saberia me dizer o por que o R não faz o pesos da matriz de vizinhança, mesmo usando o argumento Zero.policy= TRUE. Tudo consiste em: 1- Crio o arquivo municipios.gal no Geoda.2- importo esse arquivo em formato .gal para o R: este.nb<- read.gal("municipios.gal" , override.id=TRUE) 3- tento criar a matriz de pesos de vizinhança : nb2listw(este.nb,style="W", zero.policy=TRUE) No entanto o R me retorna a mensagem abaixo, sendo que no Geoda eu consigo essa matriz de pesos. O fato é que não entendo por que o R não calcula essa matriz de pesos, já que todas as equações que terei de fazer serão no R, tendo o vista que o Geoda não tem essas equações, só o R Error in print.listw(x) : regions with no neighbours found, use zero.policy=TRUE

2 2

delineamento vs repetição
by Rafael Tieppo 10 Jul '17

10 Jul '17

Após consultar algumas referências (livros, apostila, papers, etc), não consegui a chegar a uma conclusão. O que encontrei de caso mais próximo foi uma solução do Walmes, que é o que tenho em mente até o momento. http://r-br.2285057.n4.nabble.com/R-br-Repeticao-dentro-da-parcela-td466187… Para acabar com as dúvidas, por favor, se alguém puder ajudar, segue a situação: No planejamento do experimento a ser realizado pretende-se verificar o espaçamento entre plantas em função do tipo de dosador de semente. - Var Dependente: espaçamento entre plantas (quantitativo) - Var Independente: tipo de dosador (qualitativo) Pretende-se usar DBC. A ideia é utilizar 4 blocos, com 4 tratamentos (dosadores). Sendo assim, cada bloco teria quatro parcelas (uma com cada tipo de sistema dosador), totalizando 4 repetições. Contudo, em cada parcela não pretende-se realizar apenas uma mensuração de distância entre plantas, deseja-se realizar pelo menos 10. O que gera minha dúvida são dois casos: CASO 1 O número de amostras é: 4 blocos * 4 tratamentos * 10 amostras (por parcela) = 160 CASO 2 Das 10 mensurações realizadas em uma parcela, calcula-se a média para a obtenção de uma única amostra por parcela, dessa forma tem-se: O número de amostras é 4 blocos * 4 tratamentos 1 amostra (por parcela)= 16 obs: O CASO 2 é o mais comum em alguns trabalhos, contudo não estou convicto que é a melhor opção. ### Perguntas: Há algum CASO correto? Se sim, em algum dos casos pode-se assumir como DBC ou DQL (quadrado latino)? Há outro delineamento que se enquadre melhor à situação? Obrigado. Rafael Tieppo

5 4

Re: [R-br] delineamento vs repetição
by Rafael Tieppo 10 Jul '17

10 Jul '17

Prezad@s (Gisele, João Henrique e Walmes), Muito obrigado pelas considerações. Gisele, Sobre o Arranjo fatorial, é interessante, vamos ver a possibilidade da execução. João Henrique, o conceito "a média das amostras simples é o resultado da parcela", ajudou a elucidar a dúvida maior sobre a repetição. Walmes, te respondendo, vc entendeu corretamente. Obrigado pela aula, sua disposição é sinônimo de motivação para eu seguir na carreira universitária. Grato pelo tempo despendido. Boa semana à todos. Rafael Tieppo

1 0

Instalação da lme4
by Andre Oliveira 06 Jul '17

06 Jul '17

Bom dia, alguém está com mesmo problema na instalação da lme4? compilation terminated. make: ** [external.o] Erro 1 ERROR: compilation failed for package ‘lme4’ * removing ‘/usr/local/lib/R/site-library/lme4’ The downloaded source packages are in ‘/tmp/RtmpRFLqvs/downloaded_packages’ Warning message: In install.packages("lme4") : installation of package ‘lme4’ had non-zero exit status > library(lme4) Error in library(lme4) : there is no package called ‘lme4’ Grato André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 6

Equação do modelo
by Bruce Kelly 06 Jul '17

06 Jul '17

Boa tarde listeiros (@s), Estou aplicando a análise fatorial e queria a ajuda para utilizar o pacote 'sem'. De acordo com o Quick-r (http://www.statmethods.net/advstats/factor.html) Segundo o script apresentado no site os fatores F1 e F2 são descritos em função das variáveis X1, X2,...e as cargas lam1,lam2,...ect. Bem minha dúvida é: essas posições das cargas (lam1,lam2,...) são fixas e como posso colocar na equação do modelo as mais significativas, tipo (lam3, lam7). E como posso retirar a equação usando meus dados. Segue minha rotina. af=factanal(dad, factors=8, scores=c('Bartlett'), rotation='varimax') > af Call: factanal(x = dad, factors = 8, scores = c("Bartlett"), rotation = "varimax") Uniquenesses: ï..s3 s7 s11 s15 s16 ad1 ad2 ad3 ad4 ad5 ad9 ad13 ad14 ad15 ad19 0.341 0.005 0.484 0.079 0.005 0.103 0.308 0.098 0.005 0.090 0.108 0.015 0.200 0.005 0.436 Loadings: Factor1 Factor2 Factor3 Factor4 Factor5 Factor6 Factor7 Factor8 ï..s3 0.295 -0.250 0.202 0.140 0.266 0.612 s7 -0.407 0.865 -0.135 0.163 -0.153 s11 -0.426 0.541 0.118 s15 0.930 -0.164 -0.139 s16 0.247 0.561 -0.103 0.771 ad1 0.910 -0.178 0.158 ad2 0.781 -0.229 0.101 -0.101 ad3 0.875 -0.202 0.134 0.248 ad4 -0.134 0.538 -0.267 0.777 ad5 0.887 -0.156 0.295 ad9 0.128 0.926 ad13 0.894 0.258 -0.227 0.235 ad14 0.813 0.183 0.103 0.256 -0.151 ad15 0.389 -0.155 0.397 0.744 0.305 ad19 0.622 0.132 -0.123 0.138 0.149 0.312 Factor1 Factor2 Factor3 Factor4 Factor5 Factor6 Factor7 Factor8 SS loadings 3.832 3.365 1.507 0.997 0.772 0.771 0.751 0.722 Proportion Var 0.255 0.224 0.100 0.066 0.051 0.051 0.050 0.048 Cumulative Var 0.255 0.480 0.580 0.647 0.698 0.750 0.800 0.848 Test of the hypothesis that 8 factors are sufficient. The chi square statistic is 12.7 on 13 degrees of freedom. The p-value is 0.471 Muito obrigado a todos. Bruce Quick-R: Factor Analysis - Quick-R: Home Page<http://www.statmethods.net/advstats/factor.html> www.statmethods.net Learn principal components and factor analysis in R. Factor analysis includes both exploratory and confirmatory methods.

1 0

Bug R in Debian 8,
by Fernando Souza 06 Jul '17

06 Jul '17

Caros colegas Possuo R instalado em uma máquina Debian 8. Recentemente ao utilizar o pacote nome recebi algumas mesgs estranhas de erro. Fui pesquisar e parece estar relacionado a falhas de compilação devido ao período de congelamento de debian8 para passaram para debian9. Já testei e o mesmo acontece com outros pacote s. O erro pare estar relacionado aos computadores c e fortran https://bugs.debian.org/cgi-bin/bugreport.cgi?bug=861333 Minha dúvida é. Como resolver isso? Consegui corrigir o nlme(não sei como) após reinstalar o r e após o nlme. Mas como saber quais são os pacotes problemáticos e como corrigir? As dicas do link acima não funcionaram para mim.

2 1

Random Forest
by Ary Serpa Neto 05 Jul '17

05 Jul '17

2 1

Duvida estatística com R
by Fernando Gama 04 Jul '17

04 Jul '17

Pessoal estou com uma duvida, estou com dificuldades de recordar, tenho duas amostras A e B, para ilustrar: A <- length(216) # o tamanho é variável para cada laço for B <- length(87) Suponhamos que eu tenha um sample sample <- length(15) Gostaria de obter a proporção de acerto considerando o tamanho das amostra B, ou seja: sample (15) / B (216) = 0.1724. Porém, talvez seja interessante (ou correto) levar em consideração o tamanho de A já que A varia de tamanho. E eis que surge a minha dúvida, se eu considerar o tamanho de A para medir a proporção a formula seria: sample (15) / A (216) - B(87) = 0.1162. Correto? Ou teria que considerar apenas a amostra B mesmo? Agradeço a ajuda. -- Att, | Fernando Gama da Mata | | Database Specialist | Master's Degree UFPA | | Contacts: +55 91 99150 0365 | f.fabiogama88(a)gmail.com |

3 3

Anova para amostras de tamanhos diferentes
by Fabiana Goncalves Barbosa 04 Jul '17

04 Jul '17

Olá Pessoal Como faço uma ANOVA para amostras de tamanhos diferentes?

4 3

Re: [R-br] RES: Dúvida sobre pesos amostrais[Resolvido]
by Carlos Andrade 30 Jun '17

30 Jun '17

Cesar, Obrigado pelo esclarecimento. A minha dúvida está resolvida. Em 30 de junho de 2017 21:28, Cesar Rabak via R-br <r-br(a)listas.c3sl.ufpr.br> escreveu: > Carlos, > > Essa dúvida não exatamente pertinente ao R, mas sim uma dúvida de > Estatística Básica. > > No seu caso você obteve uma proporção da preferência por uma variedade de > serviço. > > Não existe o conceito de "extrapolação" em amostragem, mas sim a questão do > intervalo de confiança -- IC dessa estimativa vis a vis com a proporção > "real" da população sob estudo¹. > > Como a sua amostra é relativamente grande com relação à população, o IC > (95%) é relativamente estreito (12,0; 19,8). A forma de calcular esse IC não > é única, razão pela qual pode haver ainda outros ICs dependendo dos costumes > no domínio do problema da ciência que você está praticando. > > No caso das suas questões recolocadas, a sua maneira de dizer que 962 com IC > 95% (727 a 1199) da população têm essa opinião. > > A interpretação do significado desse intervalo e demais questões são > extremamente fora de tópico aqui, por isso se você tiver dúvidas recomendo > que consulte um estatístico na sua instituição, ou encete o resto desta > discução via e-mail em privado². > > Quanto à questão dos pesos ela só entraria em discussão caso a amostragem > não for totalmente aleatória em toda a população, tópico esse que demandaria > também consulta a material mais avançado em Amostragem, ou uma consultoria > com um estatístico aí na sua instituição. > > HTH > -- > Cesar Rabak > > > [1] Quando ela é calculada e apresentada seguindo certo formalismos ela é > denominada "margem de erro" (sim a famosa dos levantamentos de intenção de > votos!). > > [2] Você pode também colocar a sua questão no SE de estatística onde a > questão seria mais bem vinda. > > > > 2017-06-30 13:15 GMT-03:00 Carlos Andrade via R-br > <r-br(a)listas.c3sl.ufpr.br>: >> >> Agradeço a Rodrigo Campos e pjobenil pela colaboração. >> Então, considerando as premissas levantadas por vocês, posso >> seguramente aplicar o percentual da amostra na população. >> Desta forma, 16% (os 51 respondentes) podem ser aplicados aos 6055, e >> dizer que são 962 os que têm aquela opinião? >> >> Em 30 de junho de 2017 11:35, pjobenil--- via R-br >> <r-br(a)listas.c3sl.ufpr.br> escreveu: >> > Carlos, >> > >> > >> > >> > Para os dados apresentados cabe no máximo calcular a margem de erro se >> > usou >> > alguma técnica probabilística para a seleção da amostra. >> > >> > Extrapolar os dados para a população ocorreria em dois casos: >> > >> > Existem segmentos nesta população que você não considerou quando >> > selecionou >> > os 361; >> > Os 6055 indivíduos formam um segmento de uma população maior. >> > >> > Nos dois casos sim você deve aplicar os pesos destes segmentos para o >> > calculo final do estimador. >> > >> > >> > >> > Atenciosamente >> > >> > Jobenil Júnior >> > >> > >> > >> > >> > >> > Enviado do Email para Windows 10 >> > >> > >> > >> > De: Carlos Andrade via R-br >> > Enviado:quinta-feira, 29 de junho de 2017 18:42 >> > Para: R-br(a)listas.c3sl.ufpr.br >> > Assunto: [R-br] Dúvida sobre pesos amostrais >> > >> > >> > >> > Prezados amigos da lista. >> > >> > >> > >> > Desejo obter o esclarecimento de uma dúvida. >> > >> > De uma população de 6055 indivíduos foi obtida uma amostra de 321 >> > >> > respondentes de um questionário. >> > >> > Na tabulação de uma das variáveis encontro que 51(15,9%) dos >> > >> > entrevistados prefere uma variedade de serviço. >> > >> > Dúvida: como faço para extrapolar este valor amostral para o >> > >> > correspondente do universo? >> > >> > É necessário usar pesos amostrais? >> > >> > >> > >> > Agradeço antecipadamente 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 >> > >> > _______________________________________________ >> > >> > 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 fornea >> > cdigo >> > mnimo reproduzvel. >> > >> > >> > >> > >> > _______________________________________________ >> > 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. >> >> >> >> -- >> 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 >> _______________________________________________ >> 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 fornea cdigo >> mnimo reproduzvel. > > > > _______________________________________________ > 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. -- 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 1

Re: [R-br] RES: Dúvida sobre pesos amostrais[Resolvido]
by Carlos Andrade 30 Jun '17

30 Jun '17

Cesar, Obrigado pelo esclarecimento. Minha dúvida está resolvida. Em 30 de junho de 2017 21:28, Cesar Rabak via R-br <r-br(a)listas.c3sl.ufpr.br> escreveu: > Carlos, > > Essa dúvida não exatamente pertinente ao R, mas sim uma dúvida de > Estatística Básica. > > No seu caso você obteve uma proporção da preferência por uma variedade de > serviço. > > Não existe o conceito de "extrapolação" em amostragem, mas sim a questão do > intervalo de confiança -- IC dessa estimativa vis a vis com a proporção > "real" da população sob estudo¹. > > Como a sua amostra é relativamente grande com relação à população, o IC > (95%) é relativamente estreito (12,0; 19,8). A forma de calcular esse IC não > é única, razão pela qual pode haver ainda outros ICs dependendo dos costumes > no domínio do problema da ciência que você está praticando. > > No caso das suas questões recolocadas, a sua maneira de dizer que 962 com IC > 95% (727 a 1199) da população têm essa opinião. > > A interpretação do significado desse intervalo e demais questões são > extremamente fora de tópico aqui, por isso se você tiver dúvidas recomendo > que consulte um estatístico na sua instituição, ou encete o resto desta > discução via e-mail em privado². > > Quanto à questão dos pesos ela só entraria em discussão caso a amostragem > não for totalmente aleatória em toda a população, tópico esse que demandaria > também consulta a material mais avançado em Amostragem, ou uma consultoria > com um estatístico aí na sua instituição. > > HTH > -- > Cesar Rabak > > > [1] Quando ela é calculada e apresentada seguindo certo formalismos ela é > denominada "margem de erro" (sim a famosa dos levantamentos de intenção de > votos!). > > [2] Você pode também colocar a sua questão no SE de estatística onde a > questão seria mais bem vinda. > > > > 2017-06-30 13:15 GMT-03:00 Carlos Andrade via R-br > <r-br(a)listas.c3sl.ufpr.br>: >> >> Agradeço a Rodrigo Campos e pjobenil pela colaboração. >> Então, considerando as premissas levantadas por vocês, posso >> seguramente aplicar o percentual da amostra na população. >> Desta forma, 16% (os 51 respondentes) podem ser aplicados aos 6055, e >> dizer que são 962 os que têm aquela opinião? >> >> Em 30 de junho de 2017 11:35, pjobenil--- via R-br >> <r-br(a)listas.c3sl.ufpr.br> escreveu: >> > Carlos, >> > >> > >> > >> > Para os dados apresentados cabe no máximo calcular a margem de erro se >> > usou >> > alguma técnica probabilística para a seleção da amostra. >> > >> > Extrapolar os dados para a população ocorreria em dois casos: >> > >> > Existem segmentos nesta população que você não considerou quando >> > selecionou >> > os 361; >> > Os 6055 indivíduos formam um segmento de uma população maior. >> > >> > Nos dois casos sim você deve aplicar os pesos destes segmentos para o >> > calculo final do estimador. >> > >> > >> > >> > Atenciosamente >> > >> > Jobenil Júnior >> > >> > >> > >> > >> > >> > Enviado do Email para Windows 10 >> > >> > >> > >> > De: Carlos Andrade via R-br >> > Enviado:quinta-feira, 29 de junho de 2017 18:42 >> > Para: R-br(a)listas.c3sl.ufpr.br >> > Assunto: [R-br] Dúvida sobre pesos amostrais >> > >> > >> > >> > Prezados amigos da lista. >> > >> > >> > >> > Desejo obter o esclarecimento de uma dúvida. >> > >> > De uma população de 6055 indivíduos foi obtida uma amostra de 321 >> > >> > respondentes de um questionário. >> > >> > Na tabulação de uma das variáveis encontro que 51(15,9%) dos >> > >> > entrevistados prefere uma variedade de serviço. >> > >> > Dúvida: como faço para extrapolar este valor amostral para o >> > >> > correspondente do universo? >> > >> > É necessário usar pesos amostrais? >> > >> > >> > >> > Agradeço antecipadamente 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 >> > >> > _______________________________________________ >> > >> > 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 fornea >> > cdigo >> > mnimo reproduzvel. >> > >> > >> > >> > >> > _______________________________________________ >> > 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. >> >> >> >> -- >> 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 >> _______________________________________________ >> 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 fornea cdigo >> mnimo reproduzvel. > > > > _______________________________________________ > 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. -- 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

1 0

Dúvida sobre pesos amostrais
by Carlos Andrade 30 Jun '17

30 Jun '17

Prezados amigos da lista. Desejo obter o esclarecimento de uma dúvida. De uma população de 6055 indivíduos foi obtida uma amostra de 321 respondentes de um questionário. Na tabulação de uma das variáveis encontro que 51(15,9%) dos entrevistados prefere uma variedade de serviço. Dúvida: como faço para extrapolar este valor amostral para o correspondente do universo? É necessário usar pesos amostrais? Agradeço antecipadamente 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

4 4

Cluster
by Emerson Cotta Bodevan 30 Jun '17

30 Jun '17

Oi pessoal, boa noite. Que pacote do R é mais adequado para fazer análise de cluster? Incluindo dendogramas. Algum manual? Ou tutorial? Abraços e obrigado, *Emerson*

2 2

Regressão por polinômios ortogonais
by Ana Beatriz Souza Souza 28 Jun '17

28 Jun '17

________________________________ Boa tarde! Estou reproduzindo um exercício do Livro de Banzatto e Kronka no qual temos um experimento fatorial 3x3 com fatores quantitativos. O teste F para o fator A (quantitativo) não foi significativo e mesmo assim procedeu-se o estudo de regressão. Pergunto: A regressão independe da significância do teste F? Outra coisa, procendo -se a analise de regressão apresento apenas o modelo em que a regressão foi significativa, mesmo o R^2 sendo elevado? Se eu fizer essa análise no R usando o pacote ExpDes, a análise apresentada leva em conta o F? Desde já agradeço Att., Ana Beatriz

2 1

Divulgação IV WORKSHOP DE CIÊNCIA DE DADOS EDUCACIONAIS (WCDE- 2017)
by Jorge Cavalcanti 26 Jun '17

26 Jun '17

Prezados Colegas, Pedimos ajuda com a divulgação deste CFP entre seus colaboradores e estudantes: ---------------------------------------------- *IV WORKSHOP DE CIÊNCIA DE DADOS EDUCACIONAIS (WCDE- 2017)* (antigo WORKSHOP DE MINERAÇÃO DE DADOS EDUCACIONAIS) VI Congresso Brasileiro de Informática em Educação (CBIE- 2017) Recife, 30 de Outubro a 02 de Novembro de 2017 https://sites.google.com/view/wcde/ ------------------------------------------------------------ **Sobre o evento** A análise de dados educacionais, de uma maneira geral, representa uma área de pesquisa emergente em Informática em Educação para o desenvolvimento de métodos que exploram dados oriundos de ambientes educacionais e também administrativos com a finalidade de entender melhor os estudantes e os cenários em que eles aprendem. Dentro deste contexto, nas últimas três edições do Congresso Brasileiro de Informática em Educação (CBIE) esta temática estava sendo parcialmente tratada no Workshop de Mineração de Dados Educacionais (WMDE), por meio de tópicos de interesse que se restringem a dados gerados em ambientes de ensino-aprendizagem e com a aplicação de tarefas de mineração de dados como predição (regressão, séries temporais e classificação), agrupamento ou associação de dados, resultando em descobertas de conhecimento intrínseco nos dados. Contudo, ainda há outras áreas de pesquisa nesta temática como é o caso de Learning Analytics (LA) que faz uso de técnicas de análises de dados,incluindo tarefas de mineração de dados, para confirmar hipóteses colocadas em atividades de ensino que envolvem a aprendizagem do aluno, fomentando assim recursos analíticos ao tutor para entendimento e aprimoramento do ensino-aprendizagem. Por fim, existe também a temática de Academic Analytics (AA) que tem na essência as mesmas abordagens de análise utilizadas em MDE e LA, porém com alteração na fonte de dados, nesse caso advindo de sistemas educacionais administrativos. Como proposta de englobar as três áreas de pesquisa, uma vez que o foco principal de todas estas é analisar dados educacionais para o aprimoramento do ensino e aprendizagem, neste ano se fez a alteração do nome do workshop para Ciência de Dados Educacionais, ampliando o escopo do evento para uma abordagem mais ampla da temática análise de dados e fomentando a Comunidade de Informática em Educação um fórum para atualização, discussão e aplicação de pesquisas que tem como premissa o uso de dados oriundos da educação. O WCDE é parte do Congresso Brasileiro de Informática da Educação (CBIE), um evento anual da SBC, que ocorrerá este ano em Recife entre os dias 30 de Outubro e 02 de Novembro. Informações adicionais sobre o CBIE estão em https://sabertecnologias.com.br/cbie/. **Tópicos de Interesse** A lista abaixo apresenta alguns tópicos iniciais de interesse deste workshop. Porém, ela não é estritamente limitante quanto aos assuntos que podem ser tratados em Ciência de Dados Educacionais. - Modelagem analítica de Bancos de Dados Educacionais e integração com Dados Abertos (Open Data); - Mineração de Dados Educacionais e Personalização da Aprendizagem; - Infra-estrutura para análise de dados educacionais; - Aplicação de análise de dados educacionais no ensino (Superior ou Educação Básica); - Integrando dados educacionais e Internet das Coisas (IoT); - Análise de dados em ambientes emergentes pedagógicos como jogos educacionais, MOOCs e Sistemas Tutores Inteligentes; - Análise de dados na aprendizagem social e colaborativa; - Melhores práticas para adaptar o estado da arte das técnicas analíticas de recuperação de informação, sistemas de recomendação, análise de redes sociais, mineração de opinião para o domínio educativo; - Abordagens, métodos, algoritmos e ferramentas de analítica; - Análise de dados em apoio a tarefas educacionais administrativas; - Uso de dados combinados com teorias para a compreensão da aprendizagem; - Métricas de aprendizagem e avaliação; - Aprendizagem adaptativa e personalizada; - Valores, éticas e leis sobre os dados educacionais; - Visualização de Dados Educacionais; - Processos e metodologias para análises de dados educacionais. **Datas importantes** - Data limite para submissão de artigos: 31/07/2017 - Notificação de artigos aceitos: 10/08/2017 - Entrega da versão final dos artigos (com copyright): 15/08/2017 - Data limite para inscrição de pelo menos um autor por trabalho: 15/09/2017

1 0

(sem assunto)
by Alice Peres 26 Jun '17

26 Jun '17

Bom dia pessoal, estou urgentemente precisando de ajuda! Como é o intervalo de confiança dos coeficientes do modelo de regressão linear múltiplo ? Att

2 1

Ajuda com organização de dados
by Jackson Rodrigues 25 Jun '17

25 Jun '17

Olá pessoal, Preciso muito da ajuda de vocês. Tenho que sintetizar dados de precipitação de vários dias e meses em uma única tabela. No entanto, os dados que possuo estão dispersos em dias e horas diferentes. A maior frequência de observações estão nos dias mais chuvosos enquanto os dias mais secos tem poucas. Vou colocar aqui um exemplo de tabela como os dados estão e em seguida a tabela que gostaria que se tornasse. Os dados estão como na tabela abaixo que representa a precipitaçção de 1 único ano. Vejam que a quantidade de observações varia de um dia para outro o ano todo. Dia Hora Prec Jan 1 7:00 20 1 7:20 15 1 7:30 5 2 6:00 10 2 15:00 16 ... ... ... 31 7:00 2 Fev 1 9:00 6 1 9:45 30 2 6:45 12 2 16:00 20 2 22:00 60 ... ... ... 28 11:30 25 ... ... ... ... Dez 31 ... ... Desse jeito fica difícil trabalhar quando for comparar com a precipitação dos outros anos. Assim gostaria de sintetizar toda essa informação em uma única tabela. Somar tudo que for de Janeiro 1, somar tudo que for de janeiro 2 e assim por diante até o fim do ano gerando a tabela a seguir. Mês Dia-> 1 2 3 4 ... 31 Jan 40 16 50 66 ... 90 Total precipitado por dia Fev 60 120 25 45 ... 87 ... ... ... ... ... ... ... Dez 120 30 26 90 ... 5 A grande questão é como fazê-lo. Preciso de um comando que diga some tudo que for do dia 1 de Jan, some tudo de 2 de Jan e assim sucessivamente. Alguém poderia me ajudar com isso? Obrigado Jackson -

4 3

Re: [R-br] Distribuições Gamma e Inversa Gaussiana
by mauricio.soares＠globo.com 22 Jun '17

22 Jun '17

Texto muito interessante! Obrigado! <quote author='R-br mailing list'> http://www.sbea.org.br/conbea/2014/anais/R0273-1.pdf FCosta Em 21/06/2017 09:26, Maurício Soares via R-br escreveu: > > Bom dia! > > Tenho um conjunto de dados e preciso estimar os parâmetros de uma > Gamma (μ e ν) e também estimar os parâmetros de uma Inversa Gaussiana > (μ e σ^2). > > Nunca tinha mexido nessa parte de distribuições no R. > > Se alguém puder me dar uma luz, agradeço. > > Obrigado! > > Maurício > > > > _______________________________________________ > 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. _______________________________________________ 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. </quote> Quoted from: http://r-br.2285057.n4.nabble.com/R-br-Distribuicoes-Gamma-e-Inversa-Gaussi… _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

1 0

Re: [R-br] Distribuições Gamma e Inversa Gaussiana
by mauricio.soares＠globo.com 22 Jun '17

22 Jun '17

Felipe, obrigado! Eu preciso rodar algum pacote antes de executar essa função? <quote author='R-br mailing list'> Pode usar a função optim: ### Gamma G <- function(par, x) { sum(dgamma(x, shape = par[1], rate = par[2], log = TRUE)) } emv.G <- optim(par = c(1, 10), fn = G, x = dados, ,control=list("fnscale"=-1), hessian=TRUE) emv.G$par ##### standard error sqrt(diag(solve(-emv.G$hessian))) -- Atenciosamente, ================================== Felipe E. Barletta Mendes Estatístico(UFPR) - Conre3 9766-A Mestrando em Bioestatística(UEM) +55 (41)-92077191 +55 (41)-33287216 =================================== _______________________________________________ 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 fornea cdigo mnimo reproduzvel. </quote> Quoted from: http://r-br.2285057.n4.nabble.com/R-br-Distribuicoes-Gamma-e-Inversa-Gaussi… _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

1 0

Somas de Quadrados Tipo III
by Eduardo Bearzoti 22 Jun '17

22 Jun '17

Bom dia, pessoal. Minha dúvida é particularmente sobre o cálculo de somas de quadrados tipo III (mas também do tipo II) no R com conjuntos de dados desbalanceados. Estou ciente da polêmica em torno do uso destas somas de quadrados na comunidade estatística, mas gostaria muito de poder calculá-las no R, e, a não ser que eu esteja enganado, nem sempre o R as calcula corretamente. Pesquisando na rede, descobri que em geral são apontadas duas maneiras de se calcular somas de quadrados tipo III no R: a) Utilizando a função drop1() em modelos ajustados com a especificação de contrastes contr.sum (ou contr.helmert), e b) utilizando a função Anova() (com "A" maiúsculo), do pacote car. Abaixo, eu apresento um exemplo simulado de um ensaio fatorial altamente desbalanceado, faltando várias combinações entre os níveis dos dois fatores. A interação entre eles, por exemplo, apresenta um único grau de liberdade, tamanho o desbalanceamento. dados <- data.frame(y=c(35.44627,43.43353,40.20125,38.1061,39.96211,44.21636,19.78237,23.58636,31.73537, 36.49977,38.71573,41.00456,19.6787,16.43436,32.05407,33.3784,18.58266,17.74051), factorA=as.factor(c(2,2,4,4,4,4,6,6,6,6,6,6,8,8,8,8,10,10)), factorB=as.factor(c(2,2,1.5,1.5,2.5,2.5,1,1,2,2,3,3,1.5,1.5,2.5,2.5,2,2))) Tenho em mente o modelo contendo os fatores A e B, bem como sua interação. No SAS, utilizando o proc glm, obtenho as seguintes somas de quadrados tipo I, tipo II e tipo III: Source DF Type I SS factorA 4 964.4164933 factorB 3 500.0771638 factorA*factorB 1 68.7277880 Source DF Type II SS factorA 3 954.6730970 factorB 3 500.0771638 factorA*factorB 1 68.7277880 Source DF Type III SS factorA 3 954.6730970 factorB 3 500.0771638 factorA*factorB 1 68.7277880 Percebam que neste exemplo as somas de quadrados tipo II e III são iguais, pois as "caselas" apresentam o mesmo número de repetições (duas). No R, tentei reproduzir este resultado seguindo os dois procedimentos citados acima: 1) Utilizando a função drop1() > mod <- lm(y ~ factorA*factorB, contrasts=list(factorA=contr.sum, factorB=contr.sum), data=dados) > drop1(mod, ~.,test = "F") Single term deletions Model: y ~ factorA * factorB Df Sum of Sq RSS AIC F value Pr(>F) <none> 70.841 42.661 factorA 0 0.000 70.841 42.661 factorB 0 0.000 70.841 42.661 factorA:factorB 1 68.728 139.569 52.867 8.7315 0.01609 * --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (apenas a soma de quadrados da interação foi calculada) 2) Utilizando car::Anova() > mod <- lm(y ~ factorA*factorB, data=dados) > library(car) > Anova(mod,type = "II") Anova Table (Type II tests) Response: y Sum Sq Df F value Pr(>F) factorA 1284.10 4 40.7845 8.995e-06 *** factorB 500.08 3 21.1774 0.0002047 *** factorA:factorB 68.73 1 8.7315 0.0160919 * Residuals 70.84 9 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > Anova(mod,type = "III") Error in Anova.III.lm(mod, error, singular.ok = singular.ok, ...) : there are aliased coefficients in the model (soma de quadrados tipo II calculada errada para o fator A, e tipo III não calculadas) Ou seja, não consegui reproduzir nem as somas de quadrados tipo III nem as do tipo II, por nenhum dos dois procedimentos sugeridos na rede. Estou utilizando: > sessionInfo() R version 3.4.0 (2017-04-21) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Ubuntu 14.04.5 LTS (...) other attached packages: [1] car_2.1-4 Eu estaria fazendo alguma coisa equivocada? Qualquer observação, sugestão, ou indicação de alguma maneira alternativa para cálculo das somas de quadrados tipo III no R será muito bem vinda. Abs, Eduardo Bearzoti (UFOP)

5 5

Baixar imagens
by Nicolli Albuquerque 21 Jun '17

21 Jun '17

Estou tentando baixar fotos através da função download.file: foto<-as.character(linksvalidosespecies$link[50]) sp_baixada<-paste(50,".","jpg",sep="") download.file(foto,sp_baixada,mode="wb") Porém, quando tento abrir a foto, o visualizador de imagens diz que não pode abrir o arquivo. Já tentei outros programas e nada deu certo. Alguma alternativa para isso? Obrigada [https://outlook.live.com/owa/?path=/mail/inbox] [https://outlook.live.com/owa/?path=/mail/inbox] [https://outlook.live.com/owa/?path=/mail/inbox] Nicolli Albuquerque de Carvalho _____________________________________________________________ Acadêmica em Engenharia Ambiental e Sanitária - CTEC Membro do Laboratório de Conservação no Século XXI (LACOS XXI) - ICBS Universidade Federal de Alagoas CV Lattes: http://lattes.cnpq.br/8219241762500724 Skype: nicolli.albuquerque [http://www.criarumemail.com/assinatura-de-email-online_files/twitter_circle…]<http://www.criarumemail.com/assinatura-de-email-online.html#> [http://www.criarumemail.com/assinatura-de-email-online_files/facebook_circl…] <https://www.facebook.com/conservacao21/?fref=ts> [http://www.criarumemail.com/assinatura-de-email-online_files/google_circle_…] <http://www.criarumemail.com/assinatura-de-email-online.html#> [http://www.criarumemail.com/assinatura-de-email-online_files/linkedin_circl…] <https://www.linkedin.com/in/nicollialbuquerque/> [http://www.criarumemail.com/assinatura-de-email-online_files/instagram_circ…] <https://www.instagram.com/albuqnicolli/> [http://www.criarumemail.com/assinatura-de-email-online_files/dribbble_circl…] <http://www.criarumemail.com/assinatura-de-email-online.html#> [http://www.criarumemail.com/assinatura-de-email-online_files/youtube_circle…] <http://www.criarumemail.com/assinatura-de-email-online.html#> [http://www.criarumemail.com/assinatura-de-email-online_files/vimeo_circle_c…] <http://www.criarumemail.com/assinatura-de-email-online.html#> [http://www.criarumemail.com/assinatura-de-email-online_files/github_circle_…] <https://github.com/albuquerquenicolli> [http://www.criarumemail.com/assinatura-de-email-online_files/wordpress_circ…] <https://biogeografiadaconservacaoufal.wordpress.com/> Crie sua assinatura online aqui<http://www.criarumemail.com/assinatura-de-email-online.html>

2 1

Distribuições Gamma e Inversa Gaussiana
by Maurício Soares 21 Jun '17

21 Jun '17

Bom dia! Tenho um conjunto de dados e preciso estimar os parâmetros de uma Gamma (μ e ν) e também estimar os parâmetros de uma Inversa Gaussiana (μ e σ^2). Nunca tinha mexido nessa parte de distribuições no R. Se alguém puder me dar uma luz, agradeço. Obrigado! Maurício

3 2

Instalar o pacote BRugs no MAC OS X
by Ana Roberta dos Santos Silva 20 Jun '17

20 Jun '17

Prezado(a), boa tarde! Gostaria de saber como instalar o pacote BRugs no MAC OS X. Quando coloco no R, aparece uma mensagem na qual um dos erros é: OpenBUGS 3.2.2 or later not found. Install it, or specify its location using, for example, R CMD INSTALL BRugs --configure-args='--with-openbugs=/usr/local' No entanto já fiz isso que processo, e mesmo assim não consigo instalar. Alguém sabe como instalar esse pacote no MAC OS X? Segue abaixo a mensagem que aparece no R: > install.packages("BRugs") --- Please select a CRAN mirror for use in this session --- Package which is only available in source form, and may need compilation of C/C++/Fortran: ‘BRugs’ Do you want to attempt to install these from sources? y/n: y installing the source package ‘BRugs’ tentando a URL 'https://vps.fmvz.usp.br/CRAN/src/contrib/BRugs_0.8-6.tar.gz' Content type 'application/x-gzip' length 96175 bytes (93 KB) ================================================== downloaded 93 KB * installing *source* package ‘BRugs’ ... ** package ‘BRugs’ successfully unpacked and MD5 sums checked checking for prefix by checking for OpenBUGS... no configure: error: OpenBUGS 3.2.2 or later not found. Install it, or specify its location using, for example, R CMD INSTALL BRugs --configure-args='--with-openbugs=/usr/local' ERROR: configuration failed for package ‘BRugs’ * removing ‘/Library/Frameworks/R.framework/Versions/3.3/Resources/library/BRugs’ The downloaded source packages are in ‘/private/var/folders/yv/sy0zx8153kscltzr7z9zj7k00000gn/T/Rtmp2b3isg/downloaded_packages’ Warning message: In install.packages("BRugs") : installation of package ‘BRugs’ had non-zero exit status Grata pela atenção! -- Att. Ana Roberta dos Santos Silva

3 3

Criar novo dataframe a partir de uma estrutura
by Elias Carvalho 19 Jun '17

19 Jun '17

Boa tarde pessoal Tenho um dataframe com a seguinte esturutra: data.frame': 1999 obs. of 4 variables: $ RAC: Factor w/ 5 levels "amarela","branca",..: 1 1 1 1 1 1 2 2 2 2 ... $ ECV: Factor w/ 7 levels "casada","desquitada",..: 1 1 6 6 7 1 6 6 6 6 ... $ DCF: Factor w/ 3 levels "não","não sabe",..: NA NA 3 NA NA NA NA NA 3 3 ... $ ERG: Factor w/ 3 levels "ALTO RISCO","BAIXO RISCO",..: 2 1 3 1 1 1 2 2 2 2 ... Gostaria de criar um novo dataframe com o mesmo numero de linhas(1999) mesma quantidade de variáveis(4) e mesmo número de niveis dessas variáveis tipo fator. Achei o comando para criar a estrutura: df2 <- data.frame(matrix("", ncol = ncol(df1), nrow = nrow(df1))) mas não achei como criar com o mesmo tipo de variável e mesma quantidade de níveis, no caso factor -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

3 3

Vetor de Matrizes
by Clodoaldo José Figueredo 19 Jun '17

19 Jun '17

Caros(as) Estou replicando dados e calculando a matriz de covariância em cada reamostragem. Preciso armazenar essas matrizes e posteriormente calcular intervalos de confiança (IC) para cada elemento da matriz. Como faço para armazenar matrizes como num vetor e depois acessar o mesmo elemento de todas elas cada vez que calcular os IC's?? Obrigado pela colaboração Clodoaldo *"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." * Alexander Graham Bell <http://www.pensador.info/frase/NTQwOTE1/> Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348 Métodos Numéricos para Engenharia - Matemática Aplicada Instituto Federal Catarinense - Campus Araquari Rodovia BR 280 - km 27 - Cx. Postal 21 CEP 89245-000 - Araquari/SC Fone: (47) 3803-7240

3 2

Re: [R-br] Digest R-br, volume 78, assunto 8
by sznelwar＠uol.com.br 18 Jun '17

18 Jun '17

1 0

Instalação de pacotes do github
by Andre Oliveira 18 Jun '17

18 Jun '17

Bom dia, alguém saberia dizer sobre o possível motivo deste aviso? Parece ser sobre caminho e não sei resolver! (Ocorre com todos os pacotes vindos do github) Grato devtools::install_github("ropenscilabs/skimr") Error in curl::curl_fetch_disk(url, x$path, handle = handle) : Problem with the SSL CA cert (path? access rights?) André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espírito Santo. IFES

2 1

Re: [R-br] Digest R-br, volume 78, assunto 8
by Diogo Braga Mendes 18 Jun '17

18 Jun '17

Olá Clodoaldo! Se vc quer simular os itens de uma escala Likert com base em uma variável latente com distribuição normal, primeiramente deves informar quais valores desta variável irão limitar os itens da escala. Por exemplo, totalDeColunas <- ncol(norm) totalDeItens<-7 # quantidade de itens na escala # dividirá seus dados em 7 intervalos (com frequências aprox. iguais) faixas<-quantile(norm, seq(0, 1, length.out = totalDeItens+1)) # quantis empíricos norm.likert<-norm for(col in 1:totalDeColunas){ norm.likert[, col] <- cut(norm[, col], breaks = faixas) } Isto não a possibilidade de vc dividir as faixar com base em quantis da distribuição multivariada que vc gerou. É bom checar também modelos da TRI como os modelos de Rasch. ------------------------------------------------ Diogo Braga Mendes - Estatístico do IPAMB-Belém/PA. - Mestre em Estatística-PPGME/UFPA. - Bacharel em Estatística, CONRE 9558. - Fone: (91) (TIM) 98239-7799 / (CLARO) 98416-3098 O âmago da ciência que sempre se apresenta como a mais importante no que refere a resultados práticos, é a pesquisa altamente teórica e abstrata nascida da infatigável curiosidade, exibilidade e força da razão humana - Andrei Sakharov Em 17 de junho de 2017 15:00, <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. Limitando valores numa tabela (Clodoaldo José Figueredo) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 16 Jun 2017 15:35:32 -0300 > From: Clodoaldo José Figueredo > <clodoaldo.figueredo(a)ifc-araquari.edu.br> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Limitando valores numa tabela > Message-ID: > <CAHkeYEYW2wpjm-OUgpkEGQsT5n3frkAh7UjQEhjHV_ > zyRK_r6g(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Caros(as) > > Estou usando o comando mvrnorm para gerar uma certa quantidade de valores > multivariados a partir de um vetor de médias m e uma matriz de covariância > S. > Após estou arredondando os valores usando o comando round > Preciso que os valores criados estejam entre 1 e 7 de forma que fique na > forma de uma escala de Likert. > > Tentei usar o comando abaixo mas não da resultado > > Lim <- 1:600 > for(j in 1:16){ > Lim <- norm[,c(j)] (norm é a matriz (600x16) de valores normal > nultivariados gerados randomicamente) > for (i in 1:600){ > ifelse(Lim[i]<1,1,Lim[i]) > ifelse(Lim[i]>7,7,Lim[i]) > norm[,c(j)] <- Lim > } > } > O que estou fazendo de errado e como consigo fazer essa limitação dos > valores?? > > Obrigado pela ajuda > > Um abraço > > Clodoaldo > > > > *"Que força é esta, eu não sei; tudo o que sei é que existe, e está > disponível apenas quando alguém está num estado em que sabe exatamente o > que quer, e está totalmente determinado a não desistir até conseguir." * > Alexander Graham Bell <http://www.pensador.info/frase/NTQwOTE1/> > > Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348 > Métodos Numéricos para Engenharia - Matemática Aplicada > Instituto Federal Catarinense - Campus Araquari > Rodovia BR 280 - km 27 - Cx. Postal 21 > CEP 89245-000 - Araquari/SC > Fone: (47) 3803-7240 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20170616/b7939552/attachment-0001.html> > > ------------------------------ > > 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 78, assunto 8 > **************************************** >

1 0

Rstudio - Nova versão
by Andre Oliveira 18 Jun '17

18 Jun '17

Bom dia,algo estranho ocorre por aqui. Desde quando atualizei para a nova versão do Rstudio ocorre um problema para abrir o mesmo, carrega carrega e não abri o editor! Alguém passa pelo mesmo problema ou alguém sabe o que pode ser? Estou com Ubuntu 14.04! André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 4

Limitando valores numa tabela
by Clodoaldo José Figueredo 16 Jun '17

16 Jun '17

Caros(as) Estou usando o comando mvrnorm para gerar uma certa quantidade de valores multivariados a partir de um vetor de médias m e uma matriz de covariância S. Após estou arredondando os valores usando o comando round Preciso que os valores criados estejam entre 1 e 7 de forma que fique na forma de uma escala de Likert. Tentei usar o comando abaixo mas não da resultado Lim <- 1:600 for(j in 1:16){ Lim <- norm[,c(j)] (norm é a matriz (600x16) de valores normal nultivariados gerados randomicamente) for (i in 1:600){ ifelse(Lim[i]<1,1,Lim[i]) ifelse(Lim[i]>7,7,Lim[i]) norm[,c(j)] <- Lim } } O que estou fazendo de errado e como consigo fazer essa limitação dos valores?? Obrigado pela ajuda Um abraço Clodoaldo *"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." * Alexander Graham Bell <http://www.pensador.info/frase/NTQwOTE1/> Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348 Métodos Numéricos para Engenharia - Matemática Aplicada Instituto Federal Catarinense - Campus Araquari Rodovia BR 280 - km 27 - Cx. Postal 21 CEP 89245-000 - Araquari/SC Fone: (47) 3803-7240

1 0

Resgate de um valor de um summary
by Clodoaldo José Figueredo 16 Jun '17

16 Jun '17

Caros(as) Estou usando o pacote "sem" para rodar um modelo de equações estruturais e preciso armazenar em uma tabela alguns valores que aparecem nele. O sumário encontra-se abaixo e vou rodar algumas vezes para conseguir criar um intervalo de confiança para alguns valores. Como faço para armazenar, por exemplo, os valores de Goodness-of-fit index RMSEA index summary(cfaOut, conf.level=.95, fit.indices=c("GFI", "RMSEA", "NFI", "CFI", "RNI", "SRMR")) Model Chisquare = 1019.088 Df = 101 Pr(>Chisq) = 1.333823e-151 Goodness-of-fit index = 0.8348212 RMSEA index = 0.1231879 95% CI: (NA, NA) Bentler-Bonett NFI = 0.7700214 Bentler CFI = 0.7870473 Bentler RNI = 0.7870473 SRMR = 0.1150066 No aguardo, obrigado Clodoaldo *"Que força é esta, eu não sei; tudo o que sei é que existe, e está disponível apenas quando alguém está num estado em que sabe exatamente o que quer, e está totalmente determinado a não desistir até conseguir." * Alexander Graham Bell <http://www.pensador.info/frase/NTQwOTE1/> Prof. Clodoaldo José Figueredo Msc - SIAPE 1800348 Métodos Numéricos para Engenharia - Matemática Aplicada Instituto Federal Catarinense - Campus Araquari Rodovia BR 280 - km 27 - Cx. Postal 21 CEP 89245-000 - Araquari/SC Fone: (47) 3803-7240

3 3

Meta Análise utilizando Principal component analysis (PCA)
by Fernando Souza 16 Jun '17

16 Jun '17

Caros colegas É possivel utilizar a técnica PCA para realizar uma meta análise? Se sim como a variância entre estudos é calculada? Encontrei na cran uma pacote chamado MetaPCA que parece realizar tal análise, mas como os autores não citam nenhuma referência teórica fico em dúvida. Alguém tem algum material sobre esse assunto, ou poderia me dar alguma orientação de como proceder? Att -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

1 0

Escolha do valor que mais ocorre
by Elias Carvalho 15 Jun '17

15 Jun '17

Bom dia Preciso escolher em uma série de números, aquele que mais ocorre, só que a série pode variar com dados discretos, contínuos ou dicotômicos. Pensei em usar moda, e vi que no R não há função direta para isso, então eu vi várias sugestões, mas surgiram as seguintes dúvidas: 1. No caso de dados discretos, sem problemas, porém quando for bimodal por exemplo, há um critério para escolha de um dos valores, uma vez que preciso escolher apenas 1 ? Por exemplo, na série tenho "alto", "alto", "médio", "médio", "baixo". A moda seria "alto" e "médio", qual o melhor a se escolher? 2. Para dados contínuos vi que existem fórmulas para calcular a moda, porém isso significa que provavelmente não vou chegar no valores que tenho, então qual o melhor critério para escolher um dos valores ? Escolher o valor mais próximo da moda calculada ? 3. Existe um processo específico para calcular a moda de dados dicotômicos ou é similar a dados discretos ? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 2

CFP 4o. WCDE'2017 - Divulgação
by Jorge Cavalcanti 15 Jun '17

15 Jun '17

Prezados Colegas, Pedimos ajuda com a divulgação deste CFP entre seus colaboradores e estudantes: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ------------------------------------------------------------ ---------------------------------------------- *IV WORKSHOP DE CIÊNCIA DE DADOS EDUCACIONAIS (WCDE’2017)* (antigo WORKSHOP DE MINERAÇÃO DE DADOS EDUCACIONAIS) VI Congresso Brasileiro de Informática em Educação (CBIE’2017) Recife, 30 de Outubro a 02 de Novembro de 2017 https://sites.google.com/view/wcde/ ------------------------------------------------------------ ---------------------------------------------- *Sobre o evento* A análise de dados educacionais, de uma maneira geral, representa uma área de pesquisa emergente em Informática em Educação para o desenvolvimento de métodos que exploram dados oriundos de ambientes educacionais e também administrativos com a finalidade de entender melhor os estudantes e os cenários em que eles aprendem. Dentro deste contexto, nas últimas três edições do Congresso Brasileiro de Informática em Educação (CBIE) esta temática estava sendo parcialmente tratada no Workshop de Mineração de Dados Educacionais (WMDE), por meio de tópicos de interesse que se restringem a dados gerados em ambientes de ensino-aprendizagem e com a aplicação de tarefas de mineração de dados como predição (regressão, séries temporais e classificação), agrupamento ou associação de dados, resultando em descobertas de conhecimento intrínseco nos dados. Contudo, ainda há outras áreas de pesquisa nesta temática como é o caso de Learning Analytics (LA) que faz uso de técnicas de análises de dados,incluindo tarefas de mineração de dados, para confirmar hipóteses colocadas em atividades de ensino que envolvem a aprendizagem do aluno, fomentando assim recursos analíticos ao tutor para entendimento e aprimoramento do ensino-aprendizagem. Por fim, existe também a temática de Academic Analytics (AA) que tem na essência as mesmas abordagens de análise utilizadas em MDE e LA, porém com alteração na fonte de dados, nesse caso advindo de sistemas educacionais administrativos. Como proposta de englobar as três áreas de pesquisa, uma vez que o foco principal de todas estas é analisar dados educacionais para o aprimoramento do ensino e aprendizagem, neste ano se fez a alteração do nome do workshop para Ciência de Dados Educacionais, ampliando o escopo do evento para uma abordagem mais ampla da temática análise de dados e fomentando a Comunidade de Informática em Educação um fórum para atualização, discussão e aplicação de pesquisas que tem como premissa o uso de dados oriundos da educação. O WCDE é parte do Congresso Brasileiro de Informática da Educação (CBIE), um evento anual da SBC, que ocorrerá este ano em Recife entre os dias 30 de Outubro e 02 de Novembro. Informações adicionais sobre o CBIE estão em https://sabertecnologias.com.br/cbie/. *Tópicos de Interesse* A lista abaixo apresenta alguns tópicos iniciais de interesse deste workshop. Porém, ela não é estritamente limitante quanto aos assuntos que podem ser tratados em Ciência de Dados Educacionais. • Modelagem analítica de Bancos de Dados Educacionais e integração com Dados Abertos (Open Data); • Mineração de Dados Educacionais e Personalização da Aprendizagem; • Infra-estrutura para análise de dados educacionais; • Aplicação de análise de dados educacionais no ensino (Superior ou Educação Básica); • Integrando dados educacionais e Internet das Coisas (IoT); • Análise de dados em ambientes emergentes pedagógicos como jogos educacionais, MOOCs e Sistemas Tutores Inteligentes; • Análise de dados na aprendizagem social e colaborativa; • Melhores práticas para adaptar o estado da arte das técnicas analíticas de recuperação de informação, sistemas de recomendação, análise de redes sociais, mineração de opinião para o domínio educativo; • Abordagens, métodos, algoritmos e ferramentas de analítica; • Análise de dados em apoio a tarefas educacionais administrativas; • Uso de dados combinados com teorias para a compreensão da aprendizagem; • Métricas de aprendizagem e avaliação; • Aprendizagem adaptativa e personalizada; • Valores, éticas e leis sobre os dados educacionais; • Visualização de Dados Educacionais; • Processos e metodologias para análises de dados educacionais. *Datas importantes* • Data limite para submissão de artigos: 31/07/2017 • Notificação de artigos aceitos: 10/08/2017 • Entrega da versão final dos artigos (com copyright): 15/08/2017 • Data limite para inscrição de pelo menos um autor por trabalho: 15/09/2017 Outras informações em: https://sites.google.com/view/wcde/ e https://www. facebook.com/groups/296471344146156/ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Agradecemos antecipadamente a colaboração e aproveitamos para pedir desculpas por eventuais recebimentos duplicados desta mensagem. Saudações, Comissão Organizadora

1 0

Re: [R-br] Integral de um modelo não linear para calculo da área abaixo da curva
by sznelwar＠uol.com.br 14 Jun '17

14 Jun '17

1 0

Integral de um modelo não linear para calculo da área abaixo da curva
by ASANTOS 14 Jun '17

14 Jun '17

Caros Membros, Utilizava com frequência o Matlab para resolver integrais por não usar essas operações matemáticas rotineiramente no R. Mas agora resolvi fazer todos os cálculos em um mesmo script e estou tendo problemas para realizar a integral de um modelo não linear para calculo da área abaixo da curva, sendo meu CRM: > #Banco de dados simulado > N=100#tamanho amostral > vI<-c(0.1,0.7,1.3,2,3,4,5,6,7,13) > RES=NULL > vI<-rep(vI,N) > di<-rnorm(N,15) > ap<-rnorm(N,10) > th<-rnorm(N,50) > RES=rbind(Nsim,vI,di,ap,th) > RES_S<-as.data.frame(t(RES)) > # > > #Ajuste de modelo não linear com o pacote minpack.lm > library(minpack.lm) > d1_1_S<- nlsLM(di ~ ap*b0*(1+b1*log(1-b2*vI^b3*th^-b3)), > start=list(b0=1.3333, > b1=0.6667, > b2=-0.3333, > b3=0.3333), > data = RES_S) > summary(d1_1_S) ## Fazendo de conta que os coeficientes são significativos > > #Plotando o gráfico > newdata<-RES_S > newdata$ap <- mean(RES_S$ap) > newdata$th <- mean(RES_S$th) > newdata$yhat <- predict(d1_1_S, newdata, type = "response") > newdata <- newdata[order(newdata$vI), ] > plot(di ~ vI, data = RES_S, xlab="vI",ylab="di" ) > lines(x = newdata$vI, y = newdata$yhat ) > CI<-predict(d1_1_S, interval = "confidence") > #-- Agora eu gostaria de resolver a integral de di ~ ap*b0*(1+b1*log(1-b2*vI^b3*th^-b3), para calcular a área abaixo da curva para meus limites de vI para cada intervalo de 0,1 a 0.7 até 7 a 13. Alguém poderia me dar uma luz? 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/ ======================================================================

2 1

gráfico: fonte itálico nos grupos
by Maurício Lordêlo 10 Jun '17

10 Jun '17

Saudações a todos usuários do R!!! Gostaria de saber como fazer com que os nomes dos grupos (tratamentos) apareçam em itálico em um gráfico. Considerando o cmr abaixo, como fazer para que os nomes das espécies (Alfa e Beta) fiquem em itálico? Maurício resposta = rnorm(20) especie = c(rep("A",10),rep("B",10)) boxplot(resposta~especie,xlab="Espécie",ylab="Variável resposta",names=c("Alfa","Beta"))

3 3

Painel espacial com uso de variavel instrumental
by Edimeire Alexandra Pinto 09 Jun '17

09 Jun '17

Oi gente. Andei vendo que para trabalhar com painel espacial, devemos usar o pacote spml. Estou em duvida sobre como trabalhar com instrumentos usando o comando spgm. Para facilitar, suponha que eu tenha o modelo (resumi bastante só para simplificar, só para dar uma ideia) Y = λ W Y + X1 β+X2 + U onde U = ρ WU + ε suponha que queira usar instrumentos para X1 e WY, como eu faria? Outra coisa, tem como usar um durbin espacial (correlação espacial em X e Y) usando spgm? Obrigada

1 0

Re: [R-br] Digest R-br, volume 78, assunto 2
by Diogo Braga Mendes 09 Jun '17

09 Jun '17

Olá José. Como dito é difícil ajudá-lo sem o código. Mas, a partir das mensagens fornecidas pelo R, deve haver um problema de separabilidade ou quase separabilidade (previsão perfeita) nos seus dados que é caracteriza quando a resposta não varia entre os níveis dos fatores no modelo. Vc pode checar isso com: fn<-function (x) {c(soma = sum(x), obs = length(x))} aggregate(as.numeric(resposta) ~ fator_1+fator_2+...+fator_k, fun = fn, data = seusDados) Provavelmente a coluna *soma *será 0 ou será igual a coluna *obs* para pelo menos 1 cruzamento dos k fatores (qualitativos) em seu modelo. ------------------------------------------------ Diogo Braga Mendes - Estatístico do IPAMB-Belém/PA. - Mestre em Estatística-PPGME/UFPA. - Bacharel em Estatística, CONRE 9558. - Fone: (91) (TIM) 98239-7799 / (CLARO) 98416-3098 O âmago da ciência que sempre se apresenta como a mais importante no que refere a resultados práticos, é a pesquisa altamente teórica e abstrata nascida da infatigável curiosidade, exibilidade e força da razão humana - Andrei Sakharov Em 9 de junho de 2017 12:00, <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: erro ao rodar o modelo de regressão logística no R > (Leonard Assis) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 9 Jun 2017 09:25:04 -0300 > From: Leonard Assis <assis.leonard(a)gmail.com> > To: Jose Fernandes <josefer75(a)yahoo.com.br> > Cc: Estatisticabsb <estatisticabsb(a)yahoogrupos.com.br>, YAHOO_CONRE-3 > <conre3(a)yahoogrupos.com.br>, Matematicos Estatisticos > <stat-math(a)yahoogrupos.com>, R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] erro ao rodar o modelo de regressão logística no > R > Message-ID: > <CAEG0FK9RB+4aHup8xZSxD1yQEVHLCq-DtrpFqAfuPSiCQReYpg@mail. > gmail.com> > Content-Type: text/plain; charset="utf-8" > > José, > > quanto ao seu problema: > > é complicado te ajudar apenas com esta mensagem e sem um código que > possamos reproduzir. qualquer coisa que eu responder, baseado apenas no > warning, será um chute. > > Quanto ao numero de destinatários: > > acho desnecessário expor esta quantidade de gente que você colocou na > conversa. acredito que possa prejudicar até a intenção de respostas.Eu > mesmo recebi esta pergunta em 4 listas diferentes. > > Leonard > > lmassis <at> yahoo <dot> com <dot> br > assis.leonard <at> gmail <dot> com > > 2017-06-09 5:10 GMT-03:00 Jose Fernandes <josefer75(a)yahoo.com.br>: > > > Olá a todos, > > > > Estou a rodar um modelo de regressão logística no R e obtive o seguinte > > erro. > > There were 11 warnings (use warnings() to see them) > > glm.fit: fitted probabilities numerically 0 or 1 occurred > > > > Bom final de semana a todos, > > > > José Moniz Fernandes > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20170609/b6eb29af/attachment-0001.html> > > ------------------------------ > > 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 78, assunto 2 > **************************************** >

1 0

Re: [R-br] erro ao rodar o modelo de regressão logística no R
by Leonard Assis 09 Jun '17

09 Jun '17

José, quanto ao seu problema: é complicado te ajudar apenas com esta mensagem e sem um código que possamos reproduzir. qualquer coisa que eu responder, baseado apenas no warning, será um chute. Quanto ao numero de destinatários: acho desnecessário expor esta quantidade de gente que você colocou na conversa. acredito que possa prejudicar até a intenção de respostas.Eu mesmo recebi esta pergunta em 4 listas diferentes. Leonard lmassis <at> yahoo <dot> com <dot> br assis.leonard <at> gmail <dot> com 2017-06-09 5:10 GMT-03:00 Jose Fernandes <josefer75(a)yahoo.com.br>: > Olá a todos, > > Estou a rodar um modelo de regressão logística no R e obtive o seguinte > erro. > There were 11 warnings (use warnings() to see them) > glm.fit: fitted probabilities numerically 0 or 1 occurred > > Bom final de semana a todos, > > José Moniz Fernandes >

1 0

R + COM - Component Object Model
by Cleber N.Borges 02 Jun '17

02 Jun '17

Olá a todos, Alguém conhece algum pacote ou modos alternativos aos pacotes: 1) RDCOMServer e 2) rcom esses não são mais atualizados. Esses pacotes são para fazer comunicação com outros softwares que disponibilizam uma interface COM - Component Object Model Obrigado Cleber Borges --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

2 1

francimariodelima@gmail.com
by francimariodelima＠gmail.com 02 Jun '17

02 Jun '17

Olá... Estou com um problema na instalação do pacote "simplexreg". Abaixo o erro que estou tendo. > install.packages("simplexreg") Installing package into ‘/home/francimario/R/x86_64-pc-linux-gnu-library/3.4’ (as ‘lib’ is unspecified) tentando a URL 'https://cran.rstudio.com/src/contrib/simplexreg_1.3.tar.gz' Content type 'application/x-gzip' length 30407 bytes (29 KB) ================================================== downloaded 29 KB * installing *source* package ‘simplexreg’ ... ** package ‘simplexreg’ successfully unpacked and MD5 sums checked ** libs g++ -I/usr/share/R/include -DNDEBUG -I/gsl/win64/include -fpic -g -O2 -fstack-protector-strong -Wformat -Werror=format-security -Wdate-time -D_FORTIFY_SOURCE=2 -g -c sim.gee.cpp -o sim.gee.o sim.gee.cpp:5:27: fatal error: gsl/gsl_vector.h: No such file or directory compilation terminated. /usr/lib/R/etc/Makeconf:168: recipe for target 'sim.gee.o' failed make: *** [sim.gee.o] Error 1 ERROR: compilation failed for package ‘simplexreg’ * removing ‘/home/francimario/R/x86_64-pc-linux-gnu-library/3.4/simplexreg’ Warning in install.packages : installation of package ‘simplexreg’ had non-zero exit status The downloaded source packages are in ‘/tmp/RtmpFn17bf/downloaded_packages’ Desde já muito obrigado _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

1 0

Tenho problemas com o rattle(). Quando tento executar no Windows 10, na última versão do R, dá a seguinte mensagem: Error in method(obj, ...) : Elemento raiz invÃ¡lido: "requires" Favor indicar possível solução ppor alguém que tem o mesmo problema.
by Antonio Silva 30 Mai '17

30 Mai '17

O Deducer também tem apresentado problemas de não carregar o Java, bem atualizado. Enviado do Email<https://go.microsoft.com/fwlink/?LinkId=550986> para Windows 10

1 0