[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