- R-br - DINF UFPR

Como utilizar o método da inversa
by charlenny melo 03 Jul '12

03 Jul '12

Boa noite.! Queria saber como utilizo o método da inversa para resolver o seguinte exercício : f(x) = (nx)^n-1 , 0 < x < 1 Alguém pode me ajudar? Obrigado.!

2 1

Apostila, site
by Jodavid Ferreira 03 Jul '12

03 Jul '12

Alguém conhece alguma apostila que, -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

1 0

glm Binomial com somas de sucessos e fracassos diferentes
by Diego Pujoni 02 Jul '12

02 Jul '12

Olá pessoal, Quero fazer um modelo que descreve a probabilidade da presença de uma espécie de peixe em um lago. Tenho algumas caraterísticas dos lagos (variáveis independentes) e relatos (entrevistas) da presença, ou não, do peixe em cada lago (variável dependente). O problema é que não tenho um número fixo de relatos por lago. Ou seja, para alguns lagos, dez pessoas sabiam se tinha a espécie, ou não. Mas para outros lagos somente duas pessoas sabiam. Assim, no modelo glm (family=binomial) entro com o número de sucessos (presença da espécies) e número de fracassos (ausência da espécie), mas a soma de sucessos e fracassos não é fixa. Isto é um problema para a análise? Existe como eu especificar pesos a priori? Muito Obrigado. #Dados# com = c(0,5,3,6,1,9,0,0,4,1,0,0,4,0,2,0,0,0,2) sem = c(9,4,8,3,8,2,9,9,8,10,1,2,0,2,0,2,4,2,0) n.tot = c(340.52,344.07,147.43,411.64,264.75,263.24,1087.00,1028.32,697.57, 428.56,414.46,538.18,519.69,590.99,472.54,381.82,465.55,528.43,334.95) D.L70 = c(4.44,3.92,6.02,3.80,3.23,4.56,8.53,9.56,2.29,3.61,14.99,5.85, 3.96,15.30,3.96,1.01,1.32,1.79,4.51) GLM = glm(cbind(com,sem)~n.tot+D.L70,family = binomial) summary(GLM) Diego PJ

2 1

solidos de revolucao
by Samuel Carvalho 02 Jul '12

02 Jul '12

Boa tarde estimados membros do grupo. Alguem saberia me informar algum material ou até exemplos proprios no R sobre a geraçao de sólidos de revolução? Por exemplo como gerar um solido de revolução para um polinômio de 3º grau ou para uma função generica qualquer? Citei o polinomio de 3º grau por se tratar de um exemplo mais simples mas no meu caso preciso especificamente de algum material que eu consiga gerar volume de sólidos por invólucros cilíndricos Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

2 1

dificuldade para gerar um loop
by mauricio osvaldo moura 02 Jul '12

02 Jul '12

Caros, estou com um problema que deve ser simples de responder mas que eu fiquei empacado. Tenho uma base de dados (data frama) com 15.245 linhas por 369 colunas. Destas 369 colunas as 19 primeiras são médias de variáveis ambientais (retiradas do worldclim) e 330 representam a presená/ausência de determinada espécie nas células. a matrix ficaria assim: Rangebio2 StdDevbio2 CoefVarbio2 Minbio3 Meanbio3 Agalychniscallidryas Agalychnislitodryas Acriscrepitans 2.402 0.034 50 50.889 51.0 0 0 0 1.680 0.022 50 51.611 51.5 0 0 0 1.136 0.014 50 51.344 51.5 0 0 0 0.860 0.011 50 50.867 50.5 0 0 0 2.819 0.045 46 46.971 47.0 0 0 0 2.126 0.032 47 47.639 48.0 0 0 0 O problema é que eu preciso de uma média de cada variável climática para cada espécie, que são 330 em coluna. Minha primeira tentativa foi gerar um classificador com cada espécie e usar o nome da espécie como id e fazer as operações. isso não funcionou. O que eu fiz então, foi usar o seguinte comando: with(spp, mean(Meanbio1 [Agalychnisannae!=0])) Isso funciona bem. O problema é que são 330 nomes e mais 18 variáveis para fazer. tentei fazer um loop para chamar todos os nomes e gerar uma tabela de saida com o nome da espécie (nome da coluna) e a média da variável mas não consegui nada. Por isso gostaria de saber se alguém poderia me dar uma força sugerindo uma forma de fazer esse loop ou indicando outra possibilidade. Como apenas abri o data frame e usei esta função este é o código mínimo (e único) que eu tenho. Obrigado antecipadamente Mauricio

3 2

Comparação de modelos não lineares
by Alisson Lucrecio 02 Jul '12

02 Jul '12

Boa tarde colegas do r-br. Alguem saberia me diz se tem como comparar 3 curvas não lineares ou mais no mesmo teste de razão de verossimilhaça? Obrigado. Alisson Lucrécio da Costa

2 1

Operacionalizando a criação de múltiplos SpatialPolygons
by ASANTOS 02 Jul '12

02 Jul '12

Boa Tarde pessoal, Gostaria de criar múltiplos polígonos do tipo SpatialPolygons sem precisar ficar criando um objeto de cada vez, uma vez que disponho de muitas parcelas. Como por exemplo, tenho 4 polígonos, sendo que as etapas para a criação do SpatialPolygons, fica: ###Transforma as parcelas em objeto espacial parc1 = Polygon(coordV_gr[[1]]) parc2 = Polygon(coordV_gr[[2]]) parc3 = Polygon(coordV_gr[[3]]) parc4 = Polygon(coordV_gr[[4]]) parc.all = Polygons(list(parc1, parc2, parc3, parc4), "parcs") parcSP = SpatialPolygons(list(parc.all)) plot(parcSP) Porém, no meu caso real eu tenho centenas de polígonos (parcelas) dentro do objeto coordV_gr, então gostaria de saber se existe alguma alternativa onde eu especifique o length(coordV_gr) e sejam criados os objetos parc1 ate parc n, lançados dentro do list em Polygons(), criando um novo arquivo parcSP de uma maneira otimizada, ou tenho que fazer parcela por parcela na mão. Segue meu CRM completo abaixo, obrigado, ######################################################################## # Parcelas ##Calculado os vértices e desenhando as parcelas ID<-1:4 ## identificador xcc<-c(672036,672332,672674,672079)##coord X ycc<-c(6783650,6784105,6783470,6782354)#coord Y p.central<-cbind(ID,xcc,ycc) p.central ## Criando um arquivo único com o centro dos transectos ###Criando os vértices do transectos coordV <- vector() for (i in 1:length(p.central)) { coordV <- rbind(coordV,cbind(p.central[,2][i]+c(-20,20,20,-20,-20),p.central[,3][i]+c(20,20,-20,-20,20))) } coordV coordV2=coordV[!is.na(coordV[,1]),] coordV2 ### Criar um objeto list a cada 4 coordenadas coordV_gr<-split(as.data.frame(coordV2), rep(1:nrow(coordV2), each=5, length.out=nrow(coordV2))) coordV_gr ###Transforma as parcelas em objeto espacial parc1 = Polygon(coordV_gr[[1]]) parc2 = Polygon(coordV_gr[[2]]) parc3 = Polygon(coordV_gr[[3]]) parc4 = Polygon(coordV_gr[[4]]) parc.all = Polygons(list(parc1, parc2, parc3, parc4), "parcs") parcSP = SpatialPolygons(list(parc.all)) plot(parcSP) ######################################################################## -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

1 0

teste Skillings–Mack
by João José de Simoni Gouveia 02 Jul '12

02 Jul '12

Prezados, Gostaria da ajuda dos membros desta lista. Necessito analisar o seguinte experimento: Verificar a influência de diferentes concentrações de um aditivo na motilidade do sêmen. A situação é a seguinte, 5 animais tiveram sêmen coletado, alíquotas desse sêmen foram feitas e então a cada alíquota foi adicionada uma concentração de um aditivo (no total de 7 tratamentos). As variáveis analisadas (motilidade e vigor) não permitem análise paramétrica e não estou conseguindo ver a possibilidade de utilizar o teste de Friedman nem o teste Skillings–Mack (implementado no pacote asbio e acho que também no agricolae) pelo fato de que cada animal teve o sêmen coletado mais de uma vez, porém existe desbalanceamento no número de ejaculados/animal, o que cria blocos desbalanceados. Alguém teria uma dica? Att JJ -- ============================================ João José de Simoni Gouveia Médico Veterinário, Msc. Doutorando em Zootecnia - PDIZ (UFC/UFPB/UFRPE) Professor Assistente II - Colegiado Acadêmico de Zootecnia Universidade Federal do Vale do São Francisco (UNIVASF) ---------------------------------------------------------------------------------------- Endereço: UNIVASF (Campus Ciências Agrárias) Rodovia BR 407, Km 12, Lote 543 Projeto de Irrigação Senador Nilo Coelho, s/nº-C1 Petrolina, PE - CEP 56.300-000 Telefone: (87)3986-3800/3804-3801 Skype: joao_jose_de_simoni_gouveia E-mail alternativo: jjsgouveia(a)gmail.com ============================================

2 1

Resultados estratificados
by Edson Lira 02 Jul '12

02 Jul '12

R-Lista, boa tarde. Tenho um banco de dados de avaliação de professores, algumas linhas abaixo. professor q1 q2 q3 q4 q5 q6 q7 q8 q9 q10 1 Jair Martins bom bom bom ótimo ótimo ótimo ótimo ótimo ótimo ótimo 2 Marcos Brandão ótimo ótimo bom ótimo ótimo ótimo ótimo ótimo ótimo ótimo 3 Antonio Filho ótimo ótimo bom bom ótimo ótimo ótimo ótimo ótimo ótimo 4 Luiz Cláudio Leão insuficiente insuficiente insuficiente insuficiente insuficiente insuficiente insuficiente insuficiente insuficiente insuficiente 5 Alexandre Pereira ótimo ótimo ótimo ótimo ótimo ótimo ótimo ótimo ótimo ótimo 6 Marcio André bom ótimo bom bom bom ótimo ótimo ótimo ótimo ótimo q11 q12 1 ótimo ótimo 2 ótimo ótimo 3 ótimo ótimo 4 insuficiente insuficiente 5 ótimo ótimo 6 bom bom Como faço para obter o resultado por professor para todas as questões. []'s. Edson Lira Estatístico Manaus-Amazonas

2 1

Juntando vários objetos dentro de list em loop
by ASANTOS 02 Jul '12

02 Jul '12

