similar to: R: lags

Hi, I've written some C code for density evaluation of the non-central t distribution. It works with the R-1.7.1 source code if placed in the src/nmath directory and after appropriate changes are made to Makefiles, to the dt function in src/library/base/R/distn.R etc. I haven't read a lot of R source code so it may need some R-ification, but I've tried to use the dt.c file as a

Hi all I want to calculate certain lags of a time series and plot them simultaneously on a graph. can anyone help?

Hi all this is a silly question since i should know the answer. lm(y~x) perfroms linear regression with the intercept included. How do i estimate the equation without the intercept? cheers

> On Dec 15, 2014, at 12:21 PM, Kurt Hornik <Kurt.Hornik at wu.ac.at> wrote: > >>>>>> Spencer Graves writes: > >> Hello, All: >> What would it take to make ?iconv? portable? > > >> I ask, because I want to convert accented characters to >> vanilla ASCII, thereby converting, e.g., ?Ra?l? to ?Raul?, and >> Milan

On Dec 15, 2014, at 1:37 PM, Spencer Graves <spencer.graves at prodsyse.com> wrote: > > >> On Dec 15, 2014, at 10:13 AM, Simon Urbanek <simon.urbanek at r-project.org> wrote: >> >>> >>> On Dec 15, 2014, at 12:21 PM, Kurt Hornik <Kurt.Hornik at wu.ac.at> wrote: >>> >>>>>>>> Spencer Graves writes:

Hi, I've made a selfStart function for use with gnls and the following piece of code works nicely: check1 <- gnls(y ~ spot.shape.fct(xcord, ycord, background, spotintensity, rho, sigma, delta, mux, muy), start=getInitial(y ~ spot.shape.fct(xcord, ycord, background, spotintensity, rho,

a simple question how does one produce plots on two different graphics devices? / allan

hi all a simple question i want to run simulations in r. i however want the experiments to be repeated at a later time with exactly the same numbers by other users. can i set the random number seed for rnorm in some way? e.g. is there some arguement that goes with rnorm? please supply an example regards Allan

Greetings, I've tried to update the power.t.test function to allow for different sample sizes and sample variances in the case of a two-sample t test. I'd gladly update the corresponding help page if the code is to replace the current power.t.test function. The modified power.t.test code is included below. The changes are as follows: - Added three new arguments to the function *

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 )

hi all is there a package that undertakes subset selection but BASED ON AIC or any other information criteria. i've seen the "subselect" and the "leaps" package but i have not played around with them yet. thanx

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

hi all i have a simple question: lets assume that i can enter a variable (a) equal to some value, say 100. obviously a could be any value. i would like to create a plot with a main heading as hi all "some writing" + a how would own create a string by combining a string with a numeric variable? hope the question is not too confusing. e.g. say a=5 then i want something like

hi all hopefully some one can help. assume that i imported the following data into R (say the data frame is called a) x1 x2 x3 1 NA 3 1 2 NA 1 2 3 3 NA 6 4 5 9 7 5 6 7 8 9 NA 7 9 How can i construct a new data frame that only contains those rows that does not contain the NA's? is these a quick way? ie x1 x2 x3 1 2 3 4 5 9 7 5 6 7 8 9 in this example we can simple use

hi all why does R do this: (-8)^(1/3)=NaN the answer should be : -2 a silly question but i kept on getting errors in some of my code due to this problem. i solve the problem as follows: say we want : (-a)^(1/3) then : sign(a)*(a^(1/3)) works but there has to be a simpler way of soing such a simple mathematical operation. thanking you / allan

hi all i have another probably simple question. I have three variables say x, y and z. x and y are quite large and z is relative small. how can one plot the three variables on the same graph with two separate axis? (one for x and y and the other for z) e.g. x<-c(101,110,150,167,120) y<-c(120,135,175,95,200) z<-c(0.001, 0.15, 0.6, 0.8, 1) regards Allan

how does one numerically intergate the following: A=function(x,y) { xy } over the range: 2<x<0 4<y<10 say. ie how would one set up the integrate function? i forgot!

Hi all Another simple question. I would like to plot three graphs one the same plot with different colours. Say red, blue and black. Here are the functions. r1<-1+5*cos(2*pi*seq(1:100)/20)+rnorm(100) r2<-1+7*sin(2*pi*seq(1:100)/20)+rnorm(100) r3<-1+7*sin(2*pi*seq(1:100)/20)+5*cos(2*pi*seq(1:100)/20)+rnorm(100) Regards Allan

HI sorry to be a nuisance to all!!! how can i see the code of a particular function? e.g. nnet just as an example

hi all i have a another stats question. i would like to solve the following question: y(i)=a+b*x(i)+e(i) i.e. estimate a and b (they should be fixed) but i dont want to specify the standard density to the straight line. this can be done using kernel regression. the fitted line is however fitted locally. does anyone have a reference that will help me with my problem. i am still new to