[R-br] regressão no R
alanarocha em sapo.pt
alanarocha em sapo.pt
Quarta Maio 15 05:08:12 BRT 2013
bom dia,
porque será que esta regressão me está a dar NA
e não me mostra o r multiple nem r adjusted squared?
não percebo o erro do gráfico tambem.
mando em anexo.
envio tambem o ficheiro excel dos dados originais...
obrigada
Ana
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: gráficos_dadosPanakeia 2sem 2011 e 2012.docx
Tipo: application/vnd.openxmlformats-officedocument.wordprocessingml.document
Tamanho: 22912 bytes
Descrição: não disponÃvel
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130515/def67ec7/attachment.docx>
-------------- Próxima Parte ----------
setwd("G:/SIN/Users/Ana Rua/dados não Panakeia e Panakeia")
dados<-read.table("DadosPanakeia2sem2011.txt",comment.char = "%",dec=",",h=T)
dados
#objectivo é:
#verificar se há diferenças estatisticamente significativas entre os custos e os valores de hemoglobina para clinicas Panakeia no 2semestre de 2011.
reg<-lm(as.numeric(dados[1,-1])~as.numeric(dados[2,-1]))
reg
summary(reg)
plot(as.numeric(dados[1,-1])~as.numeric(dados[2,-1]))
#Adicionar a reta da regressão no gráfico:
abline(reg)
resultados:
X2sem2011 JUL AGO SET OUT NOV DEZ
1 custosPanakeia 25,12 27,34 23,93 23,96 33,07 22,22
2 hgbPanakeia 75.4 77.1 77.0 77.2 78.0 76.7
Call:
lm(formula = as.numeric(dados[1, -1]) ~ as.numeric(dados[2, -1]))
Coefficients:
(Intercept) as.numeric(dados[2, -1])
1 NA
Call:
lm(formula = as.numeric(dados[1, -1]) ~ as.numeric(dados[2, -1]))
Residuals:
1 2 3 4 5 6
-1.764e-16 3.528e-17 3.528e-17 3.528e-17 3.528e-17 3.528e-17
Coefficients: (1 not defined because of singularities)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.000e+00 3.528e-17 2.835e+16 <2e-16 ***
as.numeric(dados[2, -1]) NA NA NA NA
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 8.641e-17 on 5 degrees of freedom
-------------- Próxima Parte ----------
setwd("G:/SIN/Users/Ana Rua/dados não Panakeia e Panakeia")
dados<-read.table("DadosPanakeia.txt",comment.char = "%",dec=",",h=T)
dados
#objectivo é:
#verificar se há diferenças estatisticamente significativas entre os custos e os valores de hemoglobina para clinicas Panakeia no 2semestre de 2011.
reg<-lm(as.numeric(dados[1,-1])~as.numeric(dados[2,-1]))
reg
summary(reg)
plot(as.numeric(dados[1,-1])~as.numeric(dados[2,-1]))
#Adicionar a reta da regressão no gráfico:
abline(reg)
resultados:
X2012 JAN FEV MAR ABR MAI JUN JUL AGO SET OUT
1 custosPanakeia 23,39 27,90 23,15 22,65 27,42 21,20 20,81 23,76 20,48 19,36
2 hgbPanakeia 76.2 77.7 76.7 78.4 79.7 79.0 76.7 77.3 77.1 78.7
NOV DEZ
1 28,21 16,41
2 79.0 78.3
> reg<-lm(as.numeric(dados[1,-1])~as.numeric(dados[2,-1]))
> reg
Call:
lm(formula = as.numeric(dados[1, -1]) ~ as.numeric(dados[2, -1]))
Coefficients:
(Intercept) as.numeric(dados[2, -1])
1 NA
> summary(reg)
Call:
lm(formula = as.numeric(dados[1, -1]) ~ as.numeric(dados[2, -1]))
Residuals:
Min 1Q Median 3Q Max
-3.680e-16 3.346e-17 3.346e-17 3.346e-17 3.346e-17
Coefficients: (1 not defined because of singularities)
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.000e+00 3.346e-17 2.989e+16 <2e-16 ***
as.numeric(dados[2, -1]) NA NA NA NA
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 1.159e-16 on 11 degrees of freedom
-------------- Próxima Parte ----------
2sem2011 JUL AGO SET OUT NOV DEZ
custosPanakeia 25,12 27,34 23,93 23,96 33,07 22,22
hgbPanakeia 75.4 77.1 77.0 77.2 78.0 76.7
-------------- Próxima Parte ----------
2012 JAN FEV MAR ABR MAI JUN JUL AGO SET OUT NOV DEZ
custosPanakeia 23,39 27,90 23,15 22,65 27,42 21,20 20,81 23,76 20,48 19,36 28,21 16,41
hgbPanakeia 76.2 77.7 76.7 78.4 79.7 79.0 76.7 77.3 77.1 78.7 79.0 78.3
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: Dados em análise.xlsx
Tipo: application/vnd.openxmlformats-officedocument.spreadsheetml.sheet
Tamanho: 10995 bytes
Descrição: não disponÃvel
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130515/def67ec7/attachment.xlsx>
Mais detalhes sobre a lista de discussão R-br