Bom dia pessoal, Tento um objeto chamado inter que representa um data frame com as coordenada dos vértices de varias parcelas, sendo que a cada 4 linhas tenho as coordenadas de cada parcela, as quais eu gostaria de armazenar separadamente dentro de um list, mais não tenho conseguido não, alguém poderia me dar um help, suspeito que o rbind() no loop não seja a função adequada, substituí por list e não deu certo também, segue CMR: ### Coordenada X dos vértices das parcelas xV<-c(672016,672056,672056,672016,672312,672352,672352,672312,672654,672694,672694,672654,672059,672099,672099,672059) ### Coordenada Y dos vértices das parcelas yV<-c(6783670,6783670,6783630,6783630,6784125,6784125,6784085,6784085,6783490,6783490,6783450,6783450,6782374,6782374,6782334,6782334) ### cada parcela é formada pelo intervalo de a cada 4 valores e desejo criar um objeto [[]] in list para cada parcela, então: x<-seq(1,length(p.central),4)### Valor inicial do intervalo y<-seq(4,length(p.central),4)### Valor final do intervalo inter<-cbind(x,y) ##Juntado o intervalo a cada quatro linhas ### Criar um objeto list a cada 4 coordenadas coordV_gr <-NULL for(j in 1:(length(inter))){ coordV_gr <- rbind(coordV_gr, list(coordV[inter[,1]:inter[,2],])) } coordV_gr##Resultado Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 3

Função implícita
by Davi 01 Jul '12

01 Jul '12

Bom dia pessoal, O problema que estou é o seguinte: preciso encontrar uma (ou mais, não sei) soluções para uma equação do tipo f(x) - c =0 onde f(x) é uma função implícita e c uma constante qualquer. Tem algum método? Muito obrigado, Davi ___________________________________________________ Davi Butturi-Gomes Doutorando em Estatística e Experimentação Agronômica - ESALQ /USP - Piracicaba Mestre em Biometria - UNESP - Botucatu Ecólogo - UNESP - Rio Claro MSN: davi_bgomes(a)hotmail.com Skype: davibg87 (19)9809.5753

3 4

Função likfit (geoR)
by Rodrigo Sant'Ana 01 Jul '12

01 Jul '12

Caros, estou tentando ajustar um modelo geoestatístico aos meus dados com auxilio da função likfit, presente no pacote geoR, porém a mesma esta ajustando os parâmetros do modelo muito longe (fora) da distribuição observada no variograma dos dados. Achando que poderia ser alguma influência dos valores iniciais, rodei um ajuste com base na minimização dos quadrados, que ajustou bem os parâmetros a distribuição dos dados. Utilizei os valores do ajuste por minimos quadrados como valores iniciais da função de verossimilhança, e ainda sim não ajustou bem! Gostaria da ajuda de vocês para tentar entender o que esta acontecendo. Outro possível resposta que pensei, seria a transformação dos meus dados, pois os mesmos possuem uma certa assimetria. Mas a dúvida que me resultaria, se eu utiliza-se uma transformação de box-cox, como retornaria os valores ao seu normal após a modelagem e predição! []'s CMR temp <- structure(list(x = c(-42.2225631308, -41.4000579262, -41.8113105285, -42.7386448278, -41.2145910664, -43.101514771, -41.2629737255, -43.2869816308, -41.803246752, -42.569305521, -41.9322671762, -42.9966856763, -43.0128132293, -42.7708999338, -41.7790554224, -44.2949536952, -43.4563209376, -43.4805122672, -42.8757290285, -44.1901246005, -43.5692138089, -44.2143159301, -44.5449307672, -44.7465251801, -43.5530862558, -44.4884843316, -43.9885301877, -43.6014689149, -44.1901246005, -44.7949078392, -45.3996910778, -44.8029716157, -45.0932675703, -45.5287115021, -46.0044743165, -45.762561021, -46.6576402142, -46.7060228733, -45.6980508089, -46.3108978241, -46.6011937786, -45.5125839491, -46.6657039908, -46.9882550514, -45.8512625627, -46.1254309642, -46.1496222938, -46.1576860703, -47.1979132408, -46.7624693089, -47.1495305817, -47.0043826044, -47.1575943582, -47.6414209491, -47.9397806802, -47.7623775969, -47.7543138203, -47.0043826044, -47.1172754756, -47.7946327029, -47.076956593, -47.0850203696, -47.7139949378, -47.1092116991, -47.5365918544, -47.1817856877, -48.0446097749, -47.1656581347, -47.9720357863, -48.1494388696, -48.1494388696, -47.9397806802, -48.1575026461, -48.907433862, -49.0364542863, -48.3832883885, -47.9478444567, -49.0041991802, -48.9800078507, -48.9558165211, -47.9881633393, -48.1575026461, -48.9880716272, -48.4558623772, -49.7380028431, -49.7057477371, -49.9315334795, -50.0605539037, -50.0928090098, -50.3831049643, -50.1250641159, -50.1089365628, -50.3185947522, -50.2944034227, -50.4798702825, -50.9959519795, -50.8346764492, -51.6652454303, -51.7055643129, -51.7297556424), y = c(-23.2774714209, -23.3025323734, -23.4194834851, -23.4278371359, -23.4445444376, -23.4445444376, -23.4696053901, -23.4863126917, -23.5030199934, -23.5030199934, -23.5531418984, -23.6951539626, -23.7035076134, -23.7703368201, -23.8037514234, -23.8455196776, -23.8872879318, -23.8956415826, -23.9039952335, -23.9123488843, -23.9457634876, -24.1044828535, -24.129543806, -24.1713120602, -24.1963730127, -24.2047266635, -24.2464949177, -24.279909521, -24.2966168227, -24.9231406353, -24.9565552386, -25.0651526995, -25.2238720653, -25.641554607, -25.6582619087, -25.7334447662, -25.7668593696, -25.8086276237, -25.8169812746, -25.8336885762, -25.8503958779, -25.8587495287, -26.0007615929, -26.0258225454, -26.2346638163, -26.3265539755, -26.3766758805, -26.4685660397, -26.5270415955, -26.7442365172, -26.9113095339, -26.9280168356, -27.0449679472, -27.1786263606, -27.663138109, -27.6798454107, -27.863625729, -27.913747634, -27.913747634, -28.0641133491, -28.0808206507, -28.1977717624, -28.2980155724, -28.523564145, -28.5319177958, -28.5820397008, -28.6405152566, -28.6488689075, -29.1584416084, -29.3004536726, -29.4591730384, -29.5260022451, -29.6680143093, -29.8267336752, -30.1274651052, -30.2026479627, -30.2277089152, -30.2444162169, -30.3947819319, -30.4114892336, -30.4699647894, -30.4699647894, -30.5284403453, -30.6537451078, -31.1215495545, -31.3805127304, -31.4222809846, -31.4389882862, -31.7647806688, -31.8733781296, -31.9986828921, -32.0905730513, -32.2743533697, -32.3495362272, -32.8925235314, -33.519047344, -33.6527057574, -34.0119127433, -34.1455711566, -34.3042905225), biomass = c(83.0961332806, 50.7211614516, 87.2964904063, 94.9368757858, 43.9927789353, 77.2411555034, 48.3943067415, 70.9457447384, 93.3226415451, 111.558449875, 105.120137803, 121.774317746, 121.984738595, 134.194524295, 92.3472981727, 74.0563409989, 118.991685107, 118.913371773, 137.446764012, 95.1474633741, 121.714698897, 134.59315967, 101.797584806, 74.9284080818, 137.771352549, 124.913721538, 160.875499767, 138.782910412, 168.972311629, 120.529672382, 43.8330352522, 111.603138446, 82.6283167931, 102.710993825, 129.466779012, 137.200989394, 71.9554182061, 70.8317598018, 142.071147662, 151.190668965, 99.8896044376, 115.837487832, 118.971176291, 51.6652301103, 193.282510478, 250.005205316, 253.408049368, 250.625852379, 67.7379862984, 216.506386439, 124.554150793, 163.227533825, 137.090380572, 50.1620052836, 28.8243753442, 58.1501629663, 67.1120039497, 177.207650243, 189.334832544, 65.2785306408, 182.688259577, 177.288638618, 90.8691963808, 145.065202159, 127.324867544, 142.685381792, 64.4063001299, 131.548465882, 156.065568878, 169.403052205, 205.245173124, 228.321987613, 242.44647784, 40.1581317266, 43.3583613984, 197.084041126, 214.893149963, 54.8643340973, 66.3804698828, 70.0325520737, 163.330903327, 179.736007179, 70.2135583494, 123.133087815, 53.2497546853, 63.4470098487, 65.9584005227, 62.0232097497, 77.2506169312, 61.9183653079, 82.7548476928, 83.9715064102, 87.0964042733, 91.4666549661, 119.213889345, 86.1920719644, 89.9531669855, 55.7705625113, 52.9182431297, 47.4292707838)), .Names = c("x", "y", "biomass" ), row.names = c(NA, -100L), class = "data.frame") bor <- structure(list(coords.x1 = c(-52.0676766304347, -52.0039483695653, -51.956152173913, -51.7968315217391, -51.6853070652174, -51.6853070652174, -51.5419184782609, -51.4144619565217, -51.2710733695652, -51.1914130434782, -51.0002282608695, -50.8090434782609, -50.7134510869565, -50.6178586956521, -50.5541304347826, -50.5541304347826, -50.5700624999999, -50.4744701086956, -50.4585380434783, -50.410741847826, -50.4426059782608, -50.4426059782608, -50.4585380434783, -50.3310815217391, -50.2195570652173, -50.0283722826087, -49.8531195652174, -49.6619347826086, -49.47075, -49.3751576086956, -49.2636331521738, -49.1202445652174, -48.9768559782608, -48.8971956521739, -48.8334673913043, -48.7060108695651, -48.6422826086957, -48.658214673913, -48.4670298913043, -48.3077092391304, -48.1643206521739, -48.0368641304347, -48.0687282608695, -48.0209320652174, -48.0368641304347, -47.9890679347826, -47.9253396739131, -47.8297472826087, -47.7341548913043, -47.7341548913043, -47.6066983695652, -47.4951739130435, -47.4633097826087, -47.4473777173913, -47.3039891304348, -47.2083967391304, -47.0331440217391, -46.9216195652173, -46.7304347826087, -46.5073858695652, -46.2524728260869, -45.9816277173912, -45.7426467391305, -45.4558695652173, -45.2965489130434, -45.1212961956522, -44.9619755434783, -44.8504510869565, -44.7070624999999, -44.5955380434783, -44.4202853260869, -44.0857119565217, -43.894527173913, -43.6555461956522, -43.3847010869565, -43.0341956521739, -42.795214673913, -42.6040298913043, -42.4128451086956, -42.221660326087, -42.0464076086956, -41.9348831521739, -41.7436983695652, -41.4728532608696, -41.2976005434782, -41.1542119565217, -41.0108233695652, -41.0426875, -41.1542119565217, -41.2816684782609, -41.329464673913, -41.4887853260869, -41.6640380434783, -41.8233586956522, -41.9986114130435, -42.2535244565217, -42.4447092391304, -42.5243695652173, -42.7792826086957, -42.9704673913043, -43.3050407608695, -43.4802934782608, -43.6396141304347, -43.8626630434782, -44.1016440217391, -44.2928288043478, -44.4840135869565, -44.627402173913, -44.7548586956522, -44.8982472826086, -45.0097717391304, -45.1372282608695, -45.2646847826086, -45.3602771739131, -45.4558695652173, -45.5992581521739, -45.5514619565217, -45.5992581521739, -45.7426467391305, -45.9816277173912, -46.1728125, -46.3162010869565, -46.5073858695652, -46.7145027173913, -46.8419592391304, -46.9375516304348, -46.9375516304348, -47.0331440217391, -47.0331440217391, -47.0012798913043, -46.9056875, -46.8578913043478, -46.8419592391304, -46.9056875, -46.985347826087, -47.0650081521739, -47.2243288043478, -47.4155135869565, -47.6066983695652, -47.7341548913043, -47.7819510869565, -47.8297472826087, -47.7819510869565, -47.8138152173913, -47.8616114130435, -47.9253396739131, -47.9253396739131, -48.0687282608695, -48.291777173913, -48.4829619565218, -48.7219429347826, -48.8175353260869, -48.9927880434782, -49.1680407608695, -49.3114293478261, -49.5344782608695, -49.6778668478261, -49.7734592391304, -49.9487119565217, -50.0921005434783, -50.1398967391304, -50.2673532608696, -50.2832853260869, -50.2832853260869, -50.3151494565217, -50.3629456521739, -50.5541304347826, -50.6815869565217, -50.8727717391304, -51.0002282608695, -51.2073451086957, -51.3029375, -51.3348016304348, -51.4463260869564, -51.3507336956521, -51.4941222826087, -51.5737826086956, -51.6215788043478, -51.7808994565218, -51.9083559782608, -52.0676766304347), coords.x2 = c(-34.511214673913, -34.3200298913043, -34.2244375, -34.1607092391304, -34.0651168478261, -33.9217282608695, -33.7464755434783, -33.6349510869565, -33.5074945652174, -33.3800380434783, -33.2685135869565, -33.0932608695652, -33.1410570652174, -33.0613967391304, -32.9020760869565, -32.7427554347826, -32.5675027173913, -32.471910326087, -32.2647934782608, -32.0736086956522, -31.8824239130434, -31.7390353260869, -31.5956467391304, -31.372597826087, -31.1495489130434, -30.9583641304347, -30.7193831521739, -30.5600625, -30.4644701086957, -30.2732853260869, -30.0980326086956, -29.9865081521739, -29.7953233695652, -29.6041385869565, -29.4288858695652, -29.3332934782608, -29.1739728260869, -29.0305842391305, -28.8393994565217, -28.648214673913, -28.4092336956522, -28.281777173913, -28.1224565217391, -28.0268641304347, -27.8834755434783, -27.7082228260869, -27.5966983695652, -27.3895815217391, -27.2302608695651, -27.0709402173912, -26.9275516304348, -26.7363668478261, -26.6248423913044, -26.5133179347825, -26.3858614130435, -26.1946766304347, -25.9875597826087, -25.8441711956521, -25.7007826086956, -25.5573940217391, -25.3662092391304, -25.2387527173913, -25.1272282608695, -24.9997717391304, -24.8563831521739, -24.6651983695652, -24.4580815217391, -24.250964673913, -24.1394402173913, -24.0597798913043, -23.8845271739131, -23.7570706521739, -23.6614782608695, -23.5340217391304, -23.4384293478261, -23.3428369565217, -23.2950407608695, -23.2153804347826, -23.1675842391305, -23.1038559782609, -23.0560597826087, -23.0560597826087, -23.1038559782609, -23.2153804347826, -23.2153804347826, -23.2153804347826, -23.1994483695652, -23.3109728260869, -23.4702934782608, -23.5658858695652, -23.6774103260869, -23.8048668478261, -23.8207989130435, -23.8207989130435, -23.7730027173913, -23.7730027173913, -23.7730027173913, -23.8685951086956, -23.9641875, -24.0597798913043, -24.2031684782609, -24.2828288043478, -24.3465570652173, -24.4262173913044, -24.5377418478261, -24.7129945652174, -24.8404510869565, -24.9519755434783, -25.143160326087, -25.3343451086956, -25.509597826087, -25.7007826086956, -25.796375, -25.8760353260869, -25.8441711956521, -25.9716277173912, -26.0353559782608, -26.1309483695651, -26.1787445652174, -26.5133179347825, -26.8000951086956, -27.0072119565218, -27.1346684782609, -27.182464673913, -27.2461929347825, -27.3417853260869, -27.4214456521739, -27.4851739130435, -27.5807663043478, -27.7560190217391, -27.915339673913, -28.0905923913044, -28.2499130434783, -28.3933016304348, -28.5844864130434, -28.7119429347826, -28.7119429347826, -28.807535326087, -28.9031277173913, -29.0943124999999, -29.3014293478261, -29.4766820652174, -29.7475271739131, -29.9865081521739, -30.2254891304348, -30.3688777173912, -30.5122663043478, -30.5122663043478, -30.6078586956521, -30.6875190217391, -30.7353152173913, -30.7353152173913, -30.7831114130435, -30.9105679347826, -31.0380244565217, -31.2292092391304, -31.4044619565217, -31.6115788043478, -31.8346277173912, -32.1054728260869, -32.3125896739131, -32.4878423913044, -32.6949592391304, -32.8861440217391, -33.125125, -33.3003777173913, -33.4915624999999, -33.6030869565217, -33.6827472826087, -33.6986793478261, -33.7942717391304, -33.7942717391304, -33.9057961956522, -33.9057961956522, -34.0173206521739, -34.0969809782609, -34.1607092391304, -34.2881657608695, -34.4474864130434, -34.4793505434782, -34.511214673913 )), .Names = c("coords.x1", "coords.x2"), class = "data.frame", row.names = c(NA, -181L)) dados.geo <- list() dados.geo$coords <- cbind(temp$long, temp$lat) dados.geo$data <- temp$biomass dados.geo$borders <- bor class(dados.geo) <- "geodata" plot(dados.geo, lowess = TRUE) vario.exp <- variog(dados.geo, uvec = seq(0, 15, l = 100), estimator.type = "classical", max.dist = 15) ## Máxima verossimilhança - Não está convergindo corretamente!!! mod.exp <- likfit(dados.geo, ini.cov.pars = c(2354, 1.5), cov.model = "exp", trend = "cte", fix.nugget = FALSE, nugget = 0.1, lik.method = "ML") ## Minimos quadrados - OK! mod.exp <- variofit(vario.exp, ini = c(3000, 0.3), fix.nugget = FALSE, nugget = 0.1, cov.model = "exp", max.dist = 10) _______________ Rodrigo Sant'Ana Oceanógrafo

2 2

Utilizar links: log-log, log-complement e odds power para familia binomial em MLG.
by andrebvs 01 Jul '12

01 Jul '12

2 2

[Dúvida] BFGS - Chutes iniciais
by Pedro Rafael 01 Jul '12

01 Jul '12

Pessoal como escolher os chutes iniciais do Método BFGS. Sei que teoricamente os chutes podem ser quaisquer. Estou querendo utilizar o método BFGS para maximizar a função de log-verossimilança da distribuição beta. -- [ ], Pedro Rafael Diniz Marinho.

2 2

Como gerar dados de um experimento
by Edmundo Caetano 30 Jun '12

30 Jun '12

Prezados, Pretendo gerar dados de um Experimento em Quadrados Latinos com dependência espacial no R. Gostaria de saber como devo proceder para obter esses dados? Saudações. Edmundo Caetano "Ninguém é tão sábio que nada tenha para aprender, nem tão tolo que nada tenha para ensinar." Blaise Pascal ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quinta-feira, 28 de Junho de 2012 12:00 Assunto: Digest R-br, volume 16, assunto 35 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: manipulação de dados (Walmes Zeviani) 2. Re: manipulação de dados (Natalia Martins) 3. Re: manipulação de dados (Marcos Silva) 4. Re: manipulação de dados (Paulo Justiniano) 5. Re: Colorir alguns municípios em um mapa estadual com cores distintas (Wirton Macedo Coutinho) 6. Re: manipulação de dados (Walmes Zeviani) 7. Unrar (Diogo Ferrari) 8. Re: Unrar (Paulo Justiniano) 9. (sem assunto) (Taynãna César Simões) 10. (sem assunto) (Taynãna César Simões) ---------------------------------------------------------------------- Message: 1 Date: Wed, 27 Jun 2012 12:11:48 -0300 From: Walmes Zeviani <walmeszeviani(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] manipulação de dados Message-ID: <CAFU=EkazDuOZPXHiVZc8bODhEktVK3EkRGd04pq0-kQYrdOTEg(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Além da sugestão do PJ, você pode fazer isso com o comando stack(), reshape::melt() e plyr::adply(). Um exemplo disponível em http://ridiculas.wordpress.com/2011/05/22/converter-dados-em-formato-amplo-… À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120627/b67d6eae/atta…> ------------------------------ Message: 2 Date: Wed, 27 Jun 2012 12:41:40 -0300 From: Natalia Martins <nsmbarreto(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] manipulação de dados Message-ID: <CAGPGL2E2tR=xaBfH+njRW5Ft3Fw9ZRH9kY_5=-XvVya1nS24rQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Walmes mais eu nao tenho a variavel de identificação para nenhum dos casos, devo cria-la? Em 27 de junho de 2012 12:11, Walmes Zeviani <walmeszeviani(a)gmail.com>escreveu: > Além da sugestão do PJ, você pode fazer isso com o comando stack(), > reshape::melt() e plyr::adply(). Um exemplo disponível em > > > http://ridiculas.wordpress.com/2011/05/22/converter-dados-em-formato-amplo-… > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > _______________________________________________ > 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. > -- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120627/3dae60da/atta…> ------------------------------ Message: 3 Date: Wed, 27 Jun 2012 13:17:51 -0300 From: Marcos Silva <marcosfs2006(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] manipulação de dados Message-ID: <CAHKdobz_a9W8vQ_9mATPLDiHdbqcpCKcKMJ=cMuppBQLDBR9kw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Veja se a função t() para transpor os dados te ajuda. Abs. Em 27 de junho de 2012 12:41, Natalia Martins <nsmbarreto(a)gmail.com>escreveu: > Walmes mais eu nao tenho a variavel de identificação para nenhum dos > casos, devo cria-la? > > Em 27 de junho de 2012 12:11, Walmes Zeviani <walmeszeviani(a)gmail.com>escreveu: > >> Além da sugestão do PJ, você pode fazer isso com o comando stack(), >> reshape::melt() e plyr::adply(). Um exemplo disponível em >> >> >> http://ridiculas.wordpress.com/2011/05/22/converter-dados-em-formato-amplo-… >> >> À disposição. >> Walmes. >> >> ========================================================================== >> Walmes Marques Zeviani >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) >> Departamento de Estatística - Universidade Federal do Paraná >> fone: (+55) 41 3361 3573 >> VoIP: (3361 3600) 1053 1173 >> e-mail: walmes(a)ufpr.br >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> >> _______________________________________________ >> 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. >> > > > > -- > > Natália da Silva Martins > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > Contato: (19) 8306-4743 > > > > _______________________________________________ > 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. > -- Marcos F. Silva http://sites.google.com/site/marcosfs2006 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120627/ebbbd927/atta…> ------------------------------ Message: 4 Date: Wed, 27 Jun 2012 14:03:34 -0300 (BRT) From: Paulo Justiniano <paulojus(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] manipulação de dados Message-ID: <alpine.DEB.2.00.1206271402020.9577(a)pataxo.est.ufpr.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" Natalia se voce tiver uma variável apenas pode fazer algo masi simples ## simulando dados com estrutura similar M <- matrix(rpois(5*8, lam=8), nc=5) colnames(M) <- paste("D", 1:5, sep="") M o que voce quer é o seguinte: ## passando apra data-frame empilhandos as linhas de M df <- data.frame(data=colnames(M), dados=as.vector(t(M))) df On Wed, 27 Jun 2012, Natalia Martins wrote: > Walmes mais eu nao tenho a variavel de identificação para nenhum dos casos, devo cria-la? > > Em 27 de junho de 2012 12:11, Walmes Zeviani <walmeszeviani(a)gmail.com> escreveu: > Além da sugestão do PJ, você pode fazer isso com o comando stack(), reshape::melt() e plyr::adply(). Um exemplo disponível > em > > http://ridiculas.wordpress.com/2011/05/22/converter-dados-em-formato-amplo-… > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > _______________________________________________ > 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. > > > > > -- > > Natália da Silva Martins > Bacharel em Estatística - Universidade Estadual de Maringá/ UEM > Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP > Contato: (19) 8306-4743 > > > > ------------------------------ Message: 5 Date: Wed, 27 Jun 2012 12:54:19 -0700 (PDT) From: Wirton Macedo Coutinho <wirton_coutinho(a)yahoo.com.br> To: Augusto Ribas <ribas.aca(a)gmail.com>, Paulo Justiniano <paulojus(a)leg.ufpr.br>, "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Colorir alguns municípios em um mapa estadual com cores distintas Message-ID: <1340826859.51307.YahooMailNeo(a)web163805.mail.gq1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Boa tarde a todos, Com as dicas do Augusto Ribas e do prof. Paulo Justiniano, consegui, por meio de um post na internet, resolver o meu problema de colorir com cores distintas alguns municípios do Estado de Goiásem um mapa com extensão shp. O código final, a quem possa interessar, ficou assim: require(maptools) map <- readShapePoly(file.choose()) summary(map) map$Nome_Munic colors <- c("lightskyblue3", "antiquewhite1", "darksalmon", "cornflowerblue", "gray94", "darkkhaki", "darkslategray3", "bisque2") counties2 <- c("Acreúna", "Chapadão do Céu", "Cristalina", "Montividiu", "Palmeiras de Goiás", "Santa Helena de Goiás", "Silvânia") counties.all = map$Nome_Munic n.all = length(counties.all) col.map = numeric(n.all) for (i in 1:length(counties2)){ col.map[grep(counties2[i],counties.all)] = runif(1) } c2 = col.map!=0 col.map2=col.map col.map2[!c2]= gray(0.94) col.map2[c2]=colors plot(map,col=col.map2,axes=T) legend("topleft", legend=c("Acreúna", "Chapadão do Céu", "Cristalina", "Montividiu", "Palmeiras de Goiás", "Santa Helena de Goiás", "Silvânia"), cex=0.9, bty="n", fill= c("lightskyblue3", "antiquewhite1", "darksalmon", "cornflowerblue", "darkkhaki", "darkslategray3", "bisque2")) Obrigado a todos, Att., ----- Wirton Macedo Coutinho Pesquisador - Fitopatologia Embrapa Algodão Rua Oswaldo Cruz, 1143, Centenário Campina Grande, PB CEP 58429-095 >________________________________ > De: Augusto Ribas <ribas.aca(a)gmail.com> >Para: r-br(a)listas.c3sl.ufpr.br; Wirton Macedo Coutinho <wirton_coutinho(a)yahoo.com.br> >Enviadas: Segunda-feira, 25 de Junho de 2012 17:27 >Assunto: Re: [R-br] Colorir alguns municípios em um mapa estadual com cores distintas > >http://r-br.2285057.n4.nabble.com/R-br-Apostila-de-R-Jackson-Aquino-td44351… > >Olhe esse post antigo. Tem um link com uma apostila muito boa pra ter >uma ideia de como mexer com mapas, o basico. >Capitulo 11 mata sua duvida e muito mais. >Em portugues, muito bem escrito, uma sugestão de leitura que vi na >lista aqui e me ajudou muito. > >Em 25 de junho de 2012 16:14, Augusto Ribas <ribas.aca(a)gmail.com> escreveu: >> Bom, seu problema é mais ou menos isso. >> >> #abrindo pacote pra plotar mapas, assumo que vc ta usando esse, o maptools >> library(maptools) >> >> #lendo um mapa pra exemplo >> mapa<- readShapeSpatial(system.file("shapes/sids.shp", package="maptools")[1], >> IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66")) >> >> #mapa >> summary(mapa) >> #informações aqui >> mapa@data >> #vamo pega o nome das cidades >> mapa@data$NAME >> >> >> #o problema que vc ta fazendo isso >> plot(mapa, axes=TRUE, col=mapa@data$NAME=="Alamance") >> plot(mapa, col=mapa@data$NAME=="Avery", add=TRUE) >> plot(mapa, col=mapa@data$NAME=="Buncombe", add=TRUE) >> >> #no col vc ta pondo essa informação >> mapa@data$NAME=="Alamance" >> #que é sempre 1, logo vc colore de preto, 1 na palleta default é >> preto, true é igual 1, tem 0 0 0 que é branco e 1 que é preto >> sum(mapa@data$NAME=="Alamance") >> >> #vc precisa de numeros diferentes pra cores >> #ou o nome da cor pra cidade >> #ou faze uma palleta >> #mas assim vc tem numeros diferentes >> as.numeric(mapa@data$NAME) >> >> >> #e vai sair cores diferentes, olha esse plot, note o que to usando no col. >> plot(mapa, axes=TRUE, col=as.numeric(mapa@data$NAME)) >> >> #nesse exemplo vc ta usando essas cores (numeros) pra esses lugares >> data.frame(cores=as.numeric(mapa@data$NAME),localidade=mapa@data$NAME) >> >> #então faz um vetor, com o numero de cidades, e escrever "black", >> "blue", as cores que vc quer por cidade, ou usa numeros, ou usa uma >> palleta pronta, tipo as que vem no pacote Rcolorbrewer. >> >> Bem espero ter ajudado. >> >> Em 25 de junho de 2012 13:00, Wirton Macedo Coutinho >> <wirton_coutinho(a)yahoo.com.br> escreveu: >>> Boa tarde pessoal, >>> >>> Gostaria de saber se alguém tem alguma dica de como colorir alguns >>> municipios em um mapa estadual (extensão shp) com cores distintas. >>> >>> Preciso colorir 7 municipios no mapa do Estado de Goiás com cores distintas >>> (cada município com uma cor diferente), mas não estou conseguindo. >>> >>> Tenho tentado, sem sucesso, da seguinte forma: >>> >>> plot(mapa, axes=TRUE, col=mapa$Nome_Munic=="Santa Helena de Goiás") >>> plot(mapa, col=mapa$Nome_Munic=="Chapadão do Céu", add=TRUE) >>> plot(mapa, col=mapa$Nome_Munic=="Montividiu", add=TRUE) >>> plot(mapa, col=mapa$Nome_Munic=="Acreúna", add=TRUE) >>> plot(mapa, col=mapa$Nome_Munic=="Perolândia", add=TRUE) >>> plot(mapa, col=mapa$Nome_Munic=="Palmeiras de Goiás", add=TRUE) >>> plot(mapa, col=mapa$Nome_Munic=="Cristalina", add=TRUE) >>> >>> Apesar de conseguir colorir todos os 7 municipios, com os comando acima, >>> todos os 7 ficam com a mesma cor. >>> >>> Certo de contar com a boa vontade de todos, agradeço antecipadamente. >>> >>> Att., >>> >>> -- >>> Wirton Macedo Coutinho >>> Pesquisador Fitopatologia >>> Embrapa Algodão >>> Rua Oswaldo Cruz, 1143, Centenário >>> Campina Grande PB >>> CEP 28428-095 >>> >>> _______________________________________________ >>> 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. >> >> >> >> -- >> Grato >> Augusto C. A. Ribas >> >> Site Pessoal: http://augustoribas.heliohost.org >> Lattes: http://lattes.cnpq.br/7355685961127056 > > > >-- >Grato >Augusto C. A. Ribas > >Site Pessoal: http://augustoribas.heliohost.org >Lattes: http://lattes.cnpq.br/7355685961127056 > > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120627/4bde93f3/atta…> ------------------------------ Message: 6 Date: Wed, 27 Jun 2012 17:52:22 -0300 From: Walmes Zeviani <walmeszeviani(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] manipulação de dados Message-ID: <CAFU=EkYGBxeOYi4k6SVRryfepe09TJjH2DBzG8O0OV51iGuiLQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Sim, Possivelmente a sua identificação é o próprio número da linha, então você pode fazer dados$id <- 1:nrow(dados) e aí segue conforme a matéria do ridículas. À dispposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120627/34b2d3d2/atta…> ------------------------------ Message: 7 Date: Wed, 27 Jun 2012 18:57:54 -0300 From: Diogo Ferrari <diogoferrari(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Unrar Message-ID: <CAEvjS8+FxoBy2CihmgrLtdxTuy=TW9VT0zYZOzKV4uESuaRYWA(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Caros, Alguem sabe como descompactar arquivos .rar usando R. Valeu -- Diogo Ferrari "A vida é a arte de tirar conclusões suficientes de dados insuficientes" -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120627/4524123d/atta…> ------------------------------ Message: 8 Date: Wed, 27 Jun 2012 20:52:32 -0300 (BRT) From: Paulo Justiniano <paulojus(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Unrar Message-ID: <alpine.DEB.2.00.1206272051520.18519(a)pataxo.est.ufpr.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" de forma geral o comoando system() permite voce rodar comandos do seu sistema algo como system("unrar ETC ETC") On Wed, 27 Jun 2012, Diogo Ferrari wrote: > Caros, > Alguem sabe como descompactar arquivos .rar usando R. > > Valeu > > -- > Diogo Ferrari > > "A vida é a arte de tirar conclusões suficientes de dados insuficientes" > > ------------------------------ Message: 9 Date: Wed, 27 Jun 2012 18:09:21 -0700 (PDT) From: Taynãna César Simões <tay_estatistica(a)yahoo.com.br> To: kosiq(a)msn.com, luciane.velasque(a)uniriotec.br, fabiodemarqui(a)zipmail.com.br, leolataliza(a)hotmail.com, r_bruniera(a)yahoo.com.br, r-br(a)listas.c3sl.ufpr.br, tay_estatistica(a)yahoo.com, staelhorta(a)terra.com.br, pattyarantes(a)ibest.com.br Subject: [R-br] (sem assunto) Message-ID: <1340845761.48543.YahooMailNeo(a)web33905.mail.mud.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" http://beacheschurch.org/xmlrpc/googlesave.html?grtyu=at.gio&wrg=te.hkml&oe… -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120627/45500caa/atta…> ------------------------------ Message: 10 Date: Wed, 27 Jun 2012 19:36:23 -0700 (PDT) From: Taynãna César Simões <tay_estatistica(a)yahoo.com.br> To: ana.paula.domingues(a)bancoreal.com.br, saint_clair_jr(a)uol.com.br, alessandracsp(a)gmail.com, r-br(a)listas.c3sl.ufpr.br, natasha_feigel(a)yahoo.com.br, tgpimentel(a)gmail.com, carolina_okumura(a)yahoo.com.br, juninholana(a)yahoo.com.br, constantino(a)cpqam.fiocruz.br Subject: [R-br] (sem assunto) Message-ID: <1340850983.41807.YahooMailNeo(a)web33905.mail.mud.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" http://huntfordbcooper.com/wp-content/themes/sandbox/googlesave.html?ees=xs… -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120627/5eba0dfb/atta…> ------------------------------ _______________________________________________ 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 16, assunto 35 *****************************************

3 3

(sem assunto)
by Taynãna César Simões 29 Jun '12

29 Jun '12

http://huntfordbcooper.com/wp-content/themes/sandbox/googlesave.html?ees=xs…

2 1

Leitura de malhas digitais de 2010
by Lucas Barbosa 29 Jun '12

29 Jun '12

Boa noite, colegas! Estou tendo as mensagens de erro abaixo ao tentar ler as malhas digitais do Censo 2010. Eu sempre uso essas mesmas funções para importar a malha municipal de 2000, então a malha nova é que está com problema? Links para as malhas: ftp://geoftp.ibge.gov.br/malhas_digitais/censo_2010/setores_censitarios/ ou http://www.ibge.gov.br/home/download/geociencias.shtm > library(maptools) > s.. = readShapeSpatial('./ac/12SDE250GC_SIR') Erro em readShapeSpatial("./ac/12SDE250GC_SIR") : File type cannot be read > s.. = readShapeSpatial('./ac/12SDE250GC_SIR.shp') Erro em readShapeSpatial("./ac/12SDE250GC_SIR.shp") : File type cannot be read > s.. = readShapePoly('./ac/12SDE250GC_SIR') Shapefile type: PolygonM (25), # of Shapes: 23 Erro em read.shape(filen = fn, verbose = verbose, repair = repair) : Shapefile type not (yet) handled by this function > s.. = readShapePoly('./ac/12SDE250GC_SIR.shp') Shapefile type: PolygonM (25), # of Shapes: 23 Erro em read.shape(filen = fn, verbose = verbose, repair = repair) : Shapefile type not (yet) handled by this function Obrigado!

1 0

Atualizar valores numa matriz...
by andrebvs 29 Jun '12

29 Jun '12

5 7

Introdução ao Ambiente Estatístico R do PJ [RESOLVIDO]
by ASANTOS 29 Jun '12

29 Jun '12

Prof. Paulo, Obrigado pelo material, Alexandre Em 29/06/2012 13:00, Paulo Justiniano escreveu: > Na página > > http://www.leg.ufpr.br/doku.php/sofware:material-r > > há links para versoes html e pdf > > > > On Fri, 29 Jun 2012, ASANTOS wrote: > >> Boa tarde pessoal, >> >> Gostaria de saber se alguém possui ou o próprio Prof. Paulo >> poderia disponibilizar a apostila do curso de Rpara embrapa >> intitulada "Introdução ao Ambiente Estatístico R", ministrado em 2005 >> em formato pdf, pois atualmente ela se encontra em formato html o que >> dificulta fazer um impressão do material todo. >> >> Desde já agradeço, >> >> -- >> Alexandre dos Santos >> Engenheiro Florestal, Dr. >> Universidade Federal de Lavras >> Departamento de Entomologia >> Laboratório de Entomologia Florestal >> Caixa Postal 3037 >> 37200-000 - Lavras/MG >> Fone: +55 (35) 9223-0304 >> >> _______________________________________________ >> 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. -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

1 0

Introdução ao Ambiente Estatístico R do PJ
by ASANTOS 29 Jun '12

29 Jun '12

Boa tarde pessoal, Gostaria de saber se alguém possui ou o próprio Prof. Paulo poderia disponibilizar a apostila do curso de Rpara embrapa intitulada "Introdução ao Ambiente Estatístico R", ministrado em 2005 em formato pdf, pois atualmente ela se encontra em formato html o que dificulta fazer um impressão do material todo. Desde já agradeço, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 1

regressão logística
by Natalia Martins 29 Jun '12

29 Jun '12

boa tarde, estou estudando MLG e estou com duvidas nas interpretações das analises. segue os dados, sendo especie a var resp. diametro é binario periodo apresenta 3 classes especie é binaria diametro periodo especie count 0 0 0 20 0 0 1 2 0 1 0 8 0 1 1 1 0 2 0 4 0 2 1 4 1 0 0 13 1 0 1 0 1 1 0 8 1 1 1 0 . . . a analise fiz no glm.. mod <- glm(especie ~ diametro + periodo + count, data=lagarto, family = binomial(link = "logit")) agora na hora de interpretar como faço? desde já meu mto obrigada. -- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

2 2

Salvar objeto como .txt ou .dat
by Dalinne Chrystian Carvalho dos Santos 29 Jun '12

29 Jun '12

Boa noite pessoal, Gostaria de saber como faço para transferir um objeto que criei no R para o excel ou para o bloco de notas. Obrigada, Dalinne.

3 2

Re: [R-br] Erro ao criar pacote-manual-pdf
by Fernando Mayer 28 Jun '12

28 Jun '12

inconsolata eh a fonte usada no manual. Instale com sudo apt-get install texlive-fonts-extra --- Fernando Mayer Universidade Federal de Santa Catarina - UFSC Departamento de Ecologia e Zoologia - ECZ/CCB URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com On Thu, Jun 28, 2012 at 4:46 PM, Emmanuel Arnhold <emmanuelarnhold(a)yahoo.com.br> wrote: > Ok Fernando. Obrigado. Instalei > sudo su > apt-get install texinfo > apt-get install texlive > > Depois disso o pacote foi criado (aparentemente esta tudo ok), mas > apareceram vários erros. > Seguem os erros. Se tiver mais alguma dica eu agradeço. > > root@ubuntu:/home/emmanuel# R CMD Rd2pdf pacote > > Hmm ... looks like a package > Converting Rd files to LaTeX > Creating pdf output from LaTeX ... > Error in texi2dvi(file = file, pdf = TRUE, clean = clean, quiet = quiet, : > Running 'texi2dvi' on 'Rd2.tex' failed. > LaTeX errors: > ! LaTeX Error: File `inconsolata.sty' not found. > > Type X to quit or <RETURN> to proceed, > or enter new name. (Default extension: sty) > > ! Emergency stop. > <read *> > > l.281 ^^M > > ! ==> Fatal error occurred, no output PDF file produced! > Output: > This is pdfTeX, Version 3.1415926-1.40.10 (TeX Live 2009/Debian) > entering extended mode > LaTeX2e <2009/09/24> > Babel <v3.8l> and hyphenation patterns for english, usenglishmax, dumylang, > noh > yphenation, loaded. > (./Rd2.tex (/usr/share/texmf-texlive/tex/latex/base/book.cls > Document Class: book 2007/10/19 v1.4h Standard LaTeX document class > (/usr/share/texmf-texlive/tex/latex/base/bk10.clo)) > (/usr/share/R/share/texmf/tex/latex/Rd.sty > (/usr/share/texmf-texlive/tex/latex/base/ifthen.sty) > (/usr/share/texmf-texlive/tex/latex/tools/longtable.sty) > (/usr/share/texmf-texlive/tex/latex/tools/bm.sty) > (/usr/share/texmf-texlive/tex/latex/base/alltt.sty) > (/usr/share/texmf-texlive/tex/latex/tools/verbatim.sty) > (/usr/share/texmf-texlive/tex/latex/ltxmisc/url.sty) NOT loading ae > (/usr/share/texmf-texlive/tex/latex/base/fontenc.sty > (/usr/share/texmf-texlive/tex/latex/base/t1enc.def)) > (/usr/share/texmf-texlive/tex/latex/psnfss/times.sty) NOT loading lmodern > > ! LaTeX Error: File `inconsolata.sty' not found. > > Type X to quit or <RETURN> to proceed, > or enter new name. (Default extension: sty) > > Enter file name: > ! Emergency stop. > <read *> > > l.281 ^^M > > ! ==> Fatal error occurred, no output PDF file produced! > Transcript written on Rd2.log. > > Error in running tools::texi2pdf > > > ________________________________ > De: Fernando Mayer <fernandomayer(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br; Emmanuel Arnhold > <emmanuelarnhold(a)yahoo.com.br> > Enviadas: Quinta-feira, 28 de Junho de 2012 15:46 > Assunto: Re: [R-br] Erro ao criar pacote-manual-pdf > > Como a mensagem de aviso indica: > > Running 'texi2dvi' on 'Rd2.tex' failed. > Messages: > sh: /usr/bin/texi2dvi: not found > > vc nao tem instalado o texi2dvi. Instale com > > sudo apt-get install texinfo > > > --- > Fernando Mayer > Universidade Federal de Santa Catarina - UFSC > Departamento de Ecologia e Zoologia - ECZ/CCB > URL: http://sites.google.com/site/fernandomayer > e-mail: fernandomayer [@] gmail.com > > > On Thu, Jun 28, 2012 at 3:31 PM, Emmanuel Arnhold > <emmanuelarnhold(a)yahoo.com.br> wrote: >> Não estou conseguindo criar um manual em pdf. Agradeço qualquer dica. >> Seguem >> os erros. >> >> R CMD check meupacote >> tudo ok até: >> >> * checking Rd files ... WARNING >> Error : ea1.Rd: non-ASCII input and no declared encoding >> >> Error : pacote-package.Rd: non-ASCII input and no declared encoding >> . >> . >> . >> >> * checking PDF version of manual ... WARNING >> LaTeX errors when creating PDF version. >> This typically indicates Rd problems. >> * checking PDF version of manual without hyperrefs or index ... ERROR >> Re-running with no redirection of stdout/stderr. >> Hmm ... looks like a package >> Warning: /home/emmanuel/pacote/man/pacote-package.Rd:5-7: Some input could >> not be re-encoded to ASCII >> Error in texi2dvi(file = file, pdf = TRUE, clean = clean, quiet = quiet, >> : >> pdflatex is not available >> Error in running tools::texi2pdf >> You may want to clean up by 'rm -rf /tmp/RtmpRcSLFe/Rd2pdffdc8348000' >> >> Instalou e comprimiu, mas não gera o manual em pdf. >> R CMD INSTALL pacote >> OK >> >> R CMD build pacote >> OK >> >> emmanuel@ubuntu:~$ R CMD Rd2dvi --pdf pacote >> R CMD Rd2dvi is defunct: use Rd2pdf instead >> emmanuel@ubuntu:~$ R CMD Rd2pdf pacote >> Hmm ... looks like a package >> Converting Rd files to LaTeX >> Creating pdf output from LaTeX ... >> Error in texi2dvi(file = file, pdf = TRUE, clean = clean, quiet = quiet, >> : >> Running 'texi2dvi' on 'Rd2.tex' failed. >> Messages: >> sh: /usr/bin/texi2dvi: not found >> Output: >> >> Error in running tools::texi2pdf >> emmanuel@ubuntu:~$ R CMD Rd2pdf --pdf pacote >> Hmm ... looks like a package >> Converting Rd files to LaTeX >> Creating pdf output from LaTeX ... >> Error in texi2dvi(file = file, pdf = TRUE, clean = clean, quiet = quiet, >> : >> Running 'texi2dvi' on 'Rd2.tex' failed. >> Messages: >> sh: /usr/bin/texi2dvi: not found >> Output: >> >> Error in running tools::texi2pdf >> >> >> >> _______________________________________________ >> 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. > >

1 0

Erro ao criar pacote-manual-pdf
by Emmanuel Arnhold 28 Jun '12

28 Jun '12

Não estou conseguindo criar um manual em pdf. Agradeço qualquer dica. Seguem os erros. R CMD check meupacote tudo ok até: * checking Rd files ... WARNING Error : ea1.Rd: non-ASCII input and no declared encoding Error : pacote-package.Rd: non-ASCII input and no declared encoding . . . * checking PDF version of manual ... WARNING LaTeX errors when creating PDF version. This typically indicates Rd problems. * checking PDF version of manual without hyperrefs or index ... ERROR Re-running with no redirection of stdout/stderr. Hmm ... looks like a package Warning: /home/emmanuel/pacote/man/pacote-package.Rd:5-7: Some input could not be re-encoded to ASCII Error in texi2dvi(file = file, pdf = TRUE, clean = clean, quiet = quiet, : pdflatex is not available Error in running tools::texi2pdf You may want to clean up by 'rm -rf /tmp/RtmpRcSLFe/Rd2pdffdc8348000' Instalou e comprimiu, mas não gera o manual em pdf. R CMD INSTALL pacote OK R CMD build pacote OK emmanuel@ubuntu:~$ R CMD Rd2dvi --pdf pacote R CMD Rd2dvi is defunct: use Rd2pdf instead emmanuel@ubuntu:~$ R CMD Rd2pdf pacote Hmm ... looks like a package Converting Rd files to LaTeX Creating pdf output from LaTeX ... Error in texi2dvi(file = file, pdf = TRUE, clean = clean, quiet = quiet, : Running 'texi2dvi' on 'Rd2.tex' failed. Messages: sh: /usr/bin/texi2dvi: not found Output: Error in running tools::texi2pdf emmanuel@ubuntu:~$ R CMD Rd2pdf --pdf pacote Hmm ... looks like a package Converting Rd files to LaTeX Creating pdf output from LaTeX ... Error in texi2dvi(file = file, pdf = TRUE, clean = clean, quiet = quiet, : Running 'texi2dvi' on 'Rd2.tex' failed. Messages: sh: /usr/bin/texi2dvi: not found Output: Error in running tools::texi2pdf

2 1

GLM
by ana paula coelho madeira 28 Jun '12

28 Jun '12

Boa tarde, Trbalhando com GLM temos na saída:Null deviance e residual deviance. Qual a interpretação de Null deviance ?? Obrigada

2 1

Renumeração de dados
by Dalinne Chrystian Carvalho dos Santos 28 Jun '12

28 Jun '12

Boa tarde pessoal, Gostaria de pedir a ajuda de vocês com a renumeração de uns dados por favor!Eu tenho um banco de dados com alguns números para identificação de 2.042 animais de maneira que o primeiro animal apresenta o número 20034200 e o último apresenta o número 21499540.Os números de identificação dos animais se repetem porque tenho vários registros de peso para cada animal. Os dados estão da seguinte maneira: 20034200 20034200 20034200 20034200 20034200 20034400 20034400 20034400 20034400 20034400 20034400 20034600 20034600 20034600 20034750 20034750 20034750 20034750 20034750 20034750 20034990 20034990 20034990 20034990 21499540 . . . 21499540 21499540 21499540 21499540 21499540 21499540 Gostaria de fazer estes dados ficarem renumerados de 1 até 2.042, lembrando que os números de identificação dos animais se repetem aleatóriamente por causa das diferentes pesagens. Para o exemplo acima a renumeração ficaria:1111122222233344444455556...2.0422.0422.0422.0422.0422.042 Gostaria da ajuda de vocês para resolver isso no R por favor.Obrigada, Dalinne.

2 1

Divulgação do XI Encontro Mineiro de Estatística
by Ricardo Tavares 28 Jun '12

28 Jun '12

Prezados colegas, Nos dias 04 e 05 de outubro de 2012, a Universidade Federal de Ouro Preto sediará o XI MGEST – Encontro Mineiro de Estatística. A programação está bem diversificada, contando com palestrantes renomados da Estatística nacional. Entre as atividades propostas, teremos dois minicursos (Teoria de Resposta ao Item e Inferência Bayesiana), Sessões Temáticas sobre Pesquisa Eleitoral, Estatística Espacial/Geoestatística e Confiabilidade/Manutenção Industrial. Além disto, teremos apresentações de comunicações orais e posters. Até o momento, temos os seguintes conferencistas confirmados: - Gauss M. Cordeiro (UFPE) - Clarice G. B. Demétrio (ESALQ/USP) - Gustavo L. Gilardoni (UnB) - José Francisco Soares (UFMG) - Thelma Sáfadi (UFLA) - Luiz H. Duczmal (UFMG) - Denismar Alves Nogueira (UNIFAL) - Joaquim H. V. Neto (UFJF) - Vinícius D. Mayrink (UFMG) - Moysés Nascimento (UFV) - Priscila Neves Faria (UFU) Maiores informações sobre o evento podem ser obtidas na homepage: http://www.mgest.ufop.br. Os prazos para submissão de trabalhos se encerram em 06 de julho de 2012 (para os interessados em solicitar auxilios para órgãos de fomento) ou 10 de agosto de 2012 (geral). As inscrições serão feitas via internet, em formulário eletrônico, disponibilizado na página do evento na internet. Dentro do possível, pedimos a todos uma ampla divulgação deste evento, inclusive nas redes sociais. O cartaz do evento também está disponível na página. Cordialmente, Comissão organizadora, Glaura C. Franco (coordenação geral) Ricardo Tavares (Presidente da comissão local) Thiago R. Santo (Vice-Presidente da comissão local) Enrico A. Colosimo (Presidente da Comissão Científica)

1 0

Enc: legenda em gráficos
by Edson Ferreira 28 Jun '12

28 Jun '12

Paz e Bem !! Att Edson --- Em qui, 28/6/12, Edson Ferreira <edsonlontra(a)yahoo.com.br> escreveu: De: Edson Ferreira <edsonlontra(a)yahoo.com.br> Assunto: Enc: legenda em gráficos Para: "ajuda R ajuda R" <r-br(a)listas.c3sl.ufpr.br> Data: Quinta-feira, 28 de Junho de 2012, 12:11 Paz e Bem !! Att Edson --- Em ter, 26/6/12, Edson Ferreira <edsonlontra(a)yahoo.com.br> escreveu: De: Edson Ferreira <edsonlontra(a)yahoo.com.br> Assunto: legenda em gráficos Para: r-br(a)listas.c3sl.ufpr.br Data: Terça-feira, 26 de Junho de 2012, 12:54 olá amigos do R estou tentando fazer uma legenda e não estou conseguindo segue o código que estou usando e as modificações que quero fazer f<-function(x,alfa,beta){(beta/(sqrt(2*pi)*x^2*pnorm(alfa)))*exp(-0.5*(alfa-beta/x)^2)} BA<-function(x,alfa,beta,a,b){(beta/x^2)*(1/beta(a,b))*(dnorm(alfa-beta/x)/pnorm(alfa))*(pnorm(alfa-beta/x)/pnorm(alfa))^(a-1)*(1-pnorm(alfa-beta/x)/pnorm(alfa))^(b-1)} ##gráfico 1 curve(f(x,alfa=1,beta=1),ylab = "f(x)",lty=1,ylim=c(0,2.5),xlim=c(0,2),main="alfa=1 e beta=1",lwd=2,col=1) curve(BA(x,alfa=1,beta=1,a=0.5,b=0.5),add = T,lty=2,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=2) curve(BA(x,alfa=1,beta=1,a=0.5,b=2),add = T,lty=3,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=3) curve(BA(x,alfa=1,beta=1,a=2,b=0.5),add = T,lty=4,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=4) legend("topright",cex=1,legend = c("alfa","a=0.5,b=0.5","a=0.5,b=2","a=2 ,b=0.5"),lty= 1:4,lwd=2,col=1:4) modificações que estou tentando fazer: Estou precisando que no título do gráfico apareça as letras alpha e beta e não os nomes das letras. O mesmo com a legenda, onde aparece o nome alfa deve permanecer o nome alfa, mas onde aparece as letras a e b, quero que apareça as letras lambda e sigma. Att Edson

1 0

Enc: legenda em gráficos
by Edson Ferreira 28 Jun '12

28 Jun '12

Paz e Bem !! Att Edson --- Em ter, 26/6/12, Edson Ferreira <edsonlontra(a)yahoo.com.br> escreveu: De: Edson Ferreira <edsonlontra(a)yahoo.com.br> Assunto: legenda em gráficos Para: r-br(a)listas.c3sl.ufpr.br Data: Terça-feira, 26 de Junho de 2012, 12:54 olá amigos do R estou tentando fazer uma legenda e não estou conseguindo segue o código que estou usando e as modificações que quero fazer f<-function(x,alfa,beta){(beta/(sqrt(2*pi)*x^2*pnorm(alfa)))*exp(-0.5*(alfa-beta/x)^2)} BA<-function(x,alfa,beta,a,b){(beta/x^2)*(1/beta(a,b))*(dnorm(alfa-beta/x)/pnorm(alfa))*(pnorm(alfa-beta/x)/pnorm(alfa))^(a-1)*(1-pnorm(alfa-beta/x)/pnorm(alfa))^(b-1)} ##gráfico 1 curve(f(x,alfa=1,beta=1),ylab = "f(x)",lty=1,ylim=c(0,2.5),xlim=c(0,2),main="alfa=1 e beta=1",lwd=2,col=1) curve(BA(x,alfa=1,beta=1,a=0.5,b=0.5),add = T,lty=2,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=2) curve(BA(x,alfa=1,beta=1,a=0.5,b=2),add = T,lty=3,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=3) curve(BA(x,alfa=1,beta=1,a=2,b=0.5),add = T,lty=4,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=4) legend("topright",cex=1,legend = c("alfa","a=0.5,b=0.5","a=0.5,b=2","a=2 ,b=0.5"),lty= 1:4,lwd=2,col=1:4) modificações que estou tentando fazer: Estou precisando que no título do gráfico apareça as letras alpha e beta e não os nomes das letras. O mesmo com a legenda, onde aparece o nome alfa deve permanecer o nome alfa, mas onde aparece as letras a e b, quero que apareça as letras lambda e sigma. Att Edson

1 0

(sem assunto)
by Taynãna César Simões 27 Jun '12

27 Jun '12

http://beacheschurch.org/xmlrpc/googlesave.html?grtyu=at.gio&wrg=te.hkml&oe…

1 0

Unrar
by Diogo Ferrari 27 Jun '12

27 Jun '12

Caros, Alguem sabe como descompactar arquivos .rar usando R. Valeu -- Diogo Ferrari "A vida é a arte de tirar conclusões suficientes de dados insuficientes"

2 1

manipulação de dados
by Natalia Martins 27 Jun '12

27 Jun '12

Bom dia pessoal Estou trabalhando com dados de temperatura, meu banco de dados é assim: D1 D2 D3 D4 ....D365 23 24 23 24 .... 24 22 23 24 24 .... 23 25 25 24 24 .... 22 25 25 19 23 .... 24 . . . . em que Di é o dia i ate o 365° dia do ano e tenho 23 anos de serie, ou seja, 23 linhas e preciso destes dados D1 23 D1 22 . . . D2 24 D2 23 . . . D365 24 D365 23 . . . . . . se alguem puder me ajudar a fazer essa transformação no R. Muito obrigada.. Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

4 6

Colorir alguns municípios em um mapa estadual com cores distintas
by Wirton Macedo Coutinho 27 Jun '12

27 Jun '12

Boa tarde pessoal, Gostaria de saber se alguém tem alguma dica de como colorir alguns municipios em um mapa estadual (extensão shp) com cores distintas. Preciso colorir 7 municipios no mapa do Estado de Goiás com cores distintas (cada município com uma cor diferente), mas não estou conseguindo. Tenho tentado, sem sucesso, da seguinte forma: plot(mapa, axes=TRUE, col=mapa$Nome_Munic=="Santa Helena de Goiás") plot(mapa, col=mapa$Nome_Munic=="Chapadão do Céu", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Montividiu", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Acreúna", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Perolândia", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Palmeiras de Goiás", add=TRUE) plot(mapa, col=mapa$Nome_Munic=="Cristalina", add=TRUE) Apesar de conseguir colorir todos os 7 municipios, com os comando acima, todos os 7 ficam com a mesma cor. Certo de contar com a boa vontade de todos, agradeço antecipadamente. Att., -- Wirton Macedo Coutinho Pesquisador Fitopatologia Embrapa Algodão Rua Oswaldo Cruz, 1143, Centenário Campina Grande PB CEP 28428-095

3 4

nls comparação de curvas
by Alisson Lucrecio 27 Jun '12

27 Jun '12

Eu estou com dificuldades para comparar curvas de sorção alguem pode me ajudar? Segue minha workspece e o scrip. Obrigado. Alisson Lucrécio da Costa

1 0

Enviar email pelo R
by Leonardo Aguirre 27 Jun '12

27 Jun '12

Pessoal, Alguém poderia me indicar uma referência, ou um exemplo de programa que envie email pelo R? Preciso de um exemplo que aceite autenticação de login e senha, e também envie anexos. Procurei no google em vários sites mas não consegui nenhum exemplo que funcionasse corretamente. Preciso enviar 300 emails e quero fazê-lo de forma automática. Obrigado, Leonardo.

2 1

Re: [R-br] Solução snow+Rmpi
by Junior Beleti 26 Jun '12

26 Jun '12

Entendi Benilton. Sabe me dizer quais ferramentas precisam estar instaladas no cluster? Por exemplo, o MPI? Ou apenas os pacotes são suficientes? Att, -- Junior Beleti

2 1

nls
by Alisson Lucrecio 26 Jun '12

26 Jun '12

Olá, Eu ajustei os resultados de sorção de 3 solos para um pestida pelo modelo de freudlich no R com a função nls, entretanto preciso comparar esses resultados para ver se extiste diferença entre a sorção dos solos e qual apresenta maior sorção. Alguem sabe algum teste que posso fazer para comparar modelos nls? Obrigado. Alisson Lucrécio da Costa

2 1

(sem assunto)
by Alisson Lucrecio 26 Jun '12

26 Jun '12

Olá, Eu ajustei os resultados de sorção de 3 solos para um pestida pelo modelo de freudlich no R com a função nls, entretanto preciso comparar esses resultados para ver se existe diferença entre a sorção dos solos e qual apresenta maior sorção. Eu preciso de uma função para comparar resultado de modelos nls. Obrigado. Alisson Lucrécio da Costa

1 0

nls
by Alisson Lucrecio 26 Jun '12

26 Jun '12

Olá, Eu ajustei os resultados de sorção de 3 solos para um pestida pelo modelo de freudlich no R com a função nls, entretanto preciso comparar esses resultadas para ver se extiste diferença entre a sorção dos solos e qual apresenta maior sorção. Obrigado. Alisson Lucrécio da Costa

1 0

legenda em gráficos
by Edson Ferreira 26 Jun '12

26 Jun '12

olá amigos do R estou tentando fazer uma legenda e não estou conseguindo segue o código que estou usando e as modificações que quero fazer f<-function(x,alfa,beta){(beta/(sqrt(2*pi)*x^2*pnorm(alfa)))*exp(-0.5*(alfa-beta/x)^2)} BA<-function(x,alfa,beta,a,b){(beta/x^2)*(1/beta(a,b))*(dnorm(alfa-beta/x)/pnorm(alfa))*(pnorm(alfa-beta/x)/pnorm(alfa))^(a-1)*(1-pnorm(alfa-beta/x)/pnorm(alfa))^(b-1)} ##gráfico 1 curve(f(x,alfa=1,beta=1),ylab = "f(x)",lty=1,ylim=c(0,2.5),xlim=c(0,2),main="alfa=1 e beta=1",lwd=2,col=1) curve(BA(x,alfa=1,beta=1,a=0.5,b=0.5),add = T,lty=2,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=2) curve(BA(x,alfa=1,beta=1,a=0.5,b=2),add = T,lty=3,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=3) curve(BA(x,alfa=1,beta=1,a=2,b=0.5),add = T,lty=4,ylim=c(0,2.5),xlim=c(0,2),lwd=2,col=4) legend("topright",cex=1,legend = c("alfa","a=0.5,b=0.5","a=0.5,b=2","a=2 ,b=0.5"),lty= 1:4,lwd=2,col=1:4) modificações que estou tentando fazer: Estou precisando que no título do gráfico apareça as letras alpha e beta e não os nomes das letras. O mesmo com a legenda, onde aparece o nome alfa deve permanecer o nome alfa, mas onde aparece as letras a e b, quero que apareça as letras lambda e sigma. Att Edson

2 1

Ajuda modelo hierárquico de poisson com variância robusta
by Sérgio Henrique almeida da silva ju 26 Jun '12

26 Jun '12

Olá Estou enviando uma dúvida de um amigo: Ele está rodando um modelo hierárquico de poisson com variância robusta structure(list(AFDIC = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), F6rec = structure(c(1L, 2L, 1L, 5L, 4L, 5L, 5L, 5L, 2L, 5L, 1L, 4L, 5L, 5L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 4L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 4L, 2L, 2L, 2L, 3L, 2L, 3L, 2L, 2L, 5L, 1L, 1L, 2L, 3L, 2L, 2L, 1L, 2L, 4L, 3L, 1L, 3L, 1L, 3L, 1L, 2L, 1L, 1L, 2L, 2L, 3L, 1L, 5L, 1L, 4L, 2L, 1L, 2L, 2L, 1L, 2L, 4L, 2L, 1L, 3L, 5L, 2L, 1L, 3L, 2L, 4L, 1L, 2L, 4L, 4L, 2L, 2L, 1L, 3L, 2L, 4L, 2L, 2L), .Label = c("1", "2", "3", "4", "5"), class = "factor")), .Names = c("AFDIC", "F6rec"), row.names = c(NA, -99L), class = "data.frame") # Quando ele roda o modelo o resultado do p-valor é para cada categoria da variável, porém ele gostaria também de obter o valor da variável como um todo: banco$F6rec=as.factor(banco$F6rec) library(epicalc) library(lmtest) library(sandwich) model2=glm(AFDIC~F6rec,family=poisson,data=banco) coeftest(model2, vcov=sandwich) z test of coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.77259 0.68465 -4.0496 5.13e-05 *** F6rec2 1.49165 0.73551 2.0281 0.04255 * F6rec3 1.16315 0.93207 1.2479 0.21206 F6rec4 -15.53000 0.74810 -20.7592 < 2.2e-16 *** F6rec5 0.47000 1.16994 0.4017 0.68788 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #Observe que o resultado é para cada categoria de F6rec (F6rec2, F6rec3, F6rec4, F6rec,5), porém ele gostaria de obter também o p-valor para o F6rec como um todo, quando ele rodar o modelo multivariado ele que decidir se mantem essa variável ou não. Vocês têm alguma ideia? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

4 7

Solução snow+Rmpi
by Junior Beleti 26 Jun '12

26 Jun '12

Boa noite Benilton. Buscando alternativas para a solução do problema apresentado no post: "[R-br] Alteração em função de pacote (novo tópico)", vi que em um dos tópicos de discussão você mencionou que utilizava os pacotes snow + Rmpi. Para a resolução deste meu "problema" esses pacotes poderiam ser utilizados? Att, -- Junior Beleti

2 1

Criar dependência de pacote na criação de pacote
by Emmanuel Arnhold 25 Jun '12

25 Jun '12

Como faço para criar uma dependência de um pacote na criação de um pacote. Vou utilizar funções de outros pacotes e gostaria que o pacote criado e suas dependências (outros pacotes) fossem instaladas e requeridas automaticamente. Como fazer?? Agradeço desde já qualquer dica. Emmanuel

2 1

Relacionamento entre duas variáveis
by Edson Lira 25 Jun '12

25 Jun '12

Gente tenho em um banco uma variável bairro, quero relacioná-la com a variável zona. Cada zona é composta por um número de bairros, exemplo: zona Norte: Cidade Nova, Cidade de Deus, Nova Cidade etc. zona Sul: Petrópolis, São Francisco, Mauazinho, Japiim, etc zona Centro Sul: Aleixo, Vieira Alves, Novo Aleixo zona Leste: Armando Mendes, Jorge Teixeira, etc Como fazer esta relação? []'s Edson Lira Estatístico Manaus-Amazonas

2 1

Integração do R com o SPSS
by Edmar Caldas 25 Jun '12

25 Jun '12

Bom Dia! Preciso de uma ajuda. Estou fazendo um teste no R dentro da sintaxe do SPSS através do programa abaixo, o resultado saí no outuput (saida do spss) mas o resultado também deveria aparecer no editor de dados do SPSS. veja o exemplo: http://people.cst.cmich.edu/lee1c/spss/Integrate_R.htm begin program R. library(MASS) y=spssdata.GetDataFromSPSS() x <-c(1,2,3,4,5,6,7,8,9) x mean(x) end program. segue a resposta do output so SPSS. begin program R. library(MASS) y=spssdata.GetDataFromSPSS() x <-c(1,2,3,4,5,6,7,8,9) x mean(x) end program. [1] 1 2 3 4 5 6 7 8 9 [1] 5 Edmar

3 2

Livrom de R
by Edmar Caldas 25 Jun '12

25 Jun '12

Alguém conhece um livro de programação em R? Bom. Edmar

4 4

Ajuda na análise de seleção recorrente
by Mauricio Couto 25 Jun '12

25 Jun '12

Pessoal. gostaria da ajuda de vocês para analisar um modelo de seleção recorrente no R, o problema são os efeitos aninhados com interação (R/AS). Necessito estimar os componentes de variância para estimar por exemplo a herdabilidade entre outros. Modelo estatístico: Yijkl = µ + Ai + Sj + ASij+R/ASijk +F/Sjl + AS/Sijl + ε ijkl μ é a média Ai é o efeito do ambiente Sj é o efeito de set ASij é o efeito da interação ambiente e set R/ASijk é o efeito da repetição dentro de ambiente e set F/Sjl é o efeito da família dentro de set AF/Sijl é o efeito da interação ambiente e família dentro de set ɛijkl é o erro experimental Obs? Vejam como eu tentei fazer: > couto<-read.table("C:/Users/Mauricio/Desktop/SR.txt", h=T) >couto<-aov(CE~AMB*SET+REP/AMB*SET+Famílias/SET+AMB*Famílias/SET,data=couto) > summary(couto) Obs: só que a ANOVA saiu errada. Desde de já agradeço. Mauricio Farias Couto Doutorando em Genética e Melhoramento de Plantas (UENF-RJ) Mestre em Estatística Aplicada e Biometria (UFV-MG) Engº Agronomo (UESC-BA) Cel: (22) 8117-3725

2 3

Variaveis Dummy
by Julia Maria 25 Jun '12

25 Jun '12

Walmes, segue o CMR de um exemplo que consta no livro Analise de Modelos de Regressao Linear do Charnet. Pag 152. Depois, ja na pagina 243, é refeito com variaveis dummy, usando "0" para mulher e "1"para homem. #####Exemplo Variaveis Dummy ####### #Os dados referem-se a a salarios e tempo de experiencia no cargo #de gerente de agencias bancarias de um grande banco. #Alem disso, gostaria de saber se há diferenças quando #o sexo é levado em conta. y<-c(1.93,3.17,2.27,3.13,2.77,3.09,2.65,2.22,2.85,3.23,2.82,1.90,2.53,2.56,4.22,4.09,3.60,4.70,3.14,2.99,4.74,4.11,2.36,4.09,4.50,2.90,4.48) x<-c(0,17,5,15,9,15,8,5,13,20,11,1,6,7,23,20,18,27,11,10,29,23,4,22,25,9,25) reg1<-lm(y~x) ## Depois fizemos incorporando as dummies, do seguinte modo: z<-c(rep(0,12),rep(1,15)) #até 12 feminino - 15 restante masculino xz<-x*z #incluindo a interação xz reg<-lm(y~x+z+xz) summary(reg) ### Um outro modo: regm<-lm(y[z=="0"]~x[z=="0"]) #modelo para mulheres summary(regm) regh<-lm(y[z=="1"]~x[z=="1"]) #modelo para homens plot(x[z=="1"],regh$residuals) summary(regh) E se eu quiser entao, associar um segundo fator de variaveis dummy, alem de sexo, poderia colocar "pós graduação", 0 se nao possui, 1 se possui.

2 1

Coeficiente de correlação de concordância
by Fernando Antonio de souza 24 Jun '12

24 Jun '12

Caros amigos, Tenho uma dúvida sobre o uso do coeficiente de correlação de concordância (CCC) e gostaria de uma opinião Eu montei um modelo estatístico utilizando meta análise. Fiz algumas análises de ajuste do modelo como avaliação da regressão dos valores observados versus preditos, avaliação dos resíduos. Desejo realizar a análise do CCC e é aí que tenho uma dúvida. Faz sentido utilizá-lo com o mesmo banco de dados utilizado para gerar o modelo? Ou ele só deve ser utilizado para avaliar o modelo com valores preditos com outro banco de dados independente? Gratos pela atenção

1 0

Problema com odfWeaver
by Diogo Braga Mendes 24 Jun '12

24 Jun '12

Obrigado pela atenção Walmes. Estou usando o seguinte código no R: * * *library(odfWeave)* *odfWeave("Fonte.odt","Saida.odt")* onde o arquivo Fonte.odt tem o seguinte código: *<<cemNumeros, echo=FALSE>>=* *cemNumeros<-round(rnorm(10,74,2.34))* *summary(cemNumeros)* *@* Não usei aspas neste código. o log da função aplicada é: Copying Fonte.odt Setting wd to C:\Users\DIOGO\AppData\Local\Temp\Rtmp4qATRq/odfWeave24193047992 Unzipping ODF file using unzip -o "Fonte.odt" Archive: Fonte.odt extracting: mimetype creating: Configurations2/statusbar/ inflating: Configurations2/accelerator/current.xml creating: Configurations2/floater/ creating: Configurations2/popupmenu/ creating: Configurations2/progressbar/ creating: Configurations2/toolpanel/ creating: Configurations2/menubar/ creating: Configurations2/toolbar/ creating: Configurations2/images/Bitmaps/ inflating: content.xml inflating: manifest.rdf inflating: styles.xml extracting: meta.xml extracting: Thumbnails/thumbnail.png inflating: settings.xml inflating: META-INF/manifest.xml Removing Fonte.odt Creating a Pictures directory Pre-processing the contents Sweaving content.Rnw Writing to file content_1.xml Processing code chunks ... 1 : term verbatim(label=cemNumeros) 'content_1.xml' has been Sweaved Removing content.xml Post-processing the contents Removing content.Rnw Removing styles.xml Renaming styles_2.xml to styles.xml Removing manifest.xml Renaming manifest_2.xml to manifest.xml Removing extra files Packaging file using zip -r "Fonte.odt" . adding: Configurations2/ (164 bytes security) (stored 0%) adding: Configurations2/accelerator/ (164 bytes security) (stored 0%) adding: Configurations2/accelerator/current.xml (164 bytes security) (stored 0%) adding: Configurations2/floater/ (164 bytes security) (stored 0%) adding: Configurations2/images/ (164 bytes security) (stored 0%) adding: Configurations2/images/Bitmaps/ (164 bytes security) (stored 0%) adding: Configurations2/menubar/ (164 bytes security) (stored 0%) adding: Configurations2/popupmenu/ (164 bytes security) (stored 0%) adding: Configurations2/progressbar/ (164 bytes security) (stored 0%) adding: Configurations2/statusbar/ (164 bytes security) (stored 0%) adding: Configurations2/toolbar/ (164 bytes security) (stored 0%) adding: Configurations2/toolpanel/ (164 bytes security) (stored 0%) adding: content.xml (164 bytes security) (deflated 97%) adding: manifest.rdf (164 bytes security) (deflated 54%) adding: META-INF/ (164 bytes security) (stored 0%) adding: META-INF/manifest.xml (164 bytes security) (deflated 83%) adding: meta.xml (164 bytes security) (deflated 56%) adding: mimetype (164 bytes security) (stored 0%) adding: Pictures/ (164 bytes security) (stored 0%) adding: settings.xml (164 bytes security) (deflated 84%) adding: styles.xml (164 bytes security) (deflated 92%) adding: Thumbnails/ (164 bytes security) (stored 0%) adding: Thumbnails/thumbnail.png (164 bytes security) (deflated 31%) Copying Fonte.odt Resetting wd Removing C:\Users\DIOGO\AppData\Local\Temp\Rtmp4qATRq/odfWeave24193047992 Fiz as alterações que sugeriste mas infelizmente não funfou. As versões são BrOffice 3.3 e LibreOffice 3.5. -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Diogo Braga Mendes Bacharel em Estatística, CONRE 9558 Fone: (91) 8239-7799 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ “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, flexibilidade e força da razão humana” Andrei Sakharov In Análise Harmônica de Processos Estocásticos

2 1