## calling C functions from R
gcc -c kendall.c
R CMD SHLIB -o libkendall.so kendall.o
R
dyn.load('libkendall.so')
xd <- c(44.4, 45.9, 41.9, 53.3, 44.7, 44.1, 50.7, 45.2, 60.1)
yd <- c( 2.6, 3.1, 2.5, 5.0, 3.6, 4.0, 5.2, 2.8, 3.8)
.C("kendall",as.double(xd),as.double(yd),as.integer(9),
result=as.integer(1))$result
## header
rm(list=ls())
options(contrasts=c("contr.treatment","contr.poly"))
setwd("/mnt/win_c/my documents/stat/ammonia")
## plot function
plot(fun,-5, 5, ylim = c(-.2, 1))
eval.parent(expr)
## example eval.parent(media <- 1), used inside
## functions to change global variables
media <<- 1 ## does the same
rep(x,n) replicates x n-times
mtext("texto",side=3,line=.5)
install.packages("fields","/aa/haky/R/Library")
library("fields",lib="/aa/haky/R/Library")
detach(package:fields)
setenv R_LIBS /aa/haky/R/Library (then no need to specify lib=...)
xmat <- matrix(NA,n,n)
#print
dev.print()
dev.copy2eps(file="filename.ps")
#
postscript(file="plot.ps",horizontal=F,heigth=9,width=7)
dev.off()
as.numeric(x<1)
# stop watch
tic<-Sys.time()
....do something here
toc<-Sys.time()
print(toc-tic)
paste(s1,s2,sep="")
system("mail -s something haky@galton.uchicago.edu < output")<
## strip an object's attributes:
attributes(x) <- NULL
x # now just a vector of length 6
kk[names(kk)=="site"] is same as kk$site
unlist(x)
debug(myfunction)
undebug(myfunction)
call R from emacs (ESS): alt-x R
par(mfrow=c(3,2))
Wednesday, December 8, 2004
R legacy help galton
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment