similar to: replace a row in a matrix

Dear all, I am trying to simulate from truncated Pareto distribution. I know there is a package called PtProcess for Pareto distribution...but it is not for truncated one. Can anyone please help me with this? Thanks in advance. Cassie [[alternative HTML version deleted]]

Dear all, I am trying to simulate from discrete uniform distribution. But I could not find any in-built code in R. Could anyone help me please? Thanks in advance for the time and help. Cassie [[alternative HTML version deleted]]

Hello everyone, I am wondering if there is any routine in R which can convert time given in 'seconds' unit to the 12 hour time format. For example, suppose the data set looks like x=c(36885,84000,20) #x in seconds I want to get the output as [1] " 11:14:45 AM" [2] " 11:20:00 PM" [3] "12:20:00 AM" Does anyone have any idea? Thanks in advance. Cassie

Dear All, I am trying to define a loop for a m*n matrix, where i=1:n and j=1:m. Unlike the usual for loop, i should go in the following way: For j=1, i=1,2,3,....n For j=2, i=n,n-1,n-2,......,1 For j=3, i=1,2,3,.....n etc. which means i should go in either increasing or decreasing order alternatively. Can anyone please help me in doing this? Thanks, Cassie [[alternative HTML version

Hi all! I have a list of matrices and I want to multiply the ith element of the list with the ith element of a another vector. That is, > LL <- list(A=diag(3),B=diag(3),C=diag(3)) > vec <- 1:3 > for(i in 1:3) + { + LL[[i]] <- LL[[i]]*vec[i] + } > LL $A [,1] [,2] [,3] [1,] 1 0 0 [2,] 0 1 0 [3,] 0 0 1 $B [,1] [,2] [,3] [1,] 2 0

Hi all, I have a matrix which looks like the following, x=matrix(c(1,2,3,4,3,5,4,6,4,9),5,2,byrow=T) Now, I want to replace the 2nd columns of those 4 rows which have duplicate entries in the 1st column. For example, x[2,1]=x[3,1]=3 and I want to replace x[3,2] by x[2,2]. Similarly, I want to replace x[5,2] with x[4,2], which means that for duplicate entries in the 1st column, I want to replace

Hi, fellow R users. I have a question about sapply and split combination. I have a big dataframe (40000 observations, 21 variables). First variable (factor) is "date" and it is in format "8.29.97", that is, I have monthly data. Second variable (also factor) has levels 1 to 6 (fractiles 1 to 5 and missing value with code 6). The other 19 variables are numeric. For each month I

Dear all, I am writing a program using for loop. The loop contains i which runs from 1:n. There is a feature I need to include here. That is, when i=n, i+1 would be 1, not n+1. Basically the loop should run in a circular fashion. That also means, if i=1, i-1=n. Can anyone help me with this? How can I define it? Thanks. [[alternative HTML version deleted]]

Dear all, I have a data set which has latitude, longitude (not in increasing order) and observed precipitation corresponding to those locations. I am using image.plot command for that. But I want to overlay a map of US on that plot. I know we can use map('usa') to draw it, but I dont know how can I overlay it on the plot of the precipitation data. Can anyone please help me with it?

Hi everyone, i'm not very good with R, and i know that my question is quite simple. I generated this vector for a sample of 100 and i want to know how to save these as ez1, ez1, ez3...ez100. I guess i have to do it by "for", but i didn't succeed.. i put here my output ezt [1] -0.05470228 -1.76356668 -1.68843463 -0.77283767 0.25077826 -0.65236833 1.05230056 -0.48900254

Dear all, I am trying to define an array as follows: a=array(NA,c(m,2,N,G)), Even though the numbers N and G are fixed, m is a vector. So, I can't define 'a' as above. For eg, N=3 G=6, but m=c(2,7,5) and length(m)=N. Can anyone suggest me how can I define 'a' when m is a vector? Thanks in advance for your help. Cassie [[alternative HTML version deleted]]

Does anyone know how to do a fry plot in R? I have 600-800 points per image, and I really don't want to attempt that manually. Thank you! Cassie [[alternative HTML version deleted]]

I am hoping someone can help translate some WinBUGS code into R code. I would like to use R to create the C[] matrix required for a car.proper model in WinBUGS, but I am having a difficult time negotiating the coding. The C matrix provides normalized weights for each pair of spatial areas. So the WinBUGS example is as follows: # of the weight matrix with elements Cij. The first J1 elements

I am having a hard time putting the below into a loop, where it pulls out ppt from all he stations I have versus having to go through and hard code the data to the specific stations. I tried stnID <- stnid[which(duplicated(stnid)==FALSE)] for(i in 1:length(stnID)) { ppt[i] <- ppt[which(stnid==[i])] } but it doesn't like to use the which function inside a for loop? Any idea's here is

Dear all, I need to simulate values from a Poisson distribution which is truncated at certain value 'a'. Can anyone tell me if there is in-built package in R which can simulate from a truncated Poisson? If not, what should be the steps to write a function which would do that? Thanks in advance. Cassie [[alternative HTML version deleted]]

Hello R-users, I am having a problem with the 'break' command in R. I am wondering if anyone can help me out with this. My program is similar to the following. a=rep(NA,5) a[1]=0 for(i in 2:5) { a[i]=a[i-1]+runif(1,0,3) if(a[i]>5) { i=2 break } } What exactly I am trying to do is find the sequence 'a' where the individual values of 'a'

Hello R-users, I am trying to do an image plot where I have been given latitudes, longitudes and the values at the corresponding locations. A sample of the data is given as follows: values=c(0,1,0,0,0,0,2,2,0,0) lat=c(29.6660,29.6756,29.3492,29.2654,29.2827,29.4070,35.3510,35.6590,35.7587,38.2794)

Dear all, I am trying to plot 3 histograms on the same graph using the following command. hist(x,xlim=c(0,100)) hist(y,add=TRUE) hist(z,add=TRUE) The xlim of y is c(20,21) and that of z is c(99,99.5) , whereas the variable x has xlim at c(0.5,2). Apparently, the graph end of displaying a line only, corresponding to each histogram due to the wide range of the xlim and tightness of the

Hi everyone, I have bunch of date and time observations in the format %Y-%m-%d %I %M %S %p. I used strptime() to read this format. But the problem is some of the times are in the format of %I %M %p, so for those times, strptime is giving me NA values. For example, strptime(paste("2009-04-08","1:49:47 PM"),format="%Y-%m-%d %I:%M:%S %p") [1] "2009-04-08

Hello Folks. This must be a silly question with a (not) obvious (to me) answer. Consider this: tmp <- matrix(1:200, nrow = 20) vec <- 300:309 tmp[9,] <- vec # replacing one row works fine p <- c(3, 11, 17) tmp[p,] <- vec # replacing multple rows pastes the values down a column and recycles vec. What I want to do is replace multiple rows simultaneously at once. I suppose I can