In [1]:
# rm(list=ls()) 
options(OutDec = ",")
#==============================================================================
# Descricao
#
# Analise dos dados de medidas da forca do papel
#
#==============================================================================
silence <- suppressPackageStartupMessages # Para omitir mensagens de alertas
silence(library(SciViews))
silence(library(scatterplot3d))
In [2]:
#==============================================================================
# Lendo os dados
#==============================================================================
x <- read.table("../dados/d01_papel.txt",header=F)
colnames(x) <- c("Densidade","F_direta","F_cruzada")
print(head(x))

  Densidade F_direta F_cruzada
1     0,801   121,41     70,42
2     0,824   127,70     72,47
3     0,841   129,20     78,20
4     0,816   131,80     74,89
5     0,840   135,10     71,21
6     0,842   131,50     78,39
In [3]:
#==============================================================================
# Algumas estatisticas descritivas
#==============================================================================
summary(x)

print("==== Covariâncias ====")
print(var(x))             # matriz de covariancias amostral

print("==== Correlações ====")
print(cor(x))             # matriz de correlacoes amostral

   Densidade         F_direta       F_cruzada    
 Min.   :0,7580   Min.   :103,5   Min.   :48,93  
 1st Qu.:0,7950   1st Qu.:115,1   1st Qu.:56,53  
 Median :0,8150   Median :121,4   Median :70,70  
 Mean   :0,8119   Mean   :121,0   Mean   :67,72  
 3rd Qu.:0,8260   3rd Qu.:126,7   3rd Qu.:74,89  
 Max.   :0,9710   Max.   :135,1   Max.   :80,33  
[1] "==== Covariâncias ===="
            Densidade   F_direta F_cruzada
Densidade 0,001264578  0,1684468  0,225248
F_direta  0,168446762 59,3211480 60,992531
F_cruzada 0,225247976 60,9925314 95,856667
[1] "==== Correlações ===="
          Densidade  F_direta F_cruzada
Densidade 1,0000000 0,6150141 0,6469592
F_direta  0,6150141 1,0000000 0,8088365
F_cruzada 0,6469592 0,8088365 1,0000000
In [4]:
#==============================================================================
# Grafico de dispersao com boxplot na diagonal
#==============================================================================
pairs(x,cex=1.5,pch=16,bg="darkblue",horOdd=T,
    diag.panel=panel.boxplot,cex.labels=1.5,font.labels=1.5)
In [5]:
#==============================================================================
# Grafico de dispersao tridimensional
#==============================================================================
scatterplot3d(x,pch=16,color="red",cex.lab=1.2,cex.axis=1.2,angle=55)