Janeiro 2012 - R-br

ajuda em teste para comparação de frequências
by Flávia Machado 15 Jan '12

15 Jan '12

Galera, preciso muito da ajuda de vocês, tenho umas análises pra fazer mas não sei quais testes aplicar para responder as minhas perguntas. Tenho um conjunto de dados que correspondem a frequência com que indivíduos de uma população realizam determinados comportamentos ao longo do período de um ano. Quero verificar se há diferença significativa nas frequências dos comportamentos entre os indivíduos neste período. Além disso, se as frequências dos comportamentos de um determinado indivíduo diferem da média dos comportamentos dos demais indivíduos. Ou seja, quero ver: os indivíduos possuem mesmo comportamento? Se não, como é essa diferença e como ela varia ao longo do ano? Essa diferença aumenta, diminui ou se mantêm estável ao longo do ano? Os comportamentos do indivíduo selecionado difere dos demais? Caso haja diferença, verificar o quão diferente e se essa diferença nas freq. diminuem, aumentam, ou se mantêm constantes ao longo do ano. Imagino que para saber se existe diferença nas frequências eu poderia fazer um qui-quadrado, mas não sei fazer. Agora para o outro questionamento nem faço ideia de que teste aplicar. Segue o CMR de parte dos meus dados para vocês verem como eles estão. OBS: as letras antes do "." na função se referem ao indivíduos e as depois dele aos comportamentos. Rep - Repouso; Ali - Alimentação; D - deslocamento; S - Social; O - outros comportamentos. Cada linha corresponde a um dia amostral (71 no total). São seis indivíduos e o que eu quero comparar com os demais está nomeado "RE". # Essas são as frequências dos comportamentos de cada indivíduo: Rep.RE<-c(0.000,0.714,0.500,1.000,0.500,0.636,1.000,0.700,0.600,0.000,0.789,0.417,0.000,0.941,0.200,0.500,0.462,1.000,0.500,0.500,0.273,0.643,0.467,0.353,0.000,0.364,0.250,0.316,0.273,0.429,0.000,0.000,0.258,0.182,0.421,0.667,0.250,0.688,0.167,0.400,0.400,0.692,0.500,0.333,0.400,0.091,0.360,0.333,0.882,0.583,0.400,0.724,0.000,0.636,0.522,0.563,0.368,0.667,0.633,0.640,0.500,0.571,0.621,0.375,0.522,0.000,0.500,0.800,0.375,0.000,0.375) Ali.RE<-c(0.000,0.000,0.500,0.000,0.000,0.182,0.000,0.300,0.200,1.000,0.105,0.250,0.000,0.000,0.600,0.188,0.462,0.000,0.000,0.321,0.091,0.143,0.200,0.353,0.400,0.364,0.450,0.579,0.727,0.357,0.000,0.800,0.484,0.727,0.316,0.333,0.500,0.188,0.556,0.450,0.500,0.154,0.000,0.333,0.400,0.818,0.480,0.444,0.000,0.333,0.500,0.172,0.333,0.242,0.348,0.188,0.421,0.222,0.000,0.200,0.100,0.143,0.172,0.375,0.304,0.000,0.167,0.000,0.500,0.000,0.375) D.RE <-c(1.000,0.286,0.000,0.000,0.500,0.182,0.000,0.000,0.200,0.000,0.105,0.250,1.000,0.059,0.200,0.313,0.077,0.000,0.250,0.179,0.545,0.143,0.333,0.235,0.600,0.136,0.200,0.105,0.000,0.214,1.000,0.200,0.226,0.091,0.211,0.000,0.250,0.125,0.222,0.150,0.100,0.077,0.500,0.333,0.200,0.091,0.160,0.222,0.118,0.083,0.100,0.103,0.667,0.121,0.130,0.188,0.211,0.111,0.100,0.160,0.300,0.286,0.172,0.250,0.174,0.500,0.333,0.200,0.125,0.000,0.125) S.RE <-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.083,0.000,0.000,0.000,0.000,0.000,0.000,0.250,0.000,0.091,0.071,0.000,0.059,0.000,0.136,0.100,0.000,0.000,0.000,0.000,0.000,0.032,0.000,0.000,0.000,0.000,0.000,0.056,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.063,0.000,0.000,0.267,0.000,0.100,0.000,0.034,0.000,0.000,0.500,0.000,0.000,0.000,0.000,0.125) O.RE <-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.053,0.000,0.000,0.000,0.000,0.000,0.000,0.077,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000) Rep.SO<-c(0.667,0.833,0.750,0.500,0.000,0.789,0.889,0.500,0.250,0.625,0.000,0.733,0.727,0.500,0.563,0.429,0.455,0.750,0.500,0.600,0.857,0.000,0.500,0.550,0.444,0.308,0.200,0.500,0.300,0.200,0.333,0.000,0.235,0.000,0.278,0.154,0.429,0.333,0.091,0.393,0.619,0.692,0.522,0.500,0.000,0.333,0.154,0.429,0.824,0.619,0.143,0.385,0.000,0.368,0.182,0.714,0.471,0.464,0.400,0.500,0.000,0.545,0.560,0.643,0.167,0.200,0.833,0.000,0.615,0.250,0.000) Ali.SO<-c(0.000,0.083,0.000,0.357,0.000,0.158,0.000,0.429,0.500,0.125,0.000,0.067,0.091,0.375,0.250,0.429,0.091,0.000,0.500,0.200,0.000,0.333,0.250,0.200,0.222,0.308,0.600,0.250,0.550,0.467,0.333,0.889,0.706,0.000,0.333,0.308,0.571,0.400,0.364,0.393,0.286,0.231,0.217,0.000,0.333,0.000,0.538,0.214,0.059,0.143,0.429,0.385,0.000,0.421,0.727,0.000,0.353,0.321,0.000,0.000,0.500,0.273,0.160,0.143,0.389,0.400,0.000,0.667,0.154,0.250,0.000) D.SO<-c(0.333,0.000,0.250,0.143,1.000,0.053,0.000,0.071,0.250,0.125,0.000,0.200,0.182,0.125,0.125,0.143,0.455,0.250,0.000,0.133,0.143,0.667,0.250,0.250,0.333,0.385,0.200,0.250,0.150,0.333,0.333,0.111,0.059,1.000,0.333,0.538,0.000,0.267,0.545,0.179,0.095,0.077,0.261,0.500,0.667,0.667,0.231,0.357,0.118,0.190,0.429,0.231,0.000,0.211,0.091,0.286,0.176,0.214,0.067,0.500,0.500,0.182,0.280,0.214,0.444,0.200,0.167,0.333,0.231,0.500,0.000) S.SO<-c(0.000,0.083,0.000,0.000,0.000,0.000,0.111,0.000,0.000,0.125,0.000,0.000,0.000,0.000,0.063,0.000,0.000,0.000,0.000,0.067,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.036,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.533,0.000,0.000,0.000,0.000,0.000,0.000,0.200,0.000,0.000,0.000,0.000,0.000) O.SO<-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.056,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.077,0.000,0.000,0.048,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000) Rep.LA<-c(0.333,0.857,0.333,0.833,0.667,0.750,0.538,0.600,0.000,0.429,0.765,0.786,1.000,0.818,0.286,0.833,0.500,0.857,0.500,0.583,0.600,0.833,0.636,0.417,0.800,0.000,0.375,0.476,0.143,0.357,0.200,0.333,0.304,0.385,0.688,0.143,0.167,0.375,0.250,0.167,0.682,0.000,0.000,0.500,0.000,0.385,0.435,0.333,0.826,0.000,0.417,0.842,0.643,0.750,0.429,0.286,0.000,0.409,0.444,0.000,0.000,0.536,0.667,0.636,0.632,0.500,0.600,0.000,0.444,1.000,0.000) Ali.LA<-c(0.667,0.000,0.333,0.083,0.111,0.083,0.077,0.300,0.000,0.000,0.000,0.000,0.000,0.000,0.429,0.000,0.000,0.000,0.250,0.000,0.200,0.083,0.091,0.333,0.000,1.000,0.375,0.429,0.714,0.357,0.400,0.400,0.478,0.231,0.063,0.286,0.667,0.125,0.583,0.667,0.182,0.714,0.000,0.167,0.000,0.308,0.435,0.500,0.130,0.000,0.417,0.053,0.000,0.167,0.500,0.143,0.750,0.227,0.000,0.333,0.000,0.214,0.111,0.182,0.158,0.000,0.000,0.000,0.222,0.000,0.000) D.LA <-c(0.000,0.000,0.222,0.000,0.000,0.167,0.077,0.100,0.000,0.571,0.176,0.214,0.000,0.091,0.286,0.167,0.500,0.143,0.250,0.333,0.200,0.083,0.273,0.250,0.200,0.000,0.188,0.095,0.143,0.286,0.400,0.267,0.217,0.308,0.250,0.571,0.167,0.500,0.167,0.167,0.091,0.286,0.000,0.167,0.000,0.308,0.130,0.167,0.043,1.000,0.167,0.105,0.071,0.083,0.071,0.571,0.250,0.364,0.167,0.667,1.000,0.214,0.222,0.182,0.211,0.000,0.400,1.000,0.333,0.000,0.000) S.LA <-c(0.000,0.143,0.111,0.083,0.222,0.000,0.308,0.000,0.000,0.000,0.059,0.000,0.000,0.091,0.000,0.000,0.000,0.000,0.000,0.083,0.000,0.000,0.000,0.000,0.000,0.000,0.063,0.000,0.000,0.000,0.000,0.000,0.000,0.077,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.167,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.071,0.000,0.000,0.000,0.000,0.000,0.389,0.000,0.000,0.036,0.000,0.000,0.000,0.500,0.000,0.000,0.000,0.000,0.000) O.LA <-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.045,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.214,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000) Rep.PE<-c(0.000,0.000,0.000,0.000,0.429,0.600,0.667,0.650,0.909,0.600,0.688,0.875,0.636,0.929,0.200,0.300,1.000,0.000,0.286,0.375,0.000,0.818,0.500,0.500,0.250,0.000,0.308,0.000,0.500,0.000,0.400,0.231,0.125,0.333,0.182,0.500,0.600,0.286,0.182,0.692,0.800,0.182,0.750,0.778,0.556,0.125,0.238,0.000,0.846,0.000,0.357,0.565,0.333,0.583,0.350,0.667,0.273,0.407,0.300,0.684,0.143,0.368,0.429,0.316,0.600,0.000,1.000,0.000,0.667,0.579,0.500) Ali.PE<-c(0.000,0.000,0.000,0.000,0.143,0.050,0.000,0.300,0.000,0.000,0.000,0.125,0.227,0.071,0.600,0.600,0.000,0.000,0.143,0.250,0.500,0.091,0.333,0.200,0.250,1.000,0.308,0.000,0.500,0.500,0.600,0.538,0.250,0.444,0.545,0.500,0.000,0.429,0.545,0.154,0.133,0.455,0.000,0.111,0.222,0.250,0.476,0.667,0.077,1.000,0.429,0.217,0.000,0.167,0.400,0.000,0.455,0.333,0.000,0.158,0.571,0.263,0.333,0.368,0.250,0.333,0.000,0.000,0.222,0.211,0.250) D.PE <-c(0.000,0.000,0.000,0.000,0.143,0.050,0.000,0.000,0.091,0.200,0.000,0.000,0.000,0.000,0.200,0.100,0.000,0.000,0.429,0.375,0.250,0.091,0.167,0.200,0.500,0.000,0.154,0.000,0.000,0.500,0.000,0.154,0.375,0.222,0.000,0.000,0.200,0.286,0.273,0.154,0.067,0.364,0.250,0.000,0.222,0.625,0.190,0.333,0.077,0.000,0.214,0.174,0.000,0.250,0.200,0.333,0.273,0.259,0.000,0.158,0.286,0.316,0.190,0.211,0.150,0.333,0.000,0.667,0.111,0.211,0.000) S.PE <-c(0.000,0.000,0.000,0.000,0.000,0.300,0.333,0.050,0.000,0.200,0.313,0.000,0.136,0.000,0.000,0.000,0.000,0.000,0.143,0.000,0.250,0.000,0.000,0.100,0.000,0.000,0.231,0.000,0.000,0.000,0.000,0.077,0.250,0.000,0.273,0.000,0.200,0.000,0.000,0.000,0.000,0.000,0.000,0.111,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.043,0.000,0.000,0.050,0.000,0.000,0.000,0.700,0.000,0.000,0.053,0.048,0.105,0.000,0.333,0.000,0.333,0.000,0.000,0.250) O.PE <-c(0.000,0.000,0.000,0.000,0.286,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.095,0.000,0.000,0.000,0.000,0.000,0.667,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000) Rep.BR<-c(0.000,0.000,0.000,0.364,0.000,0.000,0.900,0.619,0.591,0.500,0.000,0.600,0.500,0.000,0.333,0.700,0.750,0.000,0.333,0.364,0.333,0.000,0.368,0.474,0.300,0.000,0.333,0.471,0.368,0.273,0.429,0.407,0.143,0.500,0.000,0.714,0.000,0.200,0.500,0.400,0.846,0.000,0.000,0.000,0.571,0.111,0.357,0.333,0.000,0.000,0.467,0.333,0.000,0.421,0.333,0.200,0.000,0.909,0.333,0.600,0.000,0.438,0.111,0.750,0.667,0.000,0.556,0.000,1.000,1.000,0.000) Ali.BR<-c(0.000,0.000,0.000,0.182,0.333,1.000,0.000,0.286,0.182,0.000,0.000,0.100,0.083,0.000,0.467,0.200,0.125,0.000,0.167,0.409,0.556,0.000,0.368,0.474,0.400,1.000,0.333,0.294,0.368,0.455,0.429,0.444,0.714,0.167,0.750,0.143,0.000,0.400,0.375,0.400,0.154,0.500,0.000,0.000,0.214,0.222,0.571,0.333,0.000,0.000,0.467,0.333,0.727,0.421,0.500,0.400,0.000,0.000,0.000,0.200,0.000,0.188,0.667,0.000,0.200,0.818,0.000,0.000,0.000,0.000,0.000) D.BR <-c(0.000,0.000,0.000,0.182,0.333,0.000,0.000,0.048,0.227,0.500,0.000,0.200,0.167,0.000,0.133,0.100,0.125,0.000,0.500,0.136,0.111,0.000,0.211,0.053,0.300,0.000,0.167,0.176,0.211,0.273,0.143,0.148,0.143,0.167,0.250,0.143,0.000,0.400,0.000,0.200,0.000,0.500,0.000,1.000,0.214,0.667,0.071,0.333,0.000,0.000,0.067,0.333,0.182,0.158,0.111,0.400,0.000,0.091,0.083,0.200,1.000,0.313,0.222,0.083,0.133,0.091,0.444,0.000,0.000,0.000,0.000) S.BR <-c(0.000,0.000,1.000,0.273,0.333,0.000,0.100,0.048,0.000,0.000,0.000,0.100,0.250,0.000,0.067,0.000,0.000,0.000,0.000,0.091,0.000,0.000,0.053,0.000,0.000,0.000,0.167,0.059,0.053,0.000,0.000,0.000,0.000,0.167,0.000,0.000,0.000,0.000,0.125,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.091,0.000,0.056,0.000,0.000,0.000,0.583,0.000,0.000,0.063,0.000,0.167,0.000,0.091,0.000,0.000,0.000,0.000,0.000) O.BR <-c(1.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000) Rep.ZU<-c(0.000,0.000,0.000,0.000,0.000,0.800,0.556,0.250,0.000,0.000,0.700,0.500,0.000,0.000,0.000,0.750,0.800,0.000,0.000,0.375,0.565,0.000,0.111,0.000,0.000,0.000,0.333,0.000,0.000,0.375,0.000,0.143,0.333,0.500,0.333,0.250,0.000,0.000,0.500,0.000,0.643,0.375,0.875,0.250,0.000,0.000,0.000,0.167,0.818,0.429,0.000,0.600,0.000,0.000,0.000,0.667,0.000,0.000,0.000,0.000,0.000,0.250,0.000,0.231,0.833,0.000,0.000,0.000,0.938,0.000,0.000) Ali.ZU<-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.750,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.250,0.174,0.000,0.556,0.000,0.000,1.000,0.167,1.000,0.500,0.375,0.000,0.857,0.067,0.400,0.333,0.750,0.000,0.500,0.250,1.000,0.143,0.000,0.000,0.000,0.000,0.000,0.500,0.167,0.000,0.143,0.000,0.000,0.000,0.667,1.000,0.000,0.000,0.000,0.000,0.000,0.000,0.125,0.000,0.308,0.000,0.500,1.000,0.000,0.000,0.000,0.000) D.ZU<-c(0.000,0.000,0.000,0.000,0.000,0.100,0.111,0.000,0.000,0.000,0.100,0.250,0.000,0.000,1.000,0.250,0.200,0.000,0.000,0.375,0.217,0.000,0.333,0.000,0.000,0.000,0.167,0.000,0.500,0.125,1.000,0.000,0.533,0.000,0.333,0.000,1.000,0.500,0.250,0.000,0.143,0.375,0.125,0.750,1.000,0.000,0.500,0.667,0.182,0.429,0.000,0.400,0.000,0.333,0.000,0.333,0.000,0.000,0.000,0.000,1.000,0.625,1.000,0.231,0.167,0.500,0.000,0.000,0.063,0.000,0.000) S.ZU<-c(0.000,0.000,0.000,0.000,0.000,0.100,0.333,0.000,0.000,0.000,0.200,0.250,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.043,0.000,0.000,0.000,0.000,0.000,0.333,0.000,0.000,0.125,0.000,0.000,0.067,0.100,0.000,0.000,0.000,0.000,0.000,0.000,0.071,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.231,0.000,0.000,0.000,0.000,0.000,0.000,0.000) O.ZU<-c(0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.250,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000) # Funções para criação das tabelas dos comportamentos para cada um dos indivíduos. tabela.repouso<-cbind(Rep.RE,Rep.SO,Rep.LA,Rep.PE,Rep.BR,Rep.ZU) tabela.alimentacao<-cbind(Ali.RE,Ali.SO,Ali.LA,Ali.PE,Ali.BR,Ali.ZU) tabela.deslocamento<-cbind(D.RE,D.SO,D.LA,D.PE,D.BR,D.ZU) tabela.social<-cbind(S.RE,S.SO,S.LA,S.PE,S.BR,S.ZU) tabela.outros<-cbind(O.RE,O.SO,O.LA,O.PE,O.BR,O.ZU) Espero que tenha dado para entender as minhas dúvidas e os meus dados. Agradeço desde já. Abraços, -- Flávia Machado

