similar to: R: method of moments

I made one typo in my previous mail.   May I ask one statistics related question please? I have one query on MLE itself. It's property says that, for large sample size it is normally distributed [i.e. asymptotically normal]. On the other hand it is "Consistent" as well. My doubt is, how this two asymptotic properties exist simultaneously? If it is consistent then asymptotically it

hi all how does one simulate a random walk process? i.e y(0)=0 y(t)=y(t-1)+ e(t) where e(t) is normal(0,1) say. Regards allan

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

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

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

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 want to calculate certain lags of a time series and plot them simultaneously on a graph. can anyone help?

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

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

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

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

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 this is a stats question. when undertaking factor analysis should the variables included always be continiuos random variables? could one include categorical data? e.g some survey type data my gut feeling is that one should not include categorical type data since factor analysis is based on the eigenanalysis of the correlation matrix of a group of variables. The introduction of

hello all i have a quick question. i have been using the RODBC library (trying to read Excel data into R but i am doing this by using Rexcel. this is probably not the correct forum - sorry for this). my code is shown below: Sub A() 'start the connection to R Call RInterface.StartRServer RInterface.RRun "library(RODBC)" RInterface.RRun "A =