similar to: problem with the image command

Dear all, I created a n*2 matrix and I want to replace the ith row with a vector. Can anyone suggest me what is the simple way to do it? Sorry for bothering you with such simple question. I appreciate any hints or advice. Thanks in advance. Cassie [[alternative HTML version deleted]]

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

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?

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

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

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

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'

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

Hi, i'm not english and i'm not very familiar to R, and i'm asking if you can help me. I'm wondering how to create a multivariate normal an then repeat this for a sample of T=1000, and the save this result. Thank you very much for your helping -- View this message in context: http://r.789695.n4.nabble.com/Multivariate-Normal-and-loops-tp4633504.html Sent from the R help mailing

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, when I run the below code,there is an error occured. could you please tell me how to treat it? > pdf('covariate.pdf') > par(mfrow=c(1,1)) > pairs(data2[,c("ID","TYPE","AGE","GNDR","HT")], + panel=function(x,y) { points(x,y); lines(lowess(x,y))}) Error in `[.data.frame`(data2, , c("ID", "TYPE",

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

Hi all, I'm trying to install the mixdist package with: install.packages("mixdist") But I'm getting the following errors: * installing *source* package ‘mixdist’ ... ** R ** data ** preparing package for lazy loading ** help *** installing help indices ** building package indices ... Warning in utils::data(list = f, package = package, lib.loc = lib.loc, envir = dataEnv) :

You probably want to use the merge() function. imaginary example for data frames named "KS" and "US," where the common identifier variable is named "ID" merged.data<-merge(KS,US,by.x="ID",by.y="ID",all.x=F,all.y=F) Note that this will retain only observations for which there is a common ID in both KS and US. If you want to keep observations