similar to: linking r and C

hello all i searched the archives and couldn't get a solution to the following question. i have the following function: F=function(z,v) { if (v==-.5) return(1) else return(2*v/z + 1/Recall(z,v-1)) } and while testing whether it works (ie F(z,v) is approx = besselK(z,1+v)/besselK(z,V). the recursion formula allows one to calculate besselK(z,1+v)/besselK(z,V) for large values of z )

hello all sorry for the following "none" R related question. does anyone know of a reference to calculate the following identity: |I + ABC| where I is an identity matrix and A, B,C may not have to be square matrices? you help will be greatly appreciated. H. V. Henderson; S. R. Searle SIAM Review, Vol. 23, No. 1. (Jan., 1981), pp. 53-60. provides a result to

hello R users strange behavior of the integrate function! i assume this occurs because of the way in which the quadriture is set up! (any comments.) f=function(x){exp(-exp(-x)-5*x)/gamma(5)} xx=seq(from=-20, to=20, length.out=1000) plot(xx,f(xx),type="l") integrate(f, lower=-Inf, upper= 1) integrate(f, lower=-Inf, upper= 10) integrate(f, lower=-10, upper= 11) integrate(f,

hello all i searched the archives and couldn't get a solution to the following question. i have the following function: F=function(z,v) { if (v==-.5) return(1) else return(2*v/z + 1/Recall(z,v-1)) } and while testing whether it works (ie F(z,v) is approx = besselK(z,1+v)/besselK(z,V). the recursion formula allows one to calculate besselK(z,1+v)/besselK(z,V) for large values of z )

hello all happy new year and hope you r having a good holiday. i would like to calculate the expectation of a particular random variable and would like to approximate it using a number of the functions contained in R. decided to do some experimentation on a trivial example. example ======== suppose x(i)~N(0,s2) where s2 = the variance the prior for s2 = p(s2)~IG(a,b) so the posterior is

I'm testing some different formulations of pairwise interaction point processes in Spatstat (version 1.11-6) using R 2.5.0 on a Windows platform and I wish to simulate them using the Metropolis-Hastings algorithm implemented with Spatstat. Spatstat utilizes Fortran77 code with the preprocessor RatFor to do the Metropolis-Hastings MCMC, but the Makefile is more complicated than any I have

Hi I know how to referenc R in a scientific paper - but is there a standardised way to reference packages? Thanks Rainer -- NEW EMAIL ADDRESS AND ADDRESS: Rainer.Krug at uct.ac.za RKrug at sun.ac.za WILL BE DISCONTINUED END OF MARCH Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation Biology (UCT) Plant Conservation Unit Department of Botany University of Cape Town Rondebosch 7701

Dear list, After running for a while, it crashes and gives the following error message: can anybody suggest how to deal with this? Error in if (ratio0[i] < log(runif(1))) { : missing value where TRUE/FALSE needed ################### original program ######## p2 <- function (Nsim=1000){ x<- c(0.301,0,-0.301,-0.602,-0.903,-1.208, -1.309,-1.807,-2.108,-2.71) # logdose

Hi Consider the following code > x <- rep(1:13, 13) > y <- 1:3 I want to select all elements in x which are equal to 1, 2 or 3. I know that I could use > sel <- x==y[1] | x==y[2] | x==y[3] > x[sel] to obtain the values, but in my analysis, the y-vector is thousands of elements long. Is there any way, that I can do that easily? Thanks Rainer -- Rainer M. Krug, Dipl.

Hi I am trying to install the package "pheno", but it needs the package "nprq" by Roger Koenker et al. which I can I find this package? It does not seem to be on CRAN and googling also doesn't give me an URL - is it still somewhere available? Thanks, Rainer -- NEW EMAIL ADDRESS AND ADDRESS: Rainer.Krug at uct.ac.za RKrug at sun.ac.za WILL BE DISCONTINUED END OF

Dear all, The new version of AdMit (version 1.01-01) is now available from CRAN. SUMMARY The package provides functions to perform the fitting of an adaptive mixture of Student-t distributions to a target density through its kernel function. The mixture approximation can then be used as the importance density in importance sampling or as the candidate density in the Metropolis-Hastings

Dear all, The new version of AdMit (version 1.01-01) is now available from CRAN. SUMMARY The package provides functions to perform the fitting of an adaptive mixture of Student-t distributions to a target density through its kernel function. The mixture approximation can then be used as the importance density in importance sampling or as the candidate density in the Metropolis-Hastings

Hi I am sure this is simple - but how can I convert one charecter vector into one string? example: x <- c("This ", "is ", "one ", "sentence.") should become "This is one entence" Thanks Rainer -- NEW EMAIL ADDRESS AND ADDRESS: Rainer.Krug at uct.ac.za RKrug at sun.ac.za WILL BE DISCONTINUED END OF MARCH Rainer M. Krug, Dipl. Phys.

Hi I created a complex layout with using layout() and it looks exactly as I need it. But I don't want to print in the order in which the subfigure are numbered, but in a different order. How can I navigate in the layout so that I can specify the subfigure in which to plot? At the moment I am using a function which is converting the number to mfg parameter for par, but it does not seem to

Hi I tried to install rkward under ubuntu hardy heron, but it tried to use the one from the cran repository which was newer, but it did not install. To be able to install rkward, I had to disable the cran repository, install rkward, lock it's version and enable the repository again. I have the feeling, that the dependencies are not correct (it is looking for libqt4GUI, but I have a

Hi all, I have some wind data (U and V components) and I would like to compute Vorticity and Divergence of these fields. Is there any R function that can easily do that? Thanks in advance for any help Igor Oliveira CSAG, Dept. Environmental & Geographical Science, University of Cape Town, Private Bag X3, Rondebosch 7701. Tel.: +27 (0)21 650 5774 South Africa Fax: +27 (0)21

Hi Sorry for crossposting, but I think this can be of interest for GRASS and R users. I am planning to write a package to make the use of GRASS from R easier. The idea is to wrap the system call to execute the GRASS command into an R command of the same name. e.g: r.to.vect <- function(..., intern=TRUE, ignore.stderr=FALSE) { comm <- paste( "r.to.vect ", ...,

Hi all, I'm trying to write a MH algorithm in R for a standard normal distribution, I've been trying for a good week or so now with multiple attempts and have finally given up trying to do it on my own as I'm beginning to run out of time for this, would somebody please tell me what is wrong with my latest attempt: n=100 mu=0 sigma=1 lik<-function(theta) exp(((theta-mu)^2)/2*sigma)

Hi I have a large dataset which follows a multimodal distribution. And I would like to identify the maxima. As the data is obtained from a stochastic simulation, not all maxima in the data are "real maxima of the dirstribution" but rather small random fluctuations. Unfortunately, it is not possible for me to run more simulations to smooth the obtained distribution. What I am doing

Hi I just installed the pastecs package and I am wondering: is there an english (or german) translation of the file pastecs.pdf? If not, is there an explanation somewhere of the object of type 'turnpoints' as a result of turnpoints(), especially the "info" field? Thanks, Rainer -- NEW EMAIL ADDRESS AND ADDRESS: Rainer.Krug at uct.ac.za RKrug at sun.ac.za WILL BE