Abril 2016 - R-br - DINF UFPR

Completar linhas de uma variável
by Amikobh 27 Abr '16

27 Abr '16

Prezados, Tenho um banco com mais de 8000 linhas e 27 variáveis. Uma delas é a "Processo". Nela, tenho várias linhas em branco as quais gostaria de completar com números sequenciais, como exemplo, de 1 até o numero total de linhas em branco (1,2,3,4,.....). De tal forma que depois do preenchimento eu tenha um novo banco com estas modificações. Agradecido pela ajuda! PROCESSO 5.155.855-8 5.153.696-8 5.154.220-6 5.155.401-1 5.155.960-6 5.157.157-7 5.157.658-4 Elerson

3 3

Como recriar um componente de PCA obtido com svyprcomp()
by Leonardo Ferreira Fontenelle 27 Abr '16

27 Abr '16

Boa tarde a todos! Estou tentando escrever os resultados de uma análise de componentes principais de forma que os leitores possam calcular o primeiro componente a partir dos dados sem a necessidade de utilizar um comando de análise de componentes principais propriamente dito. Minha ideia é informar a codificação das variáveis, as cargas, os desvios-padrão, e então orientar os leitores a multiplicar os valores das variáveis pelas cargas e dividir pelos desvios-padrão e então somar tudo. Minha expectativa é de que o resultado seja altamente correlacionado àquele obtido por uma análise de componentes principais propriamente dita. No entanto, como explico a seguir, não estou conseguindo, e de alguma forma isso tem a ver com eu estar usando svyprcomp() em vez de prcomp(). Para começar, mostro uma situação em que isso funciona: pca1 <- prcomp(USArrests, scale = TRUE) table1 <- data.frame(loadings = pca1$rotation[, "PC1"], scale = pca1$scale, coef = NA_real_, row.names = row.names(pca1$rotation)) table1$coef <- table1$loadings / table1$scale firstcomponent1 <- with(USArrests, Murder * table1["Murder", "coef"] + Assault * table1["Assault", "coef"] + UrbanPop * table1["UrbanPop", "coef"] + Rape * table1["Rape", "coef"]) cor(firstcomponent1, pca1$x[, "PC1"]) # resultado: 1,00 Agora, uma situação onde isso não funciona: library("survey") data(api) dclus2 <- svydesign(id = ~ dnum + snum, fpc = ~ fpc1 + fpc2, data = apiclus2) pca2 <- svyprcomp(~ api99 + api00 + ell, design = dclus2, scale = TRUE, scores = TRUE) table2 <- data.frame(loadings = pca2$rotation[, "PC1"], scale = pca2$scale, coef = NA_real_, row.names = row.names(pca2$rotation)) table2$coef <- table2$loadings / table2$scale firstcomponent2 <- with(apiclus2, api99 * table2["api99", "coef"] + api00 * table2["api00", "coef"] + ell * table2["ell", "coef"]) cor(firstcomponent2, pca2$x[, "PC1"]) # resultado: 0,506 Como eu consigo cargas ("loadings") e desvios-padrão ("scale") que eu possa usar de forma a conseguir calcular à mão algum vetor altamente correlacionado com o primeiro componente do svyprcomp()? Grato! Leonardo Ferreira Fontenelle[1] Links: 1. http://lattes.cnpq.br/9234772336296638

2 6

Levene test (original)
by Luiz Roberto Martins Pinto 27 Abr '16

27 Abr '16

Desejo utilizar o Levene test (o original, não o modificado). Alguém conhece alguma função? Não estou conseguindo utilizar a biblioteca 'car', que tem esta função. Luiz Roberto. Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia-Brasil luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831

3 5

Estatísticas do R
by Andre Oliveira 26 Abr '16

26 Abr '16

Pessoal,preciso de estatísticas do R tais como nº de bibliotecas oficiais, pacotes que o R já vem na instalação. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

8 8

Vale a pena tratar esses missings ?
by Elias Carvalho 26 Abr '16

26 Abr '16

Pessoal no meu caso tenho variáveis com alto indice de missing (53 A 79%) : [1] " variable: V1 number.na: 1456 percent.na: 0.73%" [1] " variable: V2 number.na: 1052 percent.na: 0.53%" [1] " variable: V3 number.na: 1577 percent.na: 0.79%" [1] " variable: V4 number.na: 1564 percent.na: 0.78%" Existe alguma regra ? Vale a pena aplicar um método pra tratar esses missing quando a quantidade é alta ?

2 1

Erro no integrate
by Daniela Recchia 26 Abr '16

26 Abr '16

