Dear Rhelp,
 
I did a daisy on 5 lifestyle variables, 3 of which were nominal and 2 were 
ordinal and assigned types “nominal” and “ordinal” for the variables, 
respectively.  I got an output indicating their types as “I” for 
interval(?). Doing it on the Rdata example “flower” gave the same types in 
the output as the types they were assigned to.  Why is this so? Below are the 
codes and outputs.
 
sfq    is a nominal variable with 5 categories pertaining to
       smoking frequency and consumption (1=none, 2=≤10
       sticks/day somedays, 3=>10 sticks/day somedays, 4=≤10
       sticks/day daily, 5=>10 sticks/day daily)
afq    is a nominal variable with 5 categories pertaining to
       alcohol frequency and consumption
pafd   is a nominal variable with 5 categories pertaining to
       physical activity frequency and duration
dietp1 is an ordinal variable with 3 categories pertaining low, 
       medium, high consumption of Western diet
dietp3 is an ordinal variable with 3 categories pertaining to
       low, medium, high consumption of prudent diet
   
 
>head(lsclusjt3)
  sfq afq pafd dietp1 dietp3
1   1   1    3      1      2
2   1   1    3      3      3
3   1   1    1      2      1
4   1   1    1      2      2
5   1   2    3      3      3
6   1   1    1      2      2
 
>dm=daisy(lsclusjt3,metric="gower", stand=FALSE,type=list(nominal=c(1,2,3), 
>ordinal=c(4,5)))
>summary(dm)
38434528 dissimilarities, summarized :
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
0.00000 0.25000 0.35000 0.36599 0.50000 1.00000 
Metric :  mixed ;  Types = I, I, I, I, I
Number of objects : 8768
 
>dfl=daisy(flower,type=list(asymm=1:3,nominal=4,ordinal=5:6,interval=7:8))
> summary(dfl)
153 dissimilarities, summarized :
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
0.15915 0.43576 0.53408 0.53473 0.62908 0.89099 
Metric :  mixed ;  Types = A, A, A, N, O, O, I, I 
Number of objects : 18
 
 
Penny
 


      
        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to