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
binc0Jqi2cXkY.bin
Description: application/vnd.openxmlformats-officedocument.wordprocessingml.document
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
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
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
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
bina3Gr3rPKDp.bin
Description: application/vnd.openxmlformats-officedocument.spreadsheetml.sheet
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