Olá pessoal, quero comparar retas desenhadas por crianças (a mão, elas viam o gráfico com pontos e tinham que desenhar uma reta que melhor se ajustasse aos pontos) com 3 métodos diferentes que resultam na "reta correta". Para isso ploto as retas sempre de duas a duas (criança com método 1, criança com método 2 e criança com método 3) e calculo a área entre as duas retas (integral). Quero no fim poder olhar para os 3 valores de cada criança e achar o mínimo (melhor método). Para isso (tentei) criei um código que infelizmente não funciona... Os valores dos coeficientes das retas dos 3 métodos de referência são fixos (b e s), tenho uma matrix com os 75 valores dos coeficientes para as crianças, os coeficientes são m e t: Quero armazenar cada área clauclada em uma linha e na coluna específica do método de referênciaarea<-function(m,t,na.rm=TRUE) { ‪#‎na‬.rm não funciona u <- matrix(nrow=length(m),ncol=3) colnames(u)<-c("Method1","Method2","Method3") for(i in 1:length(m)){ for (j in 1:3){ f1<-function(x, m, t) {m*x+t} # m e t variáveis locais mt<-function(x) {b[j]*x+s[j]} # b e s variáveis globais dif<-function(x, m,t) {abs(f1(x,m,t)-mt(x))} int[i]<-integrate(dif, -1,1,m = m[i], t = t[i])$value u[i,j]<-int[i] }} return(u)} Recebo o erro de que o objeto "int" não foi encontrado Alguém teria uma sugestão? -- Daniela Recchia “The standard of success in life is absolutely the amount of joy you feel” Esther & Jerry Hick

1 0

Nuvem de Tags
by Elias Carvalho 25 Abr '16

25 Abr '16

Bom dia pessoal Tenho um arquivo formato .csv com 1800 publicações sobre redes bayesianas eu tenho uma coluna identificando o tipo de rede bayesiana, como por exemplo: Baysian Networks (BN), Dynamic Bayesian Networks (DBN), Gaussian Bayesian Networks (GBN), etc... Eu gostaria de ler essa coluna e montar uma nuvem de tags para representar graficamente a frequencia com que esses termos aparecem. Alguem pode dar uma dica ? Obrigado -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* _____________________________________________ Prof. Elias César Araújo de Carvalho CV: http://lattes.cnpq.br/4248328961021251

3 5

Re: [R-br] Razão de prevalência
by sznelwar 25 Abr '16

25 Abr '16

3 2

Otimização - Fatorial
by Andre Oliveira 23 Abr '16

23 Abr '16

O meu e-mail foi respondido, porém, não sei o motivo de não ter recebido de volta a respostas em meu e-mail! Estou reenviando com a respostas e minha outra colocação. E-mail originalPrezad@s;estou iniciando um projeto em química onde será montado umexperimento no esquema fatorial onde tenho alguns fatores definidos a princípioe gostaria de otimizar este processo. A questão é que tenho 4respostas que serão medidas para cada tratamento. Gostaria desaber se a biblioteca rsm, lm e aov seriam a mais indicada para que eu possa analisar depois de montar o experimento e buscar o deslocamentosnecessários a fim de obter a otimização do processo nospróximos experimentos. Estou buscando material nesta linha e usando o R de preferência. Minha dúvida também se refere ao fato que tenho encontrado materiais com modelos com apenas uma resposta. Sealguém puder ajudar ficarei grato. O Ivan Bezerra Allaman respondeu da seguinte forma O é que eu entendi sua dúvida... vars=list()for(i in 1:30) vars[[i]] <- rnorm(48)names(vars) <- paste('var',1:30,sep='')vars varis <- do.call('cbind',vars) da <- expand.grid(FatA=rep(LETTERS[1:6],4), FatB=letters[1:2]) dad <- data.frame(da,varis) tabs <- apply(dad[,-c(1:2)], 2, function(x) aov(x ~ FatA*FatB,dad))lapply(tabs,summary) Oi Ivan,no caso o que procuro é uma função que faça o que você escreveu mesmo, mas no lugar de x e aov seria algo multivariado como abaixo.fit=manova(varis~da, data= dad) Neste caso a manova seria pouco informativo sobre o otimização do processo. Queria algum modelo com a estrutura do fit que escrevi, mas com mais informações sobre a otimização individual para cada y os fatores em questão. obrigado André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

1 0

[off] Ilustração de equações em uma lousa
by Listeiro 037 23 Abr '16

23 Abr '16

Olá a todos. Tenho uma dúvida que não sei onde posso esclarecer-me. Coloquei uma imagem no serviço do imgur para que possam ver. É só clicar no link. http://i.imgur.com/gBIwI1l.jpg A dúvida é se na imagem da lousa acima as equações fazem algum sentido, se no mínimo lembram alguma coisa com algum sentido, matematicamente falando. Se seriam, por exemplo, algo de Estatística ou Análise. Como temos bastante estatísticos aqui na lista creio encontrar a solução. Desde já agradeço.

1 0