The R Project SVN R

Rev

Rev 26189 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 26189 Rev 26680
Line 42... Line 42...
42 
 
 42 
 
43 
#define COV_PAIRWISE_BODY						\
 43 
#define COV_PAIRWISE_BODY						\
44 
	    xx = &x[i * n];						\
 44 
	    xx = &x[i * n];						\
45 
	    nobs = 0;							\
 45 
	    nobs = 0;							\
46 
	    if(!kendall) {						\
 46 
	    if(!kendall) {						\
- 
 
 47 
	    	xmean = ymean = 0.;					\
47 
		for (k = 0 ; k < n ; k++) {				\
 48 
		for (k = 0 ; k < n ; k++) {				\
48 
		    if(!(ISNAN(xx[k]) || ISNAN(yy[k]))) {		\
 49 
		    if(!(ISNAN(xx[k]) || ISNAN(yy[k]))) {		\
49 
			nobs ++;					\
 50 
			nobs ++;					\
50 
			xmean += xx[k];					\
 51 
			xmean += xx[k];					\
51 
			ymean += yy[k];					\
 52 
			ymean += yy[k];					\
Line 144... Line 145...
144 
 
 145 
 
145 
#define ANS(I,J)  ans[I + J * ncx]
 146 
#define ANS(I,J)  ans[I + J * ncx]
146 
 
 147 
 
147 
#define COV_init(_ny_)				\
 148 
#define COV_init(_ny_)				\
148 
    double sum, xxm, yym, *xx, *yy;		\
 149 
    double sum, xxm, yym, *xx, *yy;		\
149 
    int i, j, k, nobs, n1;			\
 150 
    int i, j, k, nobs, n1=-1;/* -Wall */	\
150 
						\
 151 
						\
151 
    /* total number of complete observations */	\
 152 
    /* total number of complete observations */	\
152 
    nobs = 0;					\
 153 
    nobs = 0;					\
153 
    for(k = 0 ; k < n ; k++) {			\
 154 
    for(k = 0 ; k < n ; k++) {			\
154 
	if (ind[k] != 0) nobs++;		\
 155 
	if (ind[k] != 0) nobs++;		\