2 1

ANOVA em regressão não linear!
by Thiago de Paula Protásio 14 Jan '12

14 Jan '12

Prezados membros da R-Br; Recorro a valiosa ajuda de vcs! Estou trabalhando com os seguintes modelos não lineares para modelar o acúmulo de biomassa com o tempo: Yi = A*exp(B*X) ln(Yi) = A * (1 - exp(-B*X) ln (Yi) = A + B*exp(C*X) Nas "Rídiculas" eu encontrei uma função e explicação do Walmes: http://ridiculas.wordpress.com/2011/05/29/analise-de-residuos-em-regressao-… Mas gostaria de saber se posso aplicar essa teoria/função para os modelos que estou trabalhando. Na verdade, estou testando modelos lineares também e gostaria de saber se é fundamental utilizar o Teste F para os modelos não lineares à semelhança dos modelos lineares.?? Cheguei a fazer alguma coisa como apresentado abaixo para um dos modelos. Mas na ANOVA que foi realizada entre o modelo não linear e o modelo reduzido para o intercepto obtive uma soma de quadrado negativa, daí a minha dúvida. Muito obrigado pela valiosa ajuda! -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

2 1

Ajuda para buscar coordenadas de um ponto qq em uma matriz com coordenadas.
by Maria Papa 13 Jan '12

13 Jan '12

Olá pessoal, boa tarde, Eu estou com outro problema! Eu tenho um conjunto de coordenadas, (X,Y,Z) uma matriz com uma coluna de valores para X, outra para Y e outra para Z. Dado uma coordenada (XYZ) qualquer, eu preciso buscar neste conjunto (matriz) os valores imediatamente maiores e menores que as coordenadas deste ponto qualquer. Eu pensei no comando ifelse, mas eu teria que fazer para cada coordenadas. Alguém saberia me dizer se existe alguma forma mais simples para fazer isso? Desde já agradeço a todos. célia

2 3

Função para somar pontuação
by Carlos Andrade 12 Jan '12

12 Jan '12

Prezados amigos da lista Que 2012 seja repleto de sucesso. Minha dúvida. Elaborei um questionário e atribuí pontos para as diferentes respostas. Como codifico uma função para gerar uma variável que será o resultado da computação da soma dos pontos de cada indivíduo de acordo com suas respostas? Abaixo segue o esquema de pontuação das questões que entram no cálculo. Pergunta Resposta Pontuação p1 1 5 2 1 3 0 p2 1 5 2 3 3 3 4 3 5 1 6 0 p3.1 1 1 0 0 P3.2 1 3 0 0 P3.3 1 3 0 0 P3.4 1 5 0 0 P3.5 1 5 0 0 P3.6 1 3 0 0 P3.7 1 0 0 0 p4 1 1 2 3 3 3 4 3 5 5 6 0 P5.1 1 3 0 0 P5.2 1 1 0 0 P5.3 1 3 0 0 P5.4 1 1 0 0 P5.5 1 1 0 0 P5.6 1 5 0 0 P5.7 1 1 0 0 P5.8 1 1 0 0 P5.9 1 1 0 0 P5.10 1 0 0 0 p6 1 1 2 3 3 3 4 3 5 5 6 0 P7.1 1 5 0 0 P7.2 1 5 0 0 P7.3 1 5 0 0 P7.4 1 5 0 0 P7.5 1 5 0 0 P7.6 1 3 0 0 P7.7 1 0 0 0 P8.1 1 5 0 0 P8.2 1 3 0 0 P8.3 1 5 0 0 P8.4 1 1 0 0 P8.5 1 0 0 0 P9.1 1 5 0 0 P9.2 1 5 0 0 P9.3 1 5 0 0 P9.4 1 5 0 0 P9.5 1 5 0 0 P9.6 1 5 0 0 P9.7 1 0 0 0 P10.1 1 5 0 0 P10.2 1 5 0 0 P10.3 1 5 0 0 P10.4 1 5 0 0 P10.5 1 5 0 0 P10.6 1 0 0 0 P10.7 1 0 0 0 p11 1 5 2 3 3 1 4 0 5 0 p15 1 5 2 3 3 3 4 1 5 1 6 0 Também: como substituir os NAs por zeros de uma vez só? Agradeço antecipadamente qualquer ajuda. -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

3 2

Gerar multinomial
by Mauro Sznelwar 12 Jan '12

12 Jan '12

Gostaria de saber onde posso achar comandos do R para geração de números aleatórios da distribuição multinomial. Não posso usar o comando rmultinom que existe, preciso implementar algo para três variáveis, teta1, teta2 e teta3. Desde já agradeço.

5 11

Re: [R-br] Superfície de Resposta (Resolvido)
by Evelyn Alecrim 11 Jan '12

11 Jan '12

Pessoal, obrigada pela ajuda. Consegui fazer o gráfico usando persp(). -- Evelyn da Fonseca Alecrim Bragion Mestranda em Ecologia Aplicada - UFLA Lavras - MG/Brasil

1 0

Ajuda com NAs
by Daniel C Bezerra 11 Jan '12

11 Jan '12

Amigos do R, Eu tenho uma variável que foi medida três vezes no tempo com NAs em pontos diferentes de cada um dos momentos em que foi medida: Eu gostaria de fazer a média destas 3 medidas mas caso uma medida fosse NA eu não gostaria de ignorar esta linha e torná-la NA e sim usar a média das medidas restantes. Alguém teria uma sugestão? Exemplo: m1<-c(1,2,3,4,NA,NA) m2<-c(1,NA,3,4,5,NA) m3<-c(NA,2,3,NA,NA,NA) Se eu faço apply(cbind(m1,m2,m3),1,mean) Eu tenho: [1] NA NA 3 NA NA NA Mas eu gostaria de: [1] 1 2 3 4 5 NA Muito obrigado, D

3 3

Anova de um modelo.
by Julia Maria 11 Jan '12

11 Jan '12

Pessoal, enviai um trabalho pra um Congresso e voltou pedindo algumas correções, dentre elas pedindo Anova dos modelos que eu utilizei. Como estou descobrindo o R agora, preciso da ajuda de voces. A partir do modelo: Y=a+b+c+ab+ac+bc, eu consigo recuperar uma anova? Tenho os valores de a, b e c e tb tenho os Y originais e os preditos. Bastaria escrever o modelo e pedir a anova?

3 5

Superfície de Resposta
by Evelyn Alecrim 11 Jan '12

11 Jan '12

Olá Pessoal, Estou precisando produzir um gráfico de superfície para o seguinte modelo : lm(y=a+bx+cz+dx*z) Alguém poderia me dizer que pacote poderia usar e qual o código? Obrigada... -- Evelyn da Fonseca Alecrim Bragion Mestranda em Ecologia Aplicada - UFLA Lavras - MG/Brasil

4 5

Resultado gráfico do teste de kruskal-Wallis
by Maria Papa 10 Jan '12

10 Jan '12

Olá Pessoal, boa tarde, Eu gostaria de saber da possibilidade de gerar um resultado gráfico do teste de Kruska-Wallis, assim como o Tukey apresenta o gráfico da diferença das médias. A comparação com o Tukey não tem sentido teórico, é apenas para exemplificar o tipo de resultado que eu preciso. Eu estou usando a função kruskal do pacote agricolae. Desde já agradeço a todos Maria Célia

2 1