The R Project SVN R

Rev

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

Rev 61168 Rev 61434
Line 165... Line 165...
165 
}
 165 
}
166 
par(mfrow = c(1, 1))
 166 
par(mfrow = c(1, 1))
167 
 
 167 
 
168 
## "Extended" Pascal triangle:
 168 
## "Extended" Pascal triangle:
169 
fN <- function(n) formatC(n, width=2)
 169 
fN <- function(n) formatC(n, width=2)
- 
 
 170 
for (n in -4:10) {
170 
for (n in -4:10) cat(fN(n),":", fN(choose(n, k = -2:max(3, n+2))), "\n")
 171 
    cat(fN(n),":", fN(choose(n, k = -2:max(3, n+2))))
- 
 
 172 
    cat("\n")
- 
 
 173 
}
171 
 
 174 
 
172 
## R code version of choose()  [simplistic; warning for k < 0]:
 175 
## R code version of choose()  [simplistic; warning for k < 0]:
173 
mychoose <- function(r, k)
 176 
mychoose <- function(r, k)
174 
    ifelse(k <= 0, (k == 0),
 177 
    ifelse(k <= 0, (k == 0),
175 
           sapply(k, function(k) prod(r:(r-k+1))) / factorial(k))
 178 
           sapply(k, function(k) prod(r:(r-k+1))) / factorial(k))