- R-br - DINF UFPR

Índice fora de limites
by Thiago de Andrade 01 Fev '13

01 Fev '13

Amigos, estou tentando rodar uns dados de um amigo com o pacote ltm, (Teoria da Resposta ao Item). Estou recebendo uma mensagem de erro "índice fora de limites". Algum de vocês pode me dar alguma luz acerca da localização deste problema? Obrigado. -- Thiago de Andrade Mestrando em Psicologia Social, PPGP-UFPB Núcleo de Pesquisas em Desenvolvimento Sócio-Moral (NPDSM) (83) 9642 4099

3 8

Gráfico dupla entrada
by Stephan M. CARVALHO 01 Fev '13

01 Fev '13

Bom dia pessoal, Preciso plotar dados em um gráfico de dupla entrada, como ocorre por exemplo com dados climatológicos (temperatura/precipitação ~ tempo). Por acaso alguém conhece alguma forma de fazer? Obrigado, Stephan -- *************************************************** STEPHAN M. CARVALHO "Honeybee protection and ecotoxicology" Universidade Estadual Paulista Centro de Estudo de Insetos Sociais Av. 24A, 1515 Bela Vista 13.506-900 Rio Claro/SP Brazil smalfitano(a)uol.com.br +55 19 3526-4284 (university) +55 19 3534-8523 (fax) +55 35 9103-9822 (mobile) +55 35 9949-0707 (mobile) Antes de imprimir, pense no ambiente! Before printing, think in the environment! Avant d'imprimer, pensez à l'environnement! ***************************************************

4 4

Informação: R como software livre...
by andrebvs＠bol.com.br 01 Fev '13

01 Fev '13

3 3

Nayara Santos is out of the office.
by Nayara Santos 31 Jan '13

31 Jan '13

I will be out of the office starting 30/01/2013 and will not return until 14/02/2013. Gentileza entrar em contato com: Flávia Novaes para testes. Marcelo Prieto e Vicente Venancio para definições, dúvidas e aprovações.

1 0

Gráfico no ggplot com escala x variável
by Alisson Lucrecio 31 Jan '13

31 Jan '13

Caros colegas, Boa noite. Classifiquei um conjunto de dados de poros por diâmetro em diferentes classes levando em consideração o numero de poros em cada classe. Os resultados me geraram outro conjunto de dados com o numero de poros em cada classe, sendo 30 classes de 0.1 a 1 mm e 1 classes de 1 a 5 mm. Preciso plotar um gráfico no ggplot2 com o diâmetro no x e numero de poros em cada classe no y, mas quando plot gera um gráfico da figura em anexo. Não consigo ter uma boa visualização dos resultado nesse plot. Então preciso fazer um corte no eixo x ou tentar uma escala x automática de 0.1 a 1 mm e escala x fixa de 1 a 5 mm (fixo no x igual a 3 mm). Alguém sabe como fazer isso? Obrigado. Alisson Lucrécio da Costa

2 2

Ajuda com o pacote randomLCA
by Victor Eduardo 30 Jan '13

30 Jan '13

Pessoal, alguém poderia me ajudar def.lcarandom <- randomLCA(definicao[,1:4], freq=definicao$x, initmodel=def.lca, random=TRUE, blocksize=1, quadpoints=190, probit=TRUE) > summary(def.lcarandom) Classes AIC BIC logLik penlogLik Link 2 3842.735 3896.336 -1911.368 -1911.373 Probit Class probabilities Class 1 Class 2 0.8213 0.1787 Conditional outcome probabilities Group.1 Group.2 Group.3 Group.4 Class 1 0.0000 0.0065 0 0.0011 Class 2 0.9612 1.0000 1 1.0000 Marginal Outcome Probabilities Group.1 Group.2 Group.3 Group.4 Class 1 0.0000 0.1987 0 0.1479 Class 2 0.7253 0.9545 1 1.0000 Loadings 2.76629> def.lcarandom$outcomep [,1] [,2] [,3] [,4] [1,] 5.653499e-50 0.006483299 5.024727e-33 0.001059331 [2,] 9.612131e-01 0.999999659 1.000000e+00 1.000000000 A primeira parte postada acima funciona normalmente, sendo que quando vou utilizar a função outcome.probs ela não funciona, ela fica "rodando" o tempo inteiro e não sai mais disso. def.lcarandom.probs<-outcome.probs(def.lcarandom, boot=T, R=100) Abraços!

1 0

geoestatistica
by Natalia Martins 30 Jan '13

30 Jan '13

Prezados, estou tentando comparar a interpolação por krigagem e pelo inverso da distancia, sempre usei somente a krigagem e agora fui questionada e resolvi procurar e e estudar mais métodos..no entanto no exemplo da cran (a seguir), nao consegui entender de onde vem os dados "meuse.grid", uma vez que estes contem as distancias. Alguem poderia me ajudar? # Inverse distance interpolation with inverse distance power set to .5: # (kriging variants need a variogram model to be specified) data(meuse) data(meuse.grid) meuse.gstat <- gstat(id = "zinc", formula = zinc ~ 1, locations = ~ x + y, data = meuse, nmax = 7, set = list(idp = .5)) meuse.gstat z <- predict(meuse.gstat, meuse.grid) library(lattice) # for levelplot levelplot(zinc.pred~x+y, z, aspect = mapasp(z)) -- Att, 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

3 2

rodar função para cada nível de um fator
by marcelo costa 30 Jan '13

30 Jan '13

Prezados membros, Peço a ajuda de vocês para o seguinte problema: Eu tenho uma planilha com dois fatores (“espécie” e “sexo”) e uma variável contínua (“tamanho”). Quero analisar se existe diferença no tamanho de machos e fêmeas para cada espécie. Assim, preciso fazer um teste T para cada espécie. O problema é que eu tenho 288 espécies (uma planilha de 6638 linhas) e, portanto, fazer um subset para cada espécie não me parece uma solução muito boa (e provavelmente não é!). Eu queria fazer um teste T para cada nível do fator "especie".Alguém tem alguma ideia de como eu poderia fazer isso sem precisar fazer um subset para cada espécie? A princípio eu achei que isso seria fácil (e provavelmente seja), no entanto o meu pouco conhecimento sobre o R não está me ajudando. Eu estou enviando um exemplo hipotético do formato da minha planilha: especie<-rep(c("a","b","c"),each=6) sexo<-rep(c("f","m"), each = 3, len = 18) tamanho<-c(1.5,1.6,1.8,1.3,1.1,1.2,2.5,2.7,2.6,2.1,2.2,2.1,3.5,3.6,3.7,3.1,3.3,3.2) data<-data.frame(especie,sexo,tamanho) Um abraço a todos e obrigado pela atenção. Marcelo Costa

5 6

Summary com todos os tipos da variável
by Fátima Lima Paula 28 Jan '13

28 Jan '13

Prezados, tenho uma variável (diagnóstico principal) que apresenta muitos tipos diferentes. Quando peço summary dessa variável, retorna dessa forma: summary(reinter$diagprinc) S720 S721 S722 T845 J189 T931 A419 N390 I500 I64 493 300 77 40 34 29 28 27 21 20 T840 J159 J158 M861 R02 S730 T814 I792 A499 I802 20 19 16 14 14 14 14 13 12 12 M966 R262 A408 L89 T813 T818 J188 J960 M840 T981 12 12 11 11 11 11 10 10 10 10 A46 I739 L989 M009 S723 T018 Y881 A418 A480 E46 9 9 9 9 9 9 9 7 7 7 E86 S525 A488 E43 G459 I509 I702 K921 K928 M171 7 7 6 6 6 6 6 6 6 6 M868 T841 T847 T857 Z470 A490 E145 I269 K566 R69 6 6 6 6 6 5 5 5 5 5 T848 Z988 A415 A498 E878 G319 I10 I743 J180 J449 5 5 4 4 4 4 4 4 4 4 K920 L97 M841 M869 N179 N189 S818 T888 E148 G20 4 4 4 4 4 4 4 4 3 3 I110 I200 I608 I635 I809 J80 K818 L031 L038 M169 3 3 3 3 3 3 3 3 3 3 M842 M848 M860 N111 N180 R33 S065 S423 S424 (Other) 3 3 3 3 3 3 3 3 3 218 Como faço para que apareçam todas? Obrigada Fátima

2 2

Re: [R-br] Método de Agrupamento Ward (bruno vidigal)
by Marcos Vital 28 Jan '13

28 Jan '13

Bruno, na verdade o método de Ward precisa sim da matriz de distância como os outros métodos de clsuter; o que muda é a forma como ele cria o dendograma. Ele calcula a soma dos quadrados das distâncias entre os objetos, dividida pelo número de objetos, buscando minimizá-la. Então é necessária uma matriz de distância como nos outros casos, e o que muda na função hclust é o argumento method="ward". Se tiver acesso, veja o excelente *Numerical Ecology with R* (Bocard, Gillet & Legendre, 2011), que é uma excelente fonte para estudar métodos multivariados no R (o foco é em Ecologia, mas com certeza pode ser útil em outros casos). Abraços Marcos -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

3 2