- 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