[R] fft and filtering puzzle

I am puzzled by a filtering problem using fft(). I don't blame R. I have a waveform y consisting of the sum of 2 sinewaves having freqs f1 and f2. I do s = fft() of y. Remove s's spike at freq=f2 Do inverse fft on s. The resulting waveform still has a lot of f2 in it! But the filtering should

Re: [R] fft and filtering puzzle

I think I've got it! I only filtered out the pos freq and not the negative! (At the far end of s). I will check now. It would be nice to forget about that bit s[n/2:n]. But I get if I get rid of it all, then when I back transform into time domain I imagine I'll get artifacts. Bill On Sat, Nov

[R] import statsoft statistica files?

Is there an option to import .sta files to R? I know you can import SPSS ones. thank you. -- View this message in context: http://n4.nabble.com/import-statsoft-statistica-files-tp875653p875653.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] fft and filtering puzzle

Yep, that was it. Needed to remove neg freq from opposite end of s[] coef- 1-(((1:n)8-1) ((1:n)8+1) | (((1:n)(n-6)-1) ((1:n)(n-6)+1)) ) Sorry to bother you. Cheers Bill __ R-help@r-project.org mailing list

Re: [R] Build of XML package failed

Thaks a lot Duncan, it worked ! Luis Selon Duncan Temple Lang dun...@wald.ucdavis.edu: Hi Luis. You can change the two lines PROBLEM buf WARN; to the one line warning(buf); That should compile. If not, please show us the compilation command for DocParse.c, i.e. all the

Re: [R] using reshape to do ANOVA mixed models

Thank you. On 28/11/2009, at 14:04, Peter Dalgaard p.dalga...@biostat.ku.dk wrote: Or Duek wrote: I'm sorry, I don't know how to copy it to the list. as for the output and functions I've used: Use the Cc:! I am NOT answering list questions in private mails. -p

[R] t-criterion calculation using means and SE

Dear all, is there any functions which allow to calculate Student t-criterion using means, their SE and sample size? I've seek for, but bulit-in t-criterion works only with initial sample... Best regards, A.Morkovin __ R-help@r-project.org mailing

[R] Plot fitted vs observed values

Dear Wiza[R]ds, # I have the following experimentally observed data: csdata - data.frame( time=c(0,1,3,9,20), conc=c(638.697,395.69,199.00,141.58,112.16) ) # weighting resp means response wt.MM- function(resp, time,A1,a1,A2,a2) { pred - A1*exp(-a1*time)+A2*exp(-a2*time) (resp - pred) /

[R] R function that duplicates Octave's poly function?

By any chance is anyone aware of an R function that duplicates Octave's poly function? Here is a description of Octave's poly function: Function File: poly (A) If A is a square N-by-N matrix, `poly (A)' is the row vector of the coefficients of `det (z * eye (N) - a)', the

Re: [R] Plot fitted vs observed values

On 28/11/2009 9:25 AM, oscar linares wrote: Dear Wiza[R]ds, # I have the following experimentally observed data: csdata - data.frame( time=c(0,1,3,9,20), conc=c(638.697,395.69,199.00,141.58,112.16) ) # weighting resp means response wt.MM- function(resp, time,A1,a1,A2,a2) { pred -

[R] Include manually an intercept in lm without breaking it?

Hi Say I want to add manually an intercept in the function lm. Even if almost all results will be identical, few stats are different as DF counting will be different as intercept will not be included in automatic case, while it will be in manual case. See: ###usual lm on freeny

Re: [R] R function that duplicates Octave's poly function?

On Nov 28, 2009, at 9:33 AM, Jason Rupert wrote: By any chance is anyone aware of an R function that duplicates Octave's poly function? Here is a description of Octave's poly function: Function File: poly (A) If A is a square N-by-N matrix, `poly (A)' is the row vector of the

Re: [R] R function that duplicates Octave's poly function?

Cool. Is one of the results the definitive function to use to mimic Octave's poly function? There were about 10 returns, and none clearly identified themselves as the definitive functionality representing the poly capability, so any insight and clarification is appreciated. -

Re: [R] R function that duplicates Octave's poly function?

On Nov 28, 2009, at 10:28 AM, Jason Rupert wrote: Cool. Is one of the results the definitive function to use to mimic Octave's poly function? yes There were about 10 returns, and none clearly identified themselves as the definitive functionality representing the poly capability, so

Re: [R] Include manually an intercept in lm without breaking it?

On 28/11/2009 10:14 AM, Matthieu Stigler wrote: Hi Say I want to add manually an intercept in the function lm. Even if almost all results will be identical, few stats are different as DF counting will be different as intercept will not be included in automatic case, while it will be in

Re: [R] Sampling dataframe

Here are some options that may help you out. First, let's put the data in a format that can be cut-and-pasted into R. myData - read.table(textConnection(var1 var2 var3 1 111 2 312 3 813 4 614 51015 6 221 7 422 8

Re: [R] Reading from Google Docs

Please oh please could someone help me or at least confirm that they are having the same problem. Why am I getting the error message from RGoogleDocs getDocs(sheets.con) Error in getDocs(sheets.con) : problems connecting to get the list of documents How do I troubleshoot? Farrel Buchinsky

[R] Plot-region not complete

Hello there, after I solved many problems in the last days with R I wanted to complete my plot. I used ylab=expression(paste(log(,mu,/,cm^2,/VS))) to have greek symbols and to change cm^2 to cm². Also cm^-3 I have in another plot, but it cuts my ^2 to _ in the y-axis. here is the code and I

[R] Whats happening to the forum???

Whats happening? everthings changing round here...plus alot less posts/users.. -- View this message in context: http://n4.nabble.com/Whats-happening-to-the-forum-tp913858p913858.html Sent from the R help mailing list archive at Nabble.com. __

[R] Exact test for count data

Hello! Bortz, Lienert, Boehnke (2008, pp. 140--142) suggest an exact polynomial test for low frequency tables. I used it recently, and thus, created the code attached. Maybe someone would use (and likely modify) it or incorporate it into their package. Sören References: Bortz, J.,

Re: [R] Whats happening to the forum???

On Nov 28, 2009, at 12:26 PM, DispersionMap wrote: Whats happening? everthings changing round here... Can you be any less vague? plus alot less posts/users.. And perhaps offer some data to accomapny your inferences? http://dir.gmane.org/gmane.comp.lang.r.general -- View this message

Re: [R] Whats happening to the forum???

DispersionMap wrote: Whats happening? everthings changing round here...plus alot less posts/users.. R-Help is a mailing list which is independent of Nabble. Nabble is only one of many front-ends that provides a forum-like interface to the mailing list. Therefore, the users of the R-Help

Re: [R] Whats happening to the forum???

im confused, everytime i try to get to the old forum i get redirected here though??? i dont know why. Sharpie wrote: DispersionMap wrote: Whats happening? everthings changing round here...plus alot less posts/users.. R-Help is a mailing list which is independent of Nabble.

Re: [R] Reading from Google Docs

Farrel Buchinsky wrote: Please oh please could someone help me or at least confirm that they are having the same problem. Why am I getting the error message from RGoogleDocs getDocs(sheets.con) Error in getDocs(sheets.con) : problems connecting to get the list of documents You are

Re: [R] Include manually an intercept in lm without breaking it?

Thanks for your answer! Actually, the reason to include it manually was that the function I wanted to create should have args like: type=c(const, trend, both, none) and so was easier to create a matrix and simply put lm(y~-1+datamat) But following your advices I think rather change that

Re: [R] Export kde object as shapefile

Tyler Dean Rudolph wrote: I am trying to estimate home range size using the plug-in method with kernel density estimation in the kernel smoothing (ks) package. Unless there is another way I am not familiar with, in order to calculate spatial area under the space I need to convert my kde

Re: [R] Whats happening to the forum???

Im lost i want to post on the other forum...the address was something like... http://old.nabble.com/R ...to get help with using R but i kept getting redirected to this one (http://n4.nabble.com), and when i tried to post a querie it made me sign up so i thought the forum was changing or

[R] r-help

Im lost i want to post on the other forum...the address was something like... http://old.nabble.com/R ...to get help with using R but i kept getting redirected to this one (http://n4.nabble.com), and when i tried to post a querie it made me sign up so i thought the forum was changing or being

Re: [R] Whats happening to the forum???

DispersionMap wrote: Im lost i want to post on the other forum...the address was something like... http://old.nabble.com/R ...to get help with using R but i kept getting redirected to this one (http://n4.nabble.com), and when i tried to post a querie it made me sign up so i thought

Re: [R] Plot-region not complete

On Nov 28, 2009, at 12:09 PM, Markus Häge wrote: Hello there, after I solved many problems in the last days with R I wanted to complete my plot. I used ylab=expression(paste(log(,mu,/,cm^2,/VS))) to have greek symbols and to change cm^2 to cm². Also cm^-3 I have in another plot, but it cuts

Re: [R] generating a matrix after a for loop..

Let me see if I understand your question. You want to run your simulation function (mclustSim) several times and collect the resulting vector from each run into a matrix. Is that correct? In that case you can just you the replicate command (as long as your function returns the same length

Re: [R] r-help

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of CHRISTOPHER FRENCH Sent: Saturday, November 28, 2009 10:27 AM To: r-help@r-project.org Subject: [R] r-help Im lost i want to post on the other forum...the address was

Re: [R] Reading from Google Docs

Thank you for the interest in my problem. I have been using the same script (see below) successfully for the past 5 months and now all of a sudden I have problems. Could R be functioning differently under 2.10? Could Google have changed their authentication procedures? In other words are you

[R] Kohonen Package

Hi All, I am still learning R, but making, IMO, great strides. I learned about Kohonen/Self-Organizing Maps in class and I would like to try to replicate some of the things we have seen in class. Below is my code. I am trying to create a u-matrix. In the documentation on page 9 it appears

[R] Orphaned R Packages (maybe this is too inside baseball?)

How do the R powers that be handle packages that are orphaned from CRAN? Recently, I was looking for a function either part of the base functionality or an add-on package that mimicked the poly functionality from Octave

[R] extract coordinates from contour lines

To everyone, Is there a way to extract coordinates or an equation from a straight contour line? Thanks. Cheers Julius Tesoro __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] extract coordinates from contour lines

On Nov 28, 2009, at 4:40 PM, Julius Tesoro wrote: To everyone, Is there a way to extract coordinates or an equation from a straight contour line? Can you set up some sort of example to offer a hint about what you want? -- David Winsemius, MD Heritage Laboratories West Hartford, CT

Re: [R] extract coordinates from contour lines

On 28-Nov-09 21:40:00, Julius Tesoro wrote: To everyone, Is there a way to extract coordinates or an equation from a straight contour line? Thanks. Cheers Julius Tesoro It's not obvious what you mean by a straight contour line! Do you mean one which is/looks like a straight line? If so,

Re: [R] Orphaned R Packages (maybe this is too inside baseball?)

A fairly obvious search strategy (again) does locate the answer but my inclination to offer is greatly diminished by your offlist response to my prior efforts. -- David On Nov 28, 2009, at 4:19 PM, Jason Rupert wrote: How do the R powers that be handle packages that are orphaned from

[R] GMM estimation

Hi there, I'd like to resort generalized methods of moment to estimate a regression model but I can't understand the help file of gmm package. If the regression model is y, the instrumental variable is x, how can I write the code? By the way, I use systemfit's 3sls with method3sls='GMM', but I

[R] 2-character plotting characters?

I am trying to make a plot using the plot command in which I would like the plotting characters to be two-character strings (they're two-letter abbreviations of country names). I've tried the pch argument and this, of course, only produces 1-character strings. Looking through Intro to R and the

Re: [R] 2-character plotting characters?

?text You will have to plot your own text at each point that want. On Sat, Nov 28, 2009 at 6:38 PM, Ben Seligman ben...@gmail.com wrote: I am trying to make a plot using the plot command in which I would like the plotting characters to be two-character strings (they're two-letter

Re: [R] 2-character plotting characters?

Jim, thank you so much; that's exactly what I needed. On Sat, Nov 28, 2009 at 3:42 PM, jim holtman jholt...@gmail.com wrote: ?text You will have to plot your own text at each point that want. On Sat, Nov 28, 2009 at 6:38 PM, Ben Seligman ben...@gmail.com wrote: I am trying to make a plot

[R] R on Large Data Sets (again)

Dear R users, I’ve search the R site for help on this topic but it is hard to find a precise answer for my questions. Which are the best options to overcome the RAM memory limitation problems when using R on “large” data sets (such as 2 or 3 million records)? - Is the free available

Re: [R] Kohonen Package

I haven't used pkg:kohonen, so this is just guessing: 1. are you using an old version of the package? 2. don't call the output of function som() 'som'. As people around here might say: would you call your iguana 'iguana'? 3. have you used str() on your object to check what it really contains?

[R] how to put ggobi display into a GUI window setup by gWidgets

Hi, I want to put a ggobi display into a GUI window setup by gWidgets, but error occur said it is not a S4 object. Does anyone have any idea about how to put it in or maybe it can not be put into a widget at all? Thanks A LOT! -- View this message in context:

Re: [R] R on Large Data Sets (again)

Hello Lars, On 2009.11.28 18:53:09, Lars Bishop wrote: Dear R users, I?ve search the R site for help on this topic but it is hard to find a precise answer for my questions. Which are the best options to overcome the RAM memory limitation problems when using R on ?large? data sets (such

[R] How to find where the source code of an R function or package is installed?

I'm wondering where is the source of an R function or a package is. For example, where is 'attributes'? attributes function (obj) .Primitive(attributes) I also do understand what .Primitive mean. Could somebody let me know how to locate source file in an R installation? Why typing 'attributes'

Re: [R] extract coordinates from contour lines

Hi I just found the answer. I used contourLines() to extract the z values. Thanks anyway. --- On Sun, 11/29/09, David Winsemius dwinsem...@comcast.net wrote: From: David Winsemius dwinsem...@comcast.net Subject: Re: [R] extract coordinates from contour lines To: Julius Tesoro

Re: [R] What is the equivalent of column.prods() from S in R?

On Mon, Nov 9, 2009 at 8:35 PM, David Winsemius dwinsem...@comcast.net wrote: On Nov 9, 2009, at 8:43 PM, Peng Yu wrote: Chambers' book Statistical Models in S mentioned 'column.prods()'. But I don't find it in R. I'm wondering if there is an equivalent in R? ??rowProds ??colProds (They

Re: [R] How to find where the source code of an R function or package is installed?

Peng Yu wrote: I'm wondering where is the source of an R function or a package is. For example, where is 'attributes'? attributes function (obj) .Primitive(attributes) I also do understand what .Primitive mean. Could somebody let me know how to locate source file in an R

Re: [R] How to find where the source code of an R function or package is installed?

On Nov 28, 2009, at 11:00 PM, Peng Yu wrote: I'm wondering where is the source of an R function or a package is. For example, where is 'attributes'? attributes function (obj) .Primitive(attributes) I also do understand what .Primitive mean. Looking below, it appears your understanding

Re: [R] What is the equivalent of column.prods() from S in R?

On Nov 28, 2009, at 11:28 PM, Peng Yu wrote: On Mon, Nov 9, 2009 at 8:35 PM, David Winsemius dwinsem...@comcast.net wrote: On Nov 9, 2009, at 8:43 PM, Peng Yu wrote: Chambers' book Statistical Models in S mentioned 'column.prods()'. But I don't find it in R. I'm wondering if there is an

[R] Plotting observed vs. fitted values

Dear Wiza[R]ds, I am very grateful to Duncan Murdoch for his assistance with this problem. His help was invaluable. However, the problem has become a little more complicated for me. Now, in each plot, I need to plot the observed and fitted values of a supine and upright posture experiment. Here

Re: [R] What is the equivalent of column.prods() from S in R?

On Sat, Nov 28, 2009 at 11:00 PM, David Winsemius dwinsem...@comcast.net wrote: On Nov 28, 2009, at 11:28 PM, Peng Yu wrote: On Mon, Nov 9, 2009 at 8:35 PM, David Winsemius dwinsem...@comcast.net wrote: On Nov 9, 2009, at 8:43 PM, Peng Yu wrote: Chambers' book Statistical Models in S

Re: [R] What is the equivalent of column.prods() from S in R?

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Peng Yu Sent: Saturday, November 28, 2009 8:29 PM To: r-h...@stat.math.ethz.ch Subject: Re: [R] What is the equivalent of column.prods() from S in R? On Mon, Nov 9, 2009 at

[R] Convergence problem with zeroinfl() and hurdle() when interaction term added

Hello, I have a data frame with 1425 observations, 539 of which are zeros. I am trying to fit the following ZINB: f3-formula(Nbr_Abs~ Zone * Year + Source) ZINB2-zeroinfl(f3, dist=negbin, link= logit, data=TheData, offset=log(trans.area), trace=TRUE) Zone is a factor with 4 levels, Year a

Re: [R] What is the equivalent of column.prods() from S in R?

On Nov 29, 2009, at 12:20 AM, Peng Yu wrote: On Sat, Nov 28, 2009 at 11:00 PM, David Winsemius dwinsem...@comcast.net wrote: On Nov 28, 2009, at 11:28 PM, Peng Yu wrote: On Mon, Nov 9, 2009 at 8:35 PM, David Winsemius dwinsem...@comcast.net wrote: On Nov 9, 2009, at 8:43 PM, Peng Yu

Re: [R] R on Large Data Sets (again)

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Jason Morgan Sent: Saturday, November 28, 2009 6:18 PM To: Lars Bishop Cc: r-help@r-project.org Subject: Re: [R] R on Large Data Sets (again) Hello Lars, On 2009.11.28

Re: [R] How to find where the source code of an R function or package is installed?

Check out: Uwe Ligges. R Help Desk: Accessing the sources. R News, 6(4):43-45, October 2006 On Sat, Nov 28, 2009 at 11:00 PM, Peng Yu pengyu...@gmail.com wrote: I'm wondering where is the source of an R function or a package is. For example, where is 'attributes'? attributes function (obj)