Re: [R] convert annual data to quarterly frequency

On Apr 5, 2013, at 2:12 PM, Ravshonbek Otojanov r.otoja...@qmul.ac.uk wrote: Hi, I was wondering if there is a package I can use to convert my 22 annual observations to quarterly time series so that I will have 88 observations. Please don't post to two lists -- and I think the package

Re: [R] convert annual data to quarterly frequency

On 05-04-2013, at 21:12, Ravshonbek Otojanov r.otoja...@qmul.ac.uk wrote: Hi, I was wondering if there is a package I can use to convert my 22 annual observations to quarterly time series so that I will have 88 observations. Also, will there be any harm in doing so? will I loose any

[R] manipulating R contingency tables

Hi Guys I am back with another thing that's puzzling me. I am creating contingency tables but then I want to filter out certain columns and also find if any entry in the table is 0. Example: gts labels A1 B2 G3 121 127 120 223 112 0 Here I want to remove B2 column

[R] How to plot several years data with date information by months?

Hi, all I have a medium sized data, 6 years. Each observation is a case with a date variable, such as '2004-08-02'. Some of the months didn't occur a case. I want to plot the 6 years data by month, and the Y_axis is the freqency of cases for each month, meaning 12*6=72 bars or points in the

[R] Graphic window dimensions

Dear all, I have a doubt: if I run windows(width=8.27, height=11.69), the size of an A4 paper, does it work correctly in all screens? Or does it depend on the inches of my screen?. I ask you about this question because I need to make the user see a graph in the graph window and he must be

Re: [R] retrieveing value from KS test

thank you. it helps! -- View this message in context: http://r.789695.n4.nabble.com/retrieveing-value-from-KS-test-tp4663439p4663499.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] ad.test parameters

Hi all, i have 2 cumulative (F(x)) distribution function that not defined in R. i want to make anderson darling goodness of fit test to first function (function 1) to check if it behaves as the other distributed function(function 2) how can i insert my 2 unknown by R function to ad.test()?

Re: [R] Text Encoding

Le vendredi 05 avril 2013 à 14:30 -0400, Emily Ottensmeyer a écrit : Dear R-Help, I am using the RDF package/ R 2.14 with the RDF package to download data from a website, and then use R to manipulate it. Text on the website is UTF-8. The RDF package's rdf_load command is converting it

Re: [R] How to plot several years data with date information by months?

Le 06/04/13 10:00, ZhaoXing a écrit : Hi, all I have a medium sized data, 6 years. Each observation is a case with a date variable, such as '2004-08-02'. Some of the months didn't occur a case. I want to plot the 6 years data by month, and the Y_axis is the freqency of cases for each month,

[R] Replace missing value within group with non-missing value

Dear List members I have a large dataset organised in choice groups see sample below +-+ | dn obs choice acid br date cdate situat~n mth year set |

Re: [R] manipulating R contingency tables

Hello, Try the following. The first function removes a column(s) from the table, and the secondd all rows and columns with zero elements in them. fun1 - function(x, col) x[, -which(colnames(x) %in% col)] fun2 - function(x){ idx - which(x == 0, arr.ind = TRUE) x[-idx[, 1],

[R] Replace missing values within a group with the non-missing value

I have a large dataset organised in choice groups see below. Each choice represents a separate occasion with 1 product chosen out of the 6 offered.            +-+      | dn   obs  

[R] Value at Risk using a volatility model?

Hi, I want to calculate the Value at Risk with using some distirbutions and a volatility model. I use the following data(http://uploadeasy.net/upload/cdm3n.rar) which are losses (negative returns) of a company of approx. the last 10 years. So I want to calculated the Value at Risk, this is nothing

Re: [R] Replace missing value within group with non-missing value

Hello, That's not a very good way of posting your data, preferably paste the output of ?dput in a post. Some thing along the lines of the following might do what you want. It seems that the groups are established by 'dn' and 'obs' numbers. If so, try # Make up some data dat - data.frame(dn

Re: [R] How to plot several years data with date information by months?

Le 06/04/13 10:00, ZhaoXing a écrit : Hi, all I have a medium sized data, 6 years. Each observation is a case with a date variable, such as '2004-08-02'. Some of the months didn't occur a case. I want to plot the 6 years data by month, and the Y_axis is the freqency of cases for each month,

Re: [R] Reversing data transformation

Sorry --failed to cc the list. -- Bert On Fri, Apr 5, 2013 at 9:29 PM, Bert Gunter bgun...@gene.com wrote: Have you tried plotting your data?! library(lattice) densityplot(yourdata) Forget tests of normality -- they are a BAD idea, despite what the statistics textbooks say (IMHO obviously)

[R] arrange data

Hello all! I have a problem to arrange data in another form. My initial data is like this: 'data.frame': 421 obs. of 58 variables: $ 01A: num NA NA NA NA NA NA NA NA NA NA ... $ 01B: num NA NA NA NA NA NA NA NA NA NA ... $ 03A: num NA NA NA NA NA NA NA NA NA NA ... $ 03B: num NA NA NA NA

Re: [R] Replace missing value within group with non-missing value

Hi Rui, Thank you for your suggestion which is very much appreciated. Unfortunately running this code produces the following error. error in '$-.data.frame' ('*tmp*', mth, value = NA_real_) : replacement has 1 rows, data has 0 I'm sure there must be an elegant solution to this problem?

Re: [R] Text Encoding

On Apr 5, 2013, at 11:30 AM, Emily Ottensmeyer wrote: Dear R-Help, I am using the RDF package/ R 2.14 with the RDF package to download data from a website, and then use R to manipulate it. Text on the website is UTF-8. The RDF package's rdf_load command is converting it into a

Re: [R] Image segmentation

Hello, Thanks for replying me! I was investigating the EBImage package and actually there is a algorithm for image segmentation. But I think this algorithm is not automated enough for me and maybe to work with multispectral remote sensing image when I want to segment the whole image. The

Re: [R] Replace missing value within group with non-missing value

Hello, Can't you post a data example? If your dataset is named 'dat' use dput(head(dat, 50)) # paste the output of this in a post Rui Barradas Em 06-04-2013 15:34, Leask, Graham escreveu: Hi Rui, Thank you for your suggestion which is very much appreciated. Unfortunately running this

Re: [R] manipulating R contingency tables

Hi, Try this: tbl1- structure(c(21L, 23L, 127L, 112L, 120L, 0L), .Dim = 2:3, .Dimnames = structure(list(     labels = c(1, 2), gts  = c(A1, B2, G3)), .Names = c(labels, gts)), class = table) dat1-as.data.frame(tbl1,stringsAsFactors=FALSE) dat2-dat1[dat1$gts!=B2 dat1$Freq!=0,]

Re: [R] arrange data

On Apr 6, 2013, at 7:28 AM, catalin roibu wrote: Hello all! I have a problem to arrange data in another form. My initial data is like this: 'data.frame': 421 obs. of 58 variables: $ 01A: num NA NA NA NA NA NA NA NA NA NA ... $ 01B: num NA NA NA NA NA NA NA NA NA NA ... $ 03A: num NA

Re: [R] Replace missing value within group with non-missing value

Hi Rui, Data as follows. These data are arranged in choice sets of 6. head(dat,50) dn obs choice br mth 1 4 1 0 1 NA 2 4 1 0 2 NA 3 4 1 0 3 NA 4 4 1 0 4 NA 5 4 1 0 5 NA 6 4 1 1 6 487 7 4 2 0 1 NA 8 4 2 0 2

Re: [R] Replace missing value within group with non-missing value

Hi Rui, Data as follows structure(list(dn = c(4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4), obs = c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5,

[R] Cannot scale data

Hello! I have this error in R: In svm.default(x, y, scale = scale, ..., na.action = na.action) : Variable(s) ‘X.11.’ constant. Cannot scale data. Of course, I have a column of 25 rows where the values are always 0. Why it is not possible to have a column with 0? I need this column because is

Re: [R] Replace missing value within group with non-missing value

On Apr 6, 2013, at 9:16 AM, Leask, Graham wrote: Hi Rui, Data as follows structure(list(dn = c(4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4), obs = c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2,

Re: [R] arrange data

If you want to reorder the column names of your data frame, try (it would have been nice if you provided a subset of your data with 'dput'): # create the order of the new columns indx - order(substring(colnnames(df), 3, 3) , substring(colnames(df), 1, 2) ) df -

Re: [R] Replace missing value within group with non-missing value

Probably not very R-ish but it works (your data in a dataframe called x), if I understand your question right: # replace NA with 0 x$mth - ifelse( is.na( x$mth ), 0, x$mth ) # loop through observation numbers and replace 0 with the month no for( i in unique( x$obs ) ) x$mth[ x$obs == i ] - max(

Re: [R] Replace missing value within group with non-missing value

On Apr 6, 2013, at 9:26 AM, David Winsemius wrote: On Apr 6, 2013, at 9:16 AM, Leask, Graham wrote: Hi Rui, Data as follows structure(list(dn = c(4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,

Re: [R] arrange data

Hi, Try this: #dat1 is dataset indx-apply(dat1,2,function(x) head(which(!is.na(x)),2)) res-as.data.frame(sapply(seq_len(ncol(indx)),function(i) dat2[indx[,i],i]))  colnames(res)- colnames(dat1)  res # 1B    2B    4B    1A    2A    4A    5B    5A  C31A  C31B  C34A  C34B  C35A #1 2.518 2.357

Re: [R] arrange data

 lst1-lapply(dat1,function(x) x[!is.na(x)])  res-as.data.frame(sapply(lst1,function(x) c(x,rep(NA,max(sapply(lst1,length))-length(x)) )))  head(res) # 1B    2B    4B    1A    2A    4A    5B    5A  C31A  C31B  C34A  C34B  C35A #1 2.518 2.357 1.499 3.647 1.890 2.249 2.896 2.175 0.452 1.177

Re: [R] Replace missing value within group with non-missing value

Hi Arun, How odd. Directly pasting the code from your email precisely repeats the error. See below. Any thoughts on the cause of this anomaly? dput(head(dat,50)) structure(list(dn = c(4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4,

Re: [R] Array analogue of row and col

slice.index() in base On 4/2/2013 6:36 AM, Enrico Bibbona wrote: Is there any function that extends to multidimentional arrays the functionalities of row and col which are just defined for matrices? Thanks, Enrico Bibbona [[alternative HTML version deleted]]

Re: [R] Array analogue of row and col

slice.index() in base On 4/2/2013 9:53 AM, Enrico Bibbona wrote: Great! Thanks a lot, Enrico 2013/4/2 Duncan Murdoch murdoch.dun...@gmail.com On 02/04/2013 6:36 AM, Enrico Bibbona wrote: Is there any function that extends to multidimentional arrays the functionalities of row and col which

Re: [R] Replace missing value within group with non-missing value

Hi David, Thank you for your suggestion unfortunately it yields an incorrect result. Within a choice set the mth should be identical. head(dat) dn obs choice br mth 1 4 1 0 1 487 2 4 1 0 2 488 3 4 1 0 3 488 4 4 1 0 4 488 5 4 1 0 5 488 6 4 1

Re: [R] Replace missing value within group with non-missing value

Hello, I've just run my code with your data and found no error. Anyway, try replacing the lapply instruction with this. tmp - lapply(sp, function(x){ idx - which(!is.na(x$mth))[1] if(length(idx) 0) x$mth - x$mth[idx] x

Re: [R] Cannot scale data

Hello, Can't you simply call svm() with scale = FALSE ? If the variable is constant, it cannot be scaled to unit variance (and zero mean). Hope this helps, Rui Barradas Em 06-04-2013 17:19, Nicolás Sánchez escreveu: Hello! I have this error in R: In svm.default(x, y, scale = scale, ...,

Re: [R] Image segmentation

I've recently had a reason to work a little with image segmentation too, and in addition to EBImage, you should look at biOps. You can learn a lot by studying these packages. Bryan On Apr 6, 2013, at 10:04 AM, Eder Paulo eder...@hotmail.com wrote: Hello, Thanks for replying me! I was

[R] Matrix reduction

Hi everyone, I asked myself if exists a way to get a rows (or columns) matrix reduction with R. My goal is for education. Thanks in advance AS __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] Replace missing value within group with non-missing value

Hi Rui, I have just pasted this direct and rerun. I still get the same error. I am running this on the full length dataset. dim [1] 255030 5 I attach the first few rows of the file. Error pasted below. I am using version 2.15.2 on a windows 7 machine. sp - split(dat, list(dat$dn,

Re: [R] Matrix reduction

Incomprehensible. Define:matrix reduction . Perhaps: ?qr ?svd ## and links therein. -- Bert On Sat, Apr 6, 2013 at 10:43 AM, Angelo Scozzarella Tiscali angeloscozzare...@tiscali.it wrote: Hi everyone, I asked myself if exists a way to get a rows (or columns) matrix reduction with R. My

Re: [R] Matrix reduction

Well, I mean to use the elimination method to transform the matrix, for example, of the coefficients of a linear equations system. AS Il giorno 06/apr/2013, alle ore 19:43, Angelo Scozzarella Tiscali ha scritto: Hi everyone, I asked myself if exists a way to get a rows (or columns) matrix

Re: [R] Matrix reduction

On 06-04-2013, at 19:58, Angelo Scozzarella Tiscali angeloscozzare...@tiscali.it wrote: Well, I mean to use the elimination method to transform the matrix, for example, of the coefficients of a linear equations system. AS Well if you only need to solve a linear equation system have a

Re: [R] Replace missing value within group with non-missing value

Hi, dat- read.csv(test1.csv,sep=,,stringsAsFactors=FALSE) sp - split(dat, list(dat$dn, dat$obs)) sp1-sp[lapply(sp,nrow)!=0] #added here names(sp1) - NULL tmp- lapply(sp1,function(x){  idx- which(!is.na(x$mth))[1]  x$mth- x$mth[idx] x  }  ) res- do.call(rbind,tmp)  row.names(res)-1:nrow(res)

Re: [R] Replace missing value within group with non-missing value

Hello, With the attached file, I could reproduce the error but I think the added line does the trick. Rui Barradas Em 06-04-2013 19:20, arun escreveu: Hi, dat- read.csv(test1.csv,sep=,,stringsAsFactors=FALSE) sp - split(dat, list(dat$dn, dat$obs)) sp1-sp[lapply(sp,nrow)!=0] #added here

Re: [R] arrange data

Hi, You could also do this with: dat2- dat1  dat2[]-lapply(dat1,function(x) c(x[!is.na(x)],x[is.na(x)])) row.names(res1)- row.names(dat2)  identical(dat2,res1) #[1] TRUE For your new question: test- data.frame(medrw1)  row.names(test)   [1] 1   2   3   4   5   6   7   8   9   10  11  12  

Re: [R] Replace missing value within group with non-missing value

Hi Arun, Thank you to yourself and Rui. The solution that you provided does indeed work. I greatly appreciate your help. Best wishes Graham -Original Message- From: arun [mailto:smartpink...@yahoo.com] Sent: 06 April 2013 19:20 To: Leask, Graham Cc: R help; Rui Barradas Subject: Re:

Re: [R] manipulating R contingency tables

Hi, You could also do: tbl1[,-match(B2,colnames(tbl1))] #  gts #labels A1  G3  #    1 21 120   #   2 23   0 #or tbl1[,!grepl(B2,colnames(tbl1))] #  gts #labels A1  G3  #    1 21 120   #   2 23   0 If you need to remove columns that contains 0 along with removing a specific column.  

[R] list of distance matrices

Dear UseRs, i want to apply mantel test by comparing a list of 124 distance matrices with a reference distance matrix q. The list of distance matrices was created by the following command.. u-lapply(el, function(x) dist(x)) where el is a data frame of 75 columns and 124 rows. Therefore, the list

Re: [R] Matrix reduction

?chol ##also -- Bert On Sat, Apr 6, 2013 at 11:12 AM, Berend Hasselman b...@xs4all.nl wrote: On 06-04-2013, at 19:58, Angelo Scozzarella Tiscali angeloscozzare...@tiscali.it wrote: Well, I mean to use the elimination method to transform the matrix, for example, of the coefficients of a

Re: [R] Matrix reduction

On 06-04-2013, at 22:01, Bert Gunter gunter.ber...@gene.com wrote: ?chol ##also Quite true and therefore: ?backsolve # forwardsolve and ?qr # qr.solve So there's a lot to choose from. Berend -- Bert On Sat, Apr 6, 2013 at 11:12 AM, Berend Hasselman b...@xs4all.nl wrote:

[R] Fw: Reversing data transformation

From: aguitatie...@hotmail.com Sent: Friday, April 05, 2013 11:47 PM To: r-help@r-project.org ; R Help Subject: Reversing data transformation Hi everybody, I would be very grateful if you could give me your thoughts on the following issue. I need to perform Box-Cox (bcPower€)

[R] Data normalization

Dear all, I’m finding difficulties to normalize this data. Could you provide some input? DATA: c(0.000103113, 0.000102948, 0.000104001, 0.000103794, 0.000104628, 9.2765e-05, 9.4296e-05, 9.5025e-05, 9.4978e-05, 9.8821e-05, 9.7586e-05, 9.6285e-05, 0.00010158, 0.000100919, 0.000103535,

[R] error message sending question to the list

Hi, I tried to send several questions to the lists (both normal R and R-Sig-Finance), but everytime I look them up in the archives my messages end up with the following An embedded and charset-unspecified text was scrubbed... for example see my post here:

Re: [R] [R-SIG-Finance] error message sending question to the list

Post in plain text (not HTML) which is a choice within your email client. And please (!!) don't post the same question to multiple lists. If your question is off-topic, the list membership will redirect as needed. Double posting simply wasted the community's energies by duplicating responses.

Re: [R] list of distance matrices

Hi, Try this: set.seed(25)  el-as.data.frame(matrix(sample(1:50,75*124,replace=TRUE),ncol=75))  lst1-lapply(el,function(x) dist(x))  set.seed(28)  q-as.data.frame(matrix(sample(1:50,75*124,replace=TRUE),ncol=75))  q-dist(q) library(ade4) res- lapply(lst1,function(x) mantel.rtest(q,x,nrepet=)

Re: [R] Data normalization

Hello, See ?scale scale(DATA) # mean 0, unit variance Hope this helps, Rui Barradas Em 06-04-2013 21:21, Beatriz González Domínguez escreveu: Dear all, I’m finding difficulties to normalize this data. Could you provide some input? DATA: c(0.000103113, 0.000102948, 0.000104001,

Re: [R] Replace missing value within group with non-missing value

Anyway, try replacing the lapply instruction with this. tmp - lapply(sp, function(x){ idx - which(!is.na(x$mth))[1] if(length(idx) 0) x$mth - x$mth[idx] x }) Note that which(anyLogicalVector)[1] always has length

Re: [R] categorized complete list of R commands?

Dear Rees Morrison, Re: Franklin, I am very impressed. I ran your code and am amazed at the output. I want to use it in my efforts to figure out which are the most widely used functions, so that I can concentrate on understanding those basics reasonably well. May I ask you two few

Re: [R] Fw: Reversing data transformation

On 06.04.2013 22:21, Beatriz González Domínguez wrote: From: aguitatie...@hotmail.com Sent: Friday, April 05, 2013 11:47 PM To: r-help@r-project.org ; R Help Subject: Reversing data transformation Hi everybody, I would be very grateful if you could give me your thoughts on the following

Re: [R] ad.test parameters

On 06.04.2013 09:45, iritgur wrote: Hi all, i have 2 cumulative (F(x)) distribution function that not defined in R. i want to make anderson darling goodness of fit test to first function (function 1) to check if it behaves as the other distributed function(function 2) You mean you have

Re: [R] Graphic window dimensions

On 06.04.2013 08:53, Eva Prieto Castro wrote: Dear all, I have a doubt: if I run windows(width=8.27, height=11.69), the size of an A4 paper, does it work correctly in all screens? Correctly: yes, but not as you expect. It keeps the aspect ration but in a smaller size if the screen is

Re: [R] Mutliple subsetting of a dataframe based on many conditions

Hi, May be this helps: input1- data.frame(answer=rep(1:4,times=18),p.number=rep(1:18,each=4),session=rep(1:2,each=36),count=rep(1:8,each=9),type=rep(1:3,each=24)) input2-

Re: [R] error message sending question to the list

There is much wisdom in the Posting Guide mentioned in the footer of every email on this list. One pearl is the request to send your email in text format, not HTML. How you do this is specific to your email program, so we cannot help you, but Google is your friend. As to your message, requests

[R] Plotting a curve for a Holling Type III Functional Response

Hey, So I have a scatter plot and I am trying to plot a curve to fit the data based on a Holling Type III functional response. My function is this: nll2-function(a,b) { conefun-(a*DBH^2)/(b^2+DBH^2) nlls2-dnbinom(x=cones ,size=DBH, mu=conefun,log=TRUE) -sum(nlls) } and my plot is this:

[R] Multiple subsetting of a dataframe based on many conditions

Hello Everybody, I'm working with a dataframe that has 18 columns. I would like to subset the data in one of these columns, present, according to combinations of data in six of the other columns within the data frame and then save this into a text file. The columns I would like to use to

Re: [R] Better way of writing R code

Dear Sir, Thanks a lot for your great help. Do appreciate it a lot. In my earlier mail, where I had attached some files, I have realized yesterday that instead of sending the R code customized by me based on your guidance, I had by mistake attached the contents of email. I do apologize to

[R] lmomco - Three-Parameter Pearson 5 Distribution

Dear R forum, I am bit confused and please guide me - (1) Is Pearson Type III Distribution as given in lmomco package same as Three Parameter Pearson 5 Distribution? If not, how do I estimate the parameters of Three Parameter Pearson 5 Distribution? (2) Is there any other R forum dealing

[R] Package ‘FAdist’ - Log-Pearson Type III Distribution

Dear Sir, I am referring to your package FAdist. I wish to know how to estimate the parameters of the distribution - Log-Pearson Type III Distribution? Will it be possible for you to guide me or inform the package in R, I can use to estimate the parameters. Regards Katherine