## 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