Re: [R] function for handling time

2011-10-15 Thread Jim Lemon
On 10/16/2011 04:13 AM, Alaios wrote: Dear all I have the following time stamps (in the following format) MeasurementSet$TimeStamps [,1] [,2] [,3] [,4] [,5] [,6] [1,] 201172 13 43 48.718 [2,] 201172 13 43 54.281 [3,] 201172 13 43 59.843 [

[R] Use of ICA for sound

2011-10-15 Thread Noah Silverman
Hi, I'm looking at the "cocktail party" classic problem. I can see how to use ICA to separate the components. But, How do I then create new wav files of the separated sounds so that they can be played? Thanks -- Noah Silverman UCLA Department of Statistics 8208 Math Sciences Building Los Ange

Re: [R] Read wav file into R

2011-10-15 Thread Noah Silverman
Nice. Thank You! -- Noah Silverman UCLA Department of Statistics 8208 Math Sciences Building Los Angeles, CA 90095 On Oct 15, 2011, at 9:41 PM, Luke Miller wrote: > Try the tuneR package. It will read in wav files and has other functions for > manipulating sound data. > On Oct 15, 2011 9:32 PM

[R] right justify right-axis tick values in lattice

2011-10-15 Thread Richard M. Heiberger
How can I right justify the right-axis tick values? They appear in the example below as left-justified. I have tried several different ways and all fail in different ways. The example below creates the right axis tick value with no attempt at adjustment. alternates I have tried are 1. formattin

Re: [R] Read wav file into R

2011-10-15 Thread Luke Miller
Try the tuneR package. It will read in wav files and has other functions for manipulating sound data. On Oct 15, 2011 9:32 PM, "Noah Silverman" wrote: > Hi, > > I'm interested in doing some sound analysis with R. > > Does anyone have any experience/methods for reading in a wav file? > > -- > Noa

[R] Read wav file into R

2011-10-15 Thread Noah Silverman
Hi, I'm interested in doing some sound analysis with R. Does anyone have any experience/methods for reading in a wav file? -- Noah Silverman UCLA Department of Statistics 8208 Math Sciences Building Los Angeles, CA 90095 [[alternative HTML version deleted]] ___

Re: [R] rpois > 0

2011-10-15 Thread Rolf Turner
On 16/10/11 03:45, peter dalgaard wrote: On Oct 15, 2011, at 01:44 , Rolf Turner wrote: On 15/10/11 10:15, knut-o wrote: Hello what is the easiest way to generate rpois(m,lambda) but only values greater than 0 and length = m. tanks, knut The rpospois() function from the VGAM package is what y

Re: [R] glmmadmb help

2011-10-15 Thread chchjames
Thanks for the reply Ben. I tried it with verbose=TRUE, and got about 7 pages of a word doc as an output, that ended with the error "Error in glmmadmb(stainp ~ beetle.ev + Caged * Section/SegmentT + (1 | : The function maximizer failed". I am not sure how I would best go about posting this, as

Re: [R] Turning Data Frame Columns into Vectors

2011-10-15 Thread jawbonemurphy
Hi Daniel, Thanks for the reply, that was the one missing piece! I had been trying to use the tapply() function, but was getting an error message about unequal length until now. I wrote: VAL_mean_xpart <- tapply(X[ ,2], PARTF, mean) and got the column 2 means by factors just as I had hoped. O

Re: [R] function for handling time

2011-10-15 Thread Daniel Malter
I find your question impossible to answer as you do not provide any description of what the matrix columns actually capture and how the variable in your proposed function relate to the columns of this matrix. Best, Daniel alaios wrote: > > Dear all > I have the following time stamps (in the fol

[R] Turning Data Frame Columns into Vectors

2011-10-15 Thread jawbonemurphy
Hi, I converted an Excel file into a .txt file "X.txt" with no header ( http://r.789695.n4.nabble.com/file/n3908157/X.txt X.txt ) and imported it with: X <- read.table("/Users/johnlogandurland/Desktop/X.txt", header=FALSE). What I would like to do is to make the first column into a factors vect

Re: [R] Irregular 3d objects with rgl

2011-10-15 Thread Greg Snow
You could use the rgl package and plot a sprite at each of your points with the color based on the concentration: plume$col <- cut(plume$conc, c(-1,0.01,0.02,0.3,0.7,1), labels=c('blue','green','yellow','orange','red')) plume2 <- plume theta <- atan2(plume2$y-mean(plume2$y), plume2$x-m

Re: [R] How to keep a coefficient fixed when using rq {quantreg}?

2011-10-15 Thread Tal Galili
Thank you for the prompt reply Roger, Your solution works great for me. With much respect, Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il

Re: [R] Turning Data Frame Columns into Vectors

2011-10-15 Thread Daniel Malter
You access columns of a data.frame by column indices as in: X[ ,1], X[ ,2], etc. The index before the comma would stand for the row if you wanted to restrict those. The index after the comma captures the column. That said, you typically would not "extract" rows from the data frame but draw directl

Re: [R] glmmadmb help

2011-10-15 Thread Ben Bolker
chchjames windowslive.com> writes: [snip] > I am using the alpha version of glmmadmb, and it works for most of the time > except for one of my models. The weird thing is that it has worked before, a > couple of months ago, and for some reason it won't now and nothing has > changed. > > The code

[R] function for handling time

2011-10-15 Thread Alaios
Dear all I have the following time stamps (in the following format) MeasurementSet$TimeStamps    [,1] [,2] [,3] [,4] [,5]   [,6]   [1,] 2011    7    2   13   43 48.718   [2,] 2011    7    2   13   43 54.281   [3,] 2011    7    2   13   43 59.843   [4,] 2011    7    2   13   44  5.390   [5,] 20

Re: [R] binomial GLM quasi separation

2011-10-15 Thread lincoln
#Uwe: I have realized that in the firstly linked post ( http://r.789695.n4.nabble.com/OT-quasi-separation-in-a-logistic-GLM-td875726.html#a3850331 OT-quasi-separation-in-a-logistic-GLM ) I have told something misleading: in fact my independent variables are not log-normally distributed since the

[R] [R-pkgs] survey 3.26

2011-10-15 Thread Thomas Lumley
Version 3.26 of the survey package is percolating through CRAN. Since the last announcement on this list, of version 3.20, about 18 months ago, the main changes are -- an option to calculate replicate-weight variances from sums of squares around the point estimate rather than from the varianc

Re: [R] how to plot two surfaces with lattice::wireframe

2011-10-15 Thread Carl Witthoft
Thank you! That bit about "Re(z) + Im(z)" was what I missed. Carl On 10/15/11 12:00 PM, Deepayan Sarkar wrote: On Fri, Oct 14, 2011 at 3:07 AM, Carl Witthoft wrote: Hi all, I'd like to plot the Real and Imaginary parts of some f(z) as two different surfaces in wireframe (the row/column axes

Re: [R] how to plot two surfaces with lattice::wireframe

2011-10-15 Thread Deepayan Sarkar
On Fri, Oct 14, 2011 at 3:07 AM, Carl Witthoft wrote: > Hi all, > I'd like to plot the Real and Imaginary parts of some f(z) as two different > surfaces in wireframe (the row/column axes are the real and imag axes).  I > know I can do it by, roughly speaking, something like > > plotz <- expand.gri

Re: [R] Multiple levelplot with title

2011-10-15 Thread Deepayan Sarkar
On Sun, Oct 9, 2011 at 10:00 PM, Richard O. Legendi wrote: > Hi all, > > I'm new to R and to the mailing list, so please bear with me :-) > > I would like to create multiple levelplots on the same chart with a nice > main title with something like this: > >  print(levelplot(matrix(c(1,2,3,4), 2, 2

Re: [R] Superposing mean line to xyplot

2011-10-15 Thread Deepayan Sarkar
2011/10/10 Niccolò Bassani : > Dear R-users, > I'm using lattice package and function xyplot for the first time so > you will excuse me for my inexperience. I'm facing quite a simple > problem but I'm having troubles on how to solve it, I've read tons of > old mails in the archives and looked at so

Re: [R] "pos" in panel.text

2011-10-15 Thread Deepayan Sarkar
On Mon, Oct 10, 2011 at 6:31 PM, Allan Sikk wrote: > Here's the code. The problem seems to be specific for lattice as I can > easily use a vector with pos in "plot". lattice::panel.text() does not support vector 'pos'. (Not very difficult to fix, and I'll put it on my TODO list). For now, you wil

Re: [R] Legend symbols (line, points) in one column.

2011-10-15 Thread Deepayan Sarkar
On Fri, Oct 14, 2011 at 10:02 AM, David Winsemius wrote: > Legends are built in columns. You need to find a graphics symbol to put in > the "points" column or you need to find something that the lines paramater > will turn into a dot (and I'm not sure what that might be.) A 'lines' component can

Re: [R] rpois > 0

2011-10-15 Thread peter dalgaard
On Oct 15, 2011, at 01:44 , Rolf Turner wrote: > On 15/10/11 10:15, knut-o wrote: >> Hello >> what is the easiest way to generate rpois(m,lambda) but only values greater >> than 0 and length = m. >> tanks, knut > > The rpospois() function from the VGAM package is what you are looking for. > Or

Re: [R] iplots

2011-10-15 Thread Simon Urbanek
On Oct 7, 2011, at 6:32 AM, wrote: > Hi guys, > > I am trying to plot a scatterplot with lines instead of with points using > function iplot in library iplots. > > The best I can come up with is > > x =(1:100) > y = rnorm(100) > iplot(x,y,ptDiam=0) > ilines(x,y) > > Do you guys know of

Re: [R] glmmadmb help

2011-10-15 Thread Uwe Ligges
On 15.10.2011 01:06, chchjames wrote: Hello everyone, I am using the alpha version of glmmadmb, Sounds like you want to contact the author of the code rather than R-help. Uwe Ligges and it works for most of the time except for one of my models. The weird thing is that it has worked befor

Re: [R] Minimum cutsets

2011-10-15 Thread malhomidi
Thanks Gabor for the help. I'll try to do it. Regards, Mohammed -- View this message in context: http://r.789695.n4.nabble.com/Minimum-cutsets-tp885346p3907214.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.

Re: [R] Irregular 3d objects with rgl

2011-10-15 Thread Duncan Murdoch
On 14/10/2011 8:14 PM, emorway wrote: Hello, While exploring if rgl is along the lines of what I need, I checked out demo("rgl") and didn't find quite what I'm looking for, and am therefore seeking additional help/suggestions. The application is geared towards displaying a 3D rendering of a con

Re: [R] strucchange Nyblom-Hansen Test?

2011-10-15 Thread Achim Zeileis
On Fri, 14 Oct 2011, buehlerman wrote: Thanks a lot for your immediate help and detailed explanation! About one thing I'm not quite clear: When the default fit = glm in gefp() is used: sctest(gefp(Employed ~ Year + GNP.deflator + GNP + Armed.Forces, data = longley, fit = lm), functional = mean