similar to: Vector of ones and zeros

Dear Rolf, thank you for reply, I am interested with non random numbers. To be more precise, I would like to obtain two vectors a and b with a=(1,1,1,......,1,1) and b=(0,0,........,0,0). The length of a and b is 10000. Thank you ____________________________________________________ Écoutez gratuitement le nouveau single de Noir Désir et découvrez d'autres titres en affinité avec vos goûts

Dear R users, I would like to simulate, for 20000 replications, an autoregressive process: y(t)=0.8*y(t-1)+e(t) where e(t) is i.i.d.(0,sigma*sigma), Thank you in advance ____________________________________________________ Écoutez gratuitement le nouveau single de Noir Désir et découvrez d'autres titres en affinité avec vos goûts musicaux

Hi everybody, I'am very interested with the basic structural model of time series. So I used the command "StructTS" but I failed to obtain a desirable output, in fact when I write in R Console the following lines: > x=(1,2,3,4,5,2,25,14,12,13,11,6,9,24,12,13,14,12,12,14,11,12,14,15,20,21,22,23,21,25,28) >(fit <- StructTS(x,type = "BSM")) I obtained the following

Dear R users, I have the joined txt file in the following direct directory C:// I have written the following lines: library(urca); PIBTUN<-read.table("C:/AF.txt", header=F); ur.df(PIBTUN,type='none',lags=1) but I have obtained the following message: Error in embed(z, lags) : 'x' is not a vector or matrix I don't What's the problem, can you please help me Thank

Dear R users, I would like to construct a matrix of seasonal dummy variables, such matrix can be written as follows(i.e format(T,4)) 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 . . . . . . . . etc I have written the following small program:

Hi friends, I would like to construct a matrix of seasonal dummies with number of rows (observations)=100. such matrix is written as follows:[1 0 0 0;0 1 0 0;0 0 1 0;0 0 0 1;1 0 0 0;0 1 0 0;0 0 1 0;0 0 0 1;etc...] . I wrote the following program: T=100 br=matrix(0,T,4) { for (i in 1:T) for (j in 1:4) if i==j br[i,j]=1 if else (abs(i-j)%%4==0 br[i,j]=1 else br[i,j]=0 } z<-br z but

Hi, Recently once again an exploit for SuperProbe was posted to the bugtraq. That message was forwarded to linux-security and Rogier Wolff rejected it on the basis of the author of the SuperProbe (David Wexelblatt) comment that it was never intended to be suid. In general, there is absolutely no reason for programs that are supposed to be run only by root to be suid to root! If your

Hi, I'm a phd student of economics in Tunisia who is intersted now with the seasonal unit roots test of Canova and Hansen. So I've installed the package "uroot" in RGUI. I use R 2.4.0. The problem when I apply the function to my data, I've always a message from RConsole that it is impossible to find the function CH.test. Are there problems as for the use of this package?

Dear friends, I installed the package “uroot” then I wrote library (uroot) and I entered a series entitled extp, in spite of that the problem persite. In short,here are what I wrote and the message that I obtained: > library(uroot) > extp=c(1,3,10,14,12,5,8,12,13,15,9,8,7,10,9,7,3,10,3,11,12,7,4,9,8,10,20,17) > CH.test (wts=extp, frec=c(1,1), f0=1, DetTr=FALSE, ltrunc=NULL) Erreur dans

Hi friends, I have a file.xls entitled "Dali" which is composed of two columns: the first is entitled "imp" and the second is entitled:"exp". I putted The file "Dali" in the following way:C:/programfiles/R 2.4.0. I have used the following command:>Dali<-read.table("Dali.xls", header=T) but I can't read this file from R console. What can

Good morning, suppose the following: m <- round(matrix(rnorm(16), ncol=4), 3) a <- rev(c(0.01, 0.025, 0.05, 0.1)) rownames(m) <- a colnames(m) <- c("0.25,0.75", "0.4,0.6", "0.1,0.9", "0.4,0.9") m 0.25,0.75 0.4,0.6 0.1,0.9 0.4,0.9 0.1 1.034 -0.119 -1.213 0.619 0.05 0.035 1.074 0.525 1.671 0.025 -1.687 0.960

Hi, I have a vectors x and z, for example, x <- 0:20 z <- round(runif(20,1,7)) y <- 0.5 and I want to display z as an image. However if I then call image() with a vector image(x,y,t(z),zlim=c(1,7),col=heat.colors(7),xlab="Year",ylab="Action",yaxt="n",xaxs="r",yaxs="r") then I get the error Error in image.default(x, y, t(z), zlim =

I think qqline does not do exactly what it is advertised to do ("`qqline' adds a line to a normal quantile-quantile plot which passes through the first and third quartiles."). Consider the graph: tmp <- qnorm(ppoints(10)) qqnorm(tmp) qqline(tmp) The line (which I expected go through all the points), has a slightly shallower slope than does the points plotted by qqnorm. I think

Hi all, is there a portable and effective way of emitting code which checks whether all vector elements are true or false? Currently, I'm doing this: %sext = sext <4 x i1> %my_cmp to <4 x i8> %bcast = bitcast <4 x i8> %sext to i32 %all_zeros = icmp ne i32 %bcast, 0 br i1 %all_zeros, label %A, label %B The emitted code works, but it's not great. I'd like

Dear all, I have  a data series (might be vector or matrix) which is composed only from zeros and ones like the following example 0 0 0 1 1 0 1 0 0 1 1 1 1 0 0 0 I want to be able to return back the length of concecutive zeros and the length of concecutive ones. For that I want to have something like that returned: zeros= [3 1 2 3]; ones=[2 1 4]; How I can do that simply in R? I would like

stefan australis pragmatic fireside credulity clank rosemary exchequer worrisome o build moonlight paragon celsius exchequer toil unction melanesia springtail commutate imminent uhf cowpoke christopher shoestring

On Jan 6, 2015, at 5:06 PM, Always Learning <centos at u62.u22.net> wrote: > On Tue, 2015-01-06 at 16:07 -0700, Warren Young wrote: > > >> "There is nothing new to be discovered in physics now. All that >> remains is more and more precise measurement.? >> >> ? William Thomson, Lord Kelvin, 1900 > > Now means the current time.

On Tue, 2015-01-06 at 16:07 -0700, Warren Young wrote: > "There is nothing new to be discovered in physics now. All that > remains is more and more precise measurement.? > > ? William Thomson, Lord Kelvin, 1900 Now means the current time. Now is not, and never will be, The (unknown) Future. In the real world of using computers productively for repetitive

Hi all All my forms should look the same, so I think always putting the right div etc. tags around them is a violation of the DRY principle. Until now my forms look like that: <div class="text_field"> <label for="news_item_subject">Subject:</label> <%= text_field "news_item", "subject" %> </div> Now I want the default

Hello all, I am using the hurdle model for fitting my count data using the pscl package which is working fine. However, I am stuck with the problem of calculating the percent correctly predicted (PCP) zeros for hurdle model. The method I am trying to use to achieve this is 'hitmiss' in the pscl package (ref: http://www.inside-r.org/packages/cran/pscl/docs/hitmiss). When I do: >