Re: [R] pulling recessions out of a hydrograph

2018-02-05 Thread stephen sefick
I have done this. I would use zoo for the time series part, and I think I wrote a function using the min max after a period to locate the peak, go to the peak, and calculate the slope 4 hours afterward. If I can locate the code I don't mind sharing, but it is buried on my PhD machine. On Feb 5, 20

Re: [R] Include pre-existing PDF files as vignettes in an R package?

2018-02-19 Thread stephen sefick
Yes, 'C-c C-e l l' I think, but follow the pop-up, and everything should be ok. On Feb 19, 2018 07:05, "Rainer Krug" wrote: > > > > On 19 Feb 2018, at 12:25, Duncan Murdoch > wrote: > > > > On 19/02/2018 5:47 AM, Michael Hannon wrote: > >> Thanks, Duncan. The f

[R] system call removes special characters from text output

2017-04-03 Thread stephen sefick
Hello, I am writing an R package, and I am using system() to call a perl script. The output of the perl script is correct except for "[A/B]" is output as "AB". Can someone explain this behavior. I would like to try and fix this. many thanks, Stephen Sefick -- Let'

Re: [R] system call removes special characters from text output

2017-04-03 Thread stephen sefick
ty. > > On April 3, 2017 8:57:07 PM PDT, stephen sefick wrote: > >Hello, > > > >I am writing an R package, and I am using system() to call a perl > >script. > >The output of the perl script is correct except for "[A/B]" is output > >as > >&qu

Re: [R] system call removes special characters from text output

2017-04-04 Thread stephen sefick
call to "echo". > Then state what you expect to happen and what happens instead. > > B. > > > > > On Apr 4, 2017, at 12:44 AM, stephen sefick wrote: > > > > Hi Jeff, > > > > My apologies for not providing enough information. The perl code

Re: [R] system call removes special characters from text output

2017-04-05 Thread stephen sefick
; > At a shell prompt: > > [72]% ./tmp.pl > [A/B] > > Inside R: > > > system(' ./tmp.pl') > [A/B] > > > -- > Don MacQueen > > Lawrence Livermore National Laboratory > 7000 East Ave., L-627 > Livermore, CA 94550 > 925-423-1062

[R] colorspace namespace problem with R CMD check --as-cran

2017-04-11 Thread stephen sefick
object ‘plot.index’ ** help *** installing help indices ** building package indices ** testing if installed package can be loaded Thank you for any help in advance. kindest regards, Stephen Sefick -- Let's not spend our time and resources thinking about things that are so little or so large

Re: [R] colorspace namespace problem with R CMD check --as-cran

2017-04-11 Thread stephen sefick
olorspace? Such a function would > have the environment getNamespace("colorspace"). > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > > On Tue, Apr 11, 2017 at 2:47 PM, stephen sefick wrote: > > I am at a loss. I do not know how to make this reproducible w

Re: [R] colorspace namespace problem with R CMD check --as-cran

2017-04-11 Thread stephen sefick
> have the environment getNamespace("colorspace"). > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > > On Tue, Apr 11, 2017 at 2:47 PM, stephen sefick wrote: > > I am at a loss. I do not know how to make this reproducible without > > providing t

Re: [R] colorspace namespace problem with R CMD check --as-cran

2017-04-13 Thread stephen sefick
-dependencies-td918313.html#a918314 ). Thanks for the help. kindest regards, Stephen Sefick On Tue, Apr 11, 2017 at 6:26 PM, stephen sefick wrote: > Here is the relevant snippett from devtools::check() logs. I searched > every file in pkg/R, and there is no call to plot.index (or plot

[R] S4 Question; show generic method

2017-05-18 Thread stephen sefick
Hello, I changed the name of a function, and updated my R packages with update.packages(ask=FALSE), and I have had a curious (to me) change in show method. I am developing a package, and a generic show method I defined for my S4 class is not autocalled when I type the object name - object_name; ho

[R] Rcpp on R 3.2.4 compile problem Linux

2016-06-21 Thread stephen sefick
grDevices utils datasets methods base loaded via a namespace (and not attached): [1] tools_3.2.4 kindest regards, -- Stephen Sefick, PhD ** Auburn University Biological Sciences

Re: [R] changing x and y ranges in a PCA plot created by library(labdsv)

2016-08-07 Thread stephen sefick
Can you provide code and data (with dput)? You will likely get an answer more quickly. On Sun, Aug 7, 2016 at 4:34 AM, Mohsen Sharafatmandrad < mohsen.sharafatmand...@gmail.com> wrote: > Hi, > > I want to change x and y ranges in a PCA plot created by library(labdsv). > When I run "plot(o.pca, xl

[R] digital online identifier (doi) CRAN manual?

2016-09-09 Thread stephen sefick
Hello all, I apologize if this is the incorrect forum for this query. I am a package maintainer, and would like to have a doi for the package manual. What might be a good way to go about this? kindest regards, Stephen Sefick -- Let's not spend our time and resources thinking about things

Re: [R] pls - No variable selection

2015-05-27 Thread stephen sefick
= > > Barry King > Associate Professor of Information Technology > Butler University > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, se

Re: [R] problems editing R console

2015-06-07 Thread stephen sefick
> >> "Many admire, few know" > >> Hippocrates > >> > >> __ > >> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > >> https://stat.ethz.ch/mailman

Re: [R] How to validate the cluster analysis?

2015-06-10 Thread stephen sefick
You could look at silhouette width in, I think, the cluster package. Please excuse my brevity; this message was sent from my telephone. On Jun 10, 2015 1:37 AM, "My List" wrote: > All, > > I am new to the world of statistics. I am interested in finding out the > validation techniques employed on

[R] RcppArmadillo compilation errors (Scientific Linux 6.5)

2014-12-03 Thread stephen sefick
_MONETARY=en_US.utf8LC_MESSAGES=en_US.utf8 [7] LC_PAPER=en_US.utf8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_US.utf8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base

Re: [R] RcppArmadillo compilation errors (Scientific Linux 6.5)

2014-12-03 Thread stephen sefick
solved; sorry for the spam. library(devtools) install_github("RcppCore/RcppArmadillo") On Wed, Dec 3, 2014 at 4:51 PM, stephen sefick wrote: > I would appreciate any help that you may be able to give. Please let me > know if any more information is required. > > I get

[R] rJava Scientific Linux 6.5 - Can not install

2015-02-02 Thread stephen sefick
a namespace (and not attached): [1] colorspace_1.2-4 digest_0.6.8 grid_3.1.0 gtable_0.1.2 [5] MASS_7.3-37 munsell_0.4.2plyr_1.8.1 proto_0.3-10 [9] Rcpp_0.11.4 reshape2_1.4.1 scales_0.2.4 stringr_0.6.2 [13] tools_3.1.0

Re: [R] rJava Scientific Linux 6.5 - Can not install

2015-02-02 Thread stephen sefick
OO#.. Dead: OO#.. Playing > Research Engineer (Solar/BatteriesO.O#. #.O#. with > /Software/Embedded Controllers) .OO#. .OO#. rocks...1k > ------- > Sent from my phone. Please excuse

Re: [R] rJava Scientific Linux 6.5 - Can not install

2015-02-03 Thread stephen sefick
JAVA_HOME="/usr/lib/jvm/java-1.7.0-openjdk-1.7.0.75.x86_64/" 5) Ensured all of the fields in the output of 4 were populated 6) in a sudo R session, I ran install.packages("rJava", library=personal_library) Stephen On Mon, Feb 2, 2015 at 2:39 PM, stephen sefick wrote: > Incl

Re: [R] Installing R on Linux Red Hat Server

2015-03-12 Thread stephen sefick
[[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, mini

Re: [R] Multiple Plots using ggplot

2015-03-30 Thread stephen sefick
; __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-cont

Re: [R] Multiple Plots using ggplot

2015-03-31 Thread stephen sefick
> > 100 125 174 >> >5 1956-01-01 123 1266.44492 >> > 92 119 119 >> >6 1957-01-01 124 1235.96492 >> > 92

Re: [R] Multiple Plots using ggplot

2015-03-31 Thread stephen sefick
uot;) > > but the error is: > > Don't know how to automatically pick scale for object of type function. > Defaulting to continuousError in data.frame(colour = function (x, ...) : > arguments imply differing number of rows: 0, 183 > > > Thanks, > > Frederic. > > >

Re: [R] PCA analysis and bootstraped loadings

2015-04-13 Thread stephen sefick
ppreciated > Cheers! > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/ > posting-guide.html > and prov

[R] R Foundation publish citable Manual for packages

2014-08-06 Thread stephen sefick
sure that everyone that makes the R ecosystem so useful gets attribution for their work. Many thanks to everyone involved. -- Stephen Sefick ** Auburn University Biological Sciences 331 Funchess Hall Auburn, Alabama 36849

Re: [R] Choosing columns by number

2015-08-25 Thread stephen sefick
_ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self

Re: [R] Fwd: ayuda- help me

2015-08-25 Thread stephen sefick
You will generally get more and better responses if you have a reproducible example. Please read the posting guide. Please excuse my brevity; this message was sent from my telephone. On Aug 25, 2015 4:14 PM, "Paulo Fonseca" wrote: > Cordial saludo, > > Muy amablemente solicito su ayuda con la si

Re: [R] reshape: melt and cast

2015-08-31 Thread stephen sefick
This is very hard to read. Please use dput to provide data. I believe the answer is in the manual. Look at the aggregation function argument. Please excuse my brevity; this message was sent from my telephone. On Sep 1, 2015 12:11 AM, "Matt Pickard" wrote: > Hi, > > I have data that looks like th

Re: [R] reshape: melt and cast

2015-09-01 Thread stephen sefick
t; "1556", "1558", "1559"), class = "factor"), RaterName = > structure(c(2L, > 2L, 2L, 2L, 2L, 2L), .Label = c("cwormhoudt", "zspeidel"), class = > "factor"), > SI1 = c(1L, 1L, 1L, 1L, 1L, 1L), SI2 = c(3L, 2L

Re: [R] reshape: melt and cast

2015-09-01 Thread stephen sefick
You are welcome. On Tue, Sep 1, 2015 at 10:44 AM, Matthew Pickard < matthew.david.pick...@gmail.com> wrote: > Yep, that works. Thanks, Stephen. I should have drawn the parallel with > Excel Pivot tables sooner. > > On Tue, Sep 1, 2015 at 9:36 AM, stephen sefick wrote: >

Re: [R] Quantile Regression without intercept

2015-10-05 Thread stephen sefick
t; and provide commented, minimal, self-contained, reproducible code. > -- Stephen Sefick ** Auburn University Biological Sciences 331 Funchess Hall Auburn, Alabama 36849 **

Re: [R] Installing pre-compiled R in Linux

2015-10-07 Thread stephen sefick
om > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented,

Re: [R] monte carlo simulations in permanova in vegan package

2015-10-27 Thread stephen sefick
gt; [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project

Re: [R] monte carlo simulations in permanova in vegan package

2015-10-27 Thread stephen sefick
easonable number of permutations. Hopefully > > this clarifies my situation and aim? > > > > > > > > I was therefore hoping there was a way of coding for the Monte-Carlo > > permutation procedure into adonis? > > > > > > > > Thanks for your h

Re: [R] Alternatives for explicit for() loops

2015-11-06 Thread stephen sefick
ppp<- sum(LED[j,])+(m-1) > >>> hlm[i]<- > >>> (choose(ppp,i))*((-1)^(i))*((i+1)^((-1)*((1/beta)+1))) > >>> } > >>> shl<-gil*(sum(hlm)+1) > >>> retu

Re: [R] FD package

2016-02-22 Thread stephen sefick
ender insists on that the respective contract is concluded only > > upon an express mutual agreement on all its aspects. > > - the sender of this e-mail informs that he/she is not authorized to > enter > > into any contracts on behalf of the company except for cases in which > > h

Re: [R] FD package

2016-02-22 Thread stephen sefick
Marine, brackish, freshwater 3.381153 Invertebrate Feeder >Browser > trait6 > com1 0.5 - 1.0 > com2 1.0 - 2.0 > com3 1.0 - 2.0 > com4 1.0 - 2.0 > com5 1.0 - 2.0 > com6 1.0 - 2.0 > com7 1.0 - 2.0 > com8 1.0 - 2.0 > com9 0.5 - 1.0 > com1

Re: [R] Melt and compute Max, Mean, Min

2016-11-18 Thread Stephen Sefick
Milu, I am unsure what you are trying to do. Please provide a minimal reproducible example, and someone will likely be able to help you with the R code. What have you tried so far? kindest regards, Stephen On 11/18/2016 06:49 AM, Miluji Sb wrote: Dear all, I have 51 years of data (1960 - 2

[R] Make sure a data frame has been "fun through" a function

2017-02-20 Thread stephen sefick
Hello, I would like to add something to a data frame that is 1) invisible to the user, 2) has no side effects, and 3) I can test for in a following function. Is this possible? I am exploring classes and attributes and I have thought about using a list (but 1 and 2 not satisfied). Any help would be

Re: [R] Make sure a data frame has been "fun through" a function

2017-02-20 Thread stephen sefick
user" then > yes, e.g., add something to the class attribute. If 1) means "can't be > discovered by the user" then no (at least not easily). Anything you can see > they can see. > > Best, > Ista > > > On Feb 20, 2017 4:21 PM, "stephen sefick" wrote:

Re: [R] Make sure a data frame has been "fun through" a function

2017-02-20 Thread stephen sefick
; Cheers, > Bert > > > Bert > > On Feb 20, 2017 1:27 PM, "Ista Zahn" wrote: > > It depends on what you mean by 1). If you mean "won't annoy the user" then > yes, e.g., add something to the class attribute. If 1) means "can't be > discover

Re: [R] Make sure a data frame has been "fun through" a function

2017-02-20 Thread stephen sefick
PM, Charles C. Berry wrote: > On Mon, 20 Feb 2017, stephen sefick wrote: > > Hello, >> >> I would like to add something to a data frame that is 1) invisible to the >> user, 2) has no side effects, and 3) I can test for in a following >> function. Is this possible

Re: [R] Make sure a data frame has been "fun through" a function

2017-02-21 Thread stephen sefick
uot;, "Alt") coded_genotype <- alleles[alleles!=impossible_genotype] return(coded_genotype) } ##stop if not pre-processed with functionA functionB(test_data) ##processed with functionA functionB(test_data_genotypes) On Tue, Feb 21, 2017 at 6:41 AM, David Winsemius wrote:

Re: [R] Make sure a data frame has been "fun through" a function

2017-02-24 Thread stephen sefick
> > wasChecked(d) > [1] TRUE > > I don't know if this would be an issue in your case. If it is, you > could subclass "data.frame" and define methods so that the operations > of interest preserve or remove the attribute in the way that you > desire. > > B

Re: [R] code reduction (if anyone feels like it)

2008-07-16 Thread stephen sefick
ot;)) > .time <- data.frame( >sunriset(.sp, study_seq, direction="sunrise", > POSIXct.out=TRUE)$time, >sunriset(.sp, study_seq, direction="sunset", POSIXct.out=TRUE)$time) > colnames(.time) <- paste(.name, c('sr', 'ss

[R] finding "chuncks" of data that have no NA s

2008-07-18 Thread stephen sefick
I have a data frame that is 122 columns and 7ish rows it is a zoo object, but could be easily converted or read in as something else. It is multiparameter multistation water quality data - there are a lot of NA s. I would like to find "chuncks" of data that are free of NA s to do some analysi

Re: [R] finding "chuncks" of data that have no NA s

2008-07-18 Thread stephen sefick
of non-NA data. > Its a zoo method of the na.contiguous generic in the core of R. > > rle(!is.na(rowSums(coredata(z will find all stretches. > > On Fri, Jul 18, 2008 at 6:47 AM, stephen sefick <[EMAIL PROTECTED]> wrote: > > I have a data frame that is 122 columns an

Re: [R] finding "chuncks" of data that have no NA s

2008-07-18 Thread stephen sefick
> > TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, > > TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, > > TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, > > TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, TRUE, FALSE, TRUE, FALSE

[R] source code functions

2008-07-18 Thread stephen sefick
>na.contiguous function (object, ...) UseMethod("na.contiguous") this is what I get when I look for the source code for some functions- Is there a way to look at the source code? -- Let's not spend our time and resources thinking about things that are so little or so large that all they really

[R] extracting colnames to label plots in a function

2008-07-19 Thread stephen sefick
#this is my little function that I would like to use the column names of the x and y arguments in the function. I would like it to read #site1-site2 how would I do this diff.temp <- function(x, y ,use="pairwise.complete.obs") { na.method <- pmatch(use, c("all.obs", "complete.obs",

Re: [R] extracting colnames to label plots in a function

2008-07-19 Thread stephen sefick
lnames(x) d <- (x[,a]-x[,b]) plot(d, main= paste(f[a], "-", f[b])) plot(density(na.omit(coredata(d))), main="") } On Sat, Jul 19, 2008 at 12:04 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > #this is my little function that I would like to use the

Re: [R] extracting colnames to label plots in a function

2008-07-19 Thread stephen sefick
hanged afterwards. See ?par. It could also benefit from the use of > ?match.arg > > Also please include drivers that call the posted function so one can > run them in > a reproducible manner. > > > On Sat, Jul 19, 2008 at 12:04 PM, stephen sefick <[EMAIL PROTECTED]> >

Re: [R] extracting colnames to label plots in a function

2008-07-19 Thread stephen sefick
quot;, "pairwise.complete.obs")) opar <- par par(mfrow=c(2,1)) f <- colnames(x) d <- (x[,a]-x[,b]) plot(d, main= paste(f[a], "-", f[b])) plot(density(na.omit(coredata(d))), main="") par(opar) } give a NULL after execution what i

Re: [R] saving plot both as jpg and pdf

2008-07-22 Thread stephen sefick
you have to go in sequence- you are opening a pdf device then a jpeg device and ploting to the jpeg device. My suggestion would be to do pdf() plot(yourplot) dev.off() jpeg plot(yourplot) dev.off() this should work On Tue, Jul 22, 2008 at 9:50 AM, Monica Pisica <[EMAIL PROTECTED]> wrote: > > H

[R] rollmean and stl

2008-07-22 Thread stephen sefick
I need to investigate how rollmean and the trend returned from stl differ. I am trying to find out exactly what the trend part of stl is (I have just started coding in R and do not know fortran). I need to extract this because it will be used in further calculations, and it needs to be verified to

[R] .asp (not really an R question) but related

2008-07-22 Thread stephen sefick
I am trying to get the paper that explains the stl procedure it is a .asp file. What can I use to read this. thanks stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We

Re: [R] Editor for Mac

2008-07-23 Thread stephen sefick
simultron On Wed, Jul 23, 2008 at 4:04 AM, Hans-Peter Suter <[EMAIL PROTECTED]> wrote: > > Thanks for the answers, > > but I'm looking for an editor like Tinn-R. > > What about VM Fusion and continue to use Tinn-R then? > > -- > Regards, > Hans-Peter > > PS: I use TextMate (and like it) > > _

Re: [R] Just 2 more questions - for now!

2008-07-24 Thread stephen sefick
lm(x~y+z) this should work- I think On Thu, Jul 24, 2008 at 7:01 AM, Williams, Robin < [EMAIL PROTECTED]> wrote: > Hi all, > Thanks for the help with my previous post. > I have just two more questions for the minute. > I think I said in a previous post that I like to use the terminal, > i.e. run

Re: [R] Question: how to build a subset to do separate calculations

2008-07-24 Thread stephen sefick
or sapply as I have recently found to be useful. On Thu, Jul 24, 2008 at 4:33 PM, Erik Iverson <[EMAIL PROTECTED]> wrote: > Spencer - > > Spencer wrote: > >> >> Dear R Experts, >> >> I am trying to create several subsets that I want to use as separate >> datasets. After separate subsets are creat

Re: [R] Question: how to build a subset to do separate calculations

2008-07-24 Thread stephen sefick
sorry disreguard On Thu, Jul 24, 2008 at 4:45 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > or sapply as I have recently found to be useful. > > > On Thu, Jul 24, 2008 at 4:33 PM, Erik Iverson <[EMAIL PROTECTED]> > wrote: > >> Spencer - >> >> Spe

Re: [R] Should this PDF render correctly without font embedding?

2008-07-24 Thread stephen sefick
Mac OS X 10.3.5 R 2.7.1 works fine in the quartz window. Screwed up in Preview with the pdf device On Thu, Jul 24, 2008 at 4:29 PM, Jeffrey Horner <[EMAIL PROTECTED]> wrote: > A professor here at Vanderbilt sent me the following code. Each of the text > strings should right justify against th

Re: [R] Rolling range and regression calculations

2008-07-25 Thread stephen sefick
how about rollapply in the zoo package? On Fri, Jul 25, 2008 at 8:37 AM, rcoder <[EMAIL PROTECTED]> wrote: > > Hi Achim, > > Thanks for your reply. rollmean and rollmax functions exist, but is there > anything for returning the minima on a rolling basis? I know there is no > rollmin in the zoo li

[R] Scatterplot matrix one column vs. the rest of the dataset

2008-07-25 Thread stephen sefick
have a data frame composed of many columns. I would like to hold one column constant and scatter plot it agianst all of the other columns? how do I do this? I have only been able to plot(x) and use indexes to create an entire scatterplot matrix, but not hold one constant. -- Let's not spend ou

Re: [R] Scatterplot matrix one column vs. the rest of the dataset

2008-07-25 Thread stephen sefick
I figured it out- always include the numbers of the columns in the plot statement that you want to scatter. On Fri, Jul 25, 2008 at 9:23 AM, stephen sefick <[EMAIL PROTECTED]> wrote: > have a data frame composed of many columns. I would like to hold one > column constant and sca

[R] "Multiple" Regressions indexed by goodness of fit

2008-07-25 Thread stephen sefick
I would like do batch linear regression and then only have the results with high R-squared's returned. Is there such a data mining package in R. thanks stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff u

[R] plot zoo custom panel help

2008-07-25 Thread stephen sefick
#the below code is the way that I would like the plot to look. I have tried to write a panel function: my.panel <- function(x, y, ..., pf = parent.frame()) { axis(side=1, at = seq(rng[1], rng[2], 1/12), labels = n, tcl = -0.3) } #but it does not work and I am at a loss and help would be apprec

Re: [R] plot zoo custom panel help

2008-07-25 Thread stephen sefick
I am still having the can't draw a plot without putting in a ylim value when trying to draw a plot. I have included reproducible code sorry for the previous post > #the below code is the way that I would like the plot to look. I have > tried > to write a panel function: > my.panel <- function(x

Re: [R] Rolling range and regression calculations

2008-07-25 Thread stephen sefick
ke to find some way of generating a matrix > containing rows that are the product of sequential rows in a preceeding > matrix: i.e. rows 1*2 -> row 1 in o/p matrix; rows2*3-> row 2; rows > 4*5->row3 etc. > > Thanks, > > rcoder > > > > stephen sefick wrote

Re: [R] create multi windows plot at the same time

2008-07-25 Thread stephen sefick
or quartz() or windows() you could also use par(mfrow=c(2,2)) On Fri, Jul 25, 2008 at 3:05 PM, <[EMAIL PROTECTED]> wrote: > Hi, you only need to add x11() in front of plot(); > such as > x11() > Plot(v1) > x11() > Plot(v2) > .. > then R will open multi windows for you > > Chunhao > > > Quotin

Re: [R] plot zoo custom panel help

2008-07-25 Thread stephen sefick
to > reposting. > > On Fri, Jul 25, 2008 at 1:37 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > > > > I am still having the can't draw a plot without putting in a ylim value > when > > trying to draw a plot. I have included reproducible code sorry for th

Re: [R] Numerical question

2008-07-25 Thread stephen sefick
If it is in a data frame with the variables as columns try f <- cbind(a*b, b*c, c*d) sapply(f, function(x)var(x)) On Fri, Jul 25, 2008 at 5:06 PM, Zhang Yanwei - Princeton-MRAm < [EMAIL PROTECTED]> wrote: > Hi all, > I have n independent variables A_1, A_2, A_3,..,A_n, and each with > known

[R] Interpolating a line and then summing there values for a diurnal oxygen curve (zoo object)

2008-07-27 Thread stephen sefick
#I would like to interpolate a straight line between 06/08/06 04:16:00 - 06/08/06 20:31:00 with values and then sum them. This is an estimate of ecosystem #respiration and I will be using this in a larger context(48 days of these diurnal curves), but for right now I am just trying to figure out h

Re: [R] writing the plots

2008-07-28 Thread stephen sefick
?pdf and copy and write whatever notes you want i a note pad On Mon, Jul 28, 2008 at 10:54 AM, Rajasekaramya <[EMAIL PROTECTED]>wrote: > > hi there, > > I want to write the plots in the pdfs and the details about the graph in a > seperate notepad. > > plot(as.numeric(lapply(resultgenes,length)),

Re: [R] Overlay of simple plots

2008-07-28 Thread stephen sefick
?lines On Mon, Jul 28, 2008 at 12:57 PM, Phillip R. Hunt <[EMAIL PROTECTED]>wrote: > I am trying to create a graph that puts three x-y plots (age vs rate for > three categories ) on a single plot. SAS does this with the overlay > function. I cannot find in the R documentation a method for doing

[R] combining zoo series with an overlapping index?

2008-07-29 Thread stephen sefick
day<-structure(c(7.7, 7.7, 7.7, 7.7, 7.7, 7.71, 7.7, 7.71, 7.71, 7.7, 7.7, 7.7, 7.7, 7.69, 7.68, 7.68, 7.67, 7.67, 7.67, 7.66, 7.65, 7.65, 7.65, 7.64, 7.64, 7.63, 7.63, 7.63, 7.62, 7.62, 7.62, 7.62, 7.63, 7.63, 7.63, 7.63, 7.63, 7.64, 7.64, 7.65, 7.65, 7.65, 7.66, 7.66, 7.67, 7.67, 7.67, 7.68, 7.68

Re: [R] combining zoo series with an overlapping index?

2008-07-29 Thread stephen sefick
ue, Jul 29, 2008 at 9:01 AM, stephen sefick <[EMAIL PROTECTED]> wrote: > > day<-structure(c(7.7, 7.7, 7.7, 7.7, 7.7, 7.71, 7.7, 7.71, 7.71, 7.7, > > 7.7, 7.7, 7.7, 7.69, 7.68, 7.68, 7.67, 7.67, 7.67, 7.66, 7.65, > > 7.65, 7.65, 7.64, 7.64, 7.63, 7.63, 7.63, 7.62, 7.62, 7.62, 7.

Re: [R] combining zoo series with an overlapping index?

2008-07-29 Thread stephen sefick
plot(replace(day, is.list(index(lin)), coredata(lin))) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide

Re: [R] combining zoo series with an overlapping index?

2008-07-29 Thread stephen sefick
My apologies it worked perfectly. On Tue, Jul 29, 2008 at 12:40 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > plot(replace(day, is.list(index(lin)), coredata(lin))) > -- Let's not spend our time and resources thinking about things that are so little or so large that all t

Re: [R] FFT - (STATS) - is this correct?

2008-07-30 Thread stephen sefick
#ts with a freq=1 d <- ts(x) #it seems that this is an oscillating signal of ~ a day plot(d) #spectrogram which shows that the main frequency component is ~ 0.7. this translates to about a day (1.3 or so) sec.pgram(d) On Wed, Jul 30, 2008 at 7:14 AM, <[EMAIL PROTECTED]> wrote: > Hello, > > I hav

Re: [R] john chambers R book

2008-07-30 Thread stephen sefick
amazon bonus bucks On Wed, Jul 30, 2008 at 4:05 PM, Marc Schwartz <[EMAIL PROTECTED]>wrote: > on 07/30/2008 01:27 PM [EMAIL PROTECTED] wrote: > >> Just to let everyone know: John Chambers R book is no longer at pre-order >> status on Amazon which is a good thing. The only problem >> is that now i

[R] S 3 generic method consistency warning please help

2008-07-30 Thread stephen sefick
I would like to include this in a package. The S3 methods on R CMD check says * checking S3 generic/method consistency ... WARNING window: function(x, ...) window.chron: function(data, day1, hour1, day2, hour2, ...) See section 'Generic functions and methods' of the 'Writing R Extensions' man

Re: [R] S 3 generic method consistency warning please help

2008-07-30 Thread stephen sefick
Mac OS X 10.5.4 R 2.7.1 the function works fine when I use it in a regular session. I can post source code and the like if it is necessary to diagnose the problem. If ya'll need anything else just tell me. thanks Stephen On Wed, Jul 30, 2008 at 10:40 PM, stephen sefick <[EMAIL P

Re: [R] S 3 generic method consistency warning please help

2008-07-30 Thread stephen sefick
changed data argument to x and it passes the check. FOR NOW. Thanks for all of the help On Wed, Jul 30, 2008 at 10:49 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > Mac OS X 10.5.4 > R 2.7.1 > the function works fine when I use it in a regular session. I can post > s

Re: [R] Is R's fast fourier transform function different from "fft2" in Matlab?

2008-07-31 Thread stephen sefick
z <- rnorm(5000) f <- fft(z) d <- fft(f, inverse=T) plot(z, d) z <- rnorm(5000) z.ts <- ts(z) f <- fft(z.ts) d <- fft(f, inverse=T) plot(z.ts, d) temp <- matrix(c(1,4,2, 20), nrow=2) d <- fft(temp) f <- fft(d, inverse=T) plot(temp, f) this, looks to me, to be the same. you have to take the

Re: [R] Is R's fast fourier transform function different from "fft2" in Matlab?

2008-07-31 Thread stephen sefick
z <- rnorm(5000) f <- fft(z) d <- fft(f, inverse=T) plot(z, d) z <- rnorm(5000) z.ts <- ts(z) f <- fft(z.ts) d <- fft(f, inverse=T) plot(z.ts, d) temp <- matrix(c(1,4,2, 20), nrow=2) d <- fft(temp) f <- fft(d, inverse=T) plot(temp, f) this, looks to me, to be the same. you have to take the

Re: [R] Is R's fast fourier transform function different from "fft2" in Matlab?

2008-08-01 Thread stephen sefick
y the signal length. plot((z.ts)-(d/5000)) Does this hold for longer signals, too? On Thu, Jul 31, 2008 at 11:30 PM, Rolf Turner <[EMAIL PROTECTED]> wrote: > > On 1/08/2008, at 2:56 PM, stephen sefick wrote: > >> z <- rnorm(5000) >> f <- fft(z) >> d <-

Re: [R] simple help request

2008-08-01 Thread stephen sefick
# I am sure there is a better way than this a <- c(1,6,10,11,3,5,10) group <- c(5,5,5,5,9,9,9) my.df <- cbind(a,group) a.5 <- subset(my.df, group==5) median(a.5[,1]) a.9 <- subset(my.df, group==9) median(a.9[,1]) # MASS 4 is a good book On Fri, Aug 1, 2008 at 1:30 PM, Lotta R <[EMAIL PROTECTED]>

Re: [R] Time series, least squares line

2008-08-05 Thread stephen sefick
do you want just a date time stamp or do you want other things. There are many ways to skin this cat, but we are going to need more information... as a first start try looking at the zoo package. stephen On Tue, Aug 5, 2008 at 5:47 PM, Gareth Campbell <[EMAIL PROTECTED]> wrote: > Hello, > > I

[R] R CMD check message that I can't figure out (.onLoad failed)

2008-08-05 Thread stephen sefick
I have been able to track and fix everything so far with R CMD check. How do I make this reproducible and I will do it. this is the message that I got: Error in fun(...) : couldn't connect to display ":0.0" Error : .onLoad failed in 'loadNamespace' for 'tcltk' -- Let's not spend our time and r

[R] Diagnosing Documentation Issues

2008-08-06 Thread stephen sefick
is part of the package to see which particular one is the bugger 4. Whatever else anyone suggests I can upload the whole package somewhere so that this is reproducible if the list would find that helpful. Thanks, Stephen Sefick Southeastern Natural Sciences Academy Research Scientist (Stream Ecology

Re: [R] Diagnosing Documentation Issues

2008-08-06 Thread stephen sefick
using R CMD Rd2dvi --pdf foo.Rd creates pdfs for everyone of the files separately. What do I do now? On Wed, Aug 6, 2008 at 1:49 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > Mac OS X 10.5.4 > intel processor > R 2.7.1 > MacTex - latest version (yesterday and checked f

Re: [R] problem running R

2008-08-07 Thread stephen sefick
supply code and then we will probbly be able to figure it out On Thu, Aug 7, 2008 at 6:12 AM, Mahesh Pal <[EMAIL PROTECTED]> wrote: > Dear all, > I am new to R and very basic idea about it. > While running some codes i am getting following error: > > "Error in m[indexMax] : invalid subscript type

Re: [R] Diagnosing Documentation Issues

2008-08-07 Thread stephen sefick
worked just fine thanks- I have been let's say frustrated a little with this package building exercise, but I think I've got it for now ... thanks to everybody that helped Stephen On Wed, Aug 6, 2008 at 8:17 PM, Rolf Turner <[EMAIL PROTECTED]> wrote: > > On 7/08/2008, at 10

Re: [R] RPro

2008-08-08 Thread stephen sefick
To the best of my knowledge no R is not going private. Search the archives and you will gain insight into this matter. Stephen Sefick On Fri, Aug 8, 2008 at 10:13 AM, <[EMAIL PROTECTED]> wrote: > I recently came across a flyer from REvolution Computing, and I wanted to > ask if thi

Re: [R] How can I create a vector of vector?

2008-08-08 Thread stephen sefick
is this what you want? x <- c(31,2,3) y <- c(1,2,5,34,5656) z <- c(211,3243,5,343,3) a <- list(x,y,z) On Fri, Aug 8, 2008 at 1:59 PM, Jim <[EMAIL PROTECTED]> wrote: > Dear All, > > How can I create a vector of vector in R? > for example, I would input the data as follows: > x[1] <- c(31,2,3) > x[

Re: [R] box-plot without a box surrounding the plot

2008-08-09 Thread stephen sefick
par(bty="l") boxplot(count ~ spray, data = InsectSprays, col = "lightgray") #is that what you want? Stephen On Sat, Aug 9, 2008 at 8:12 PM, Jörg Groß <[EMAIL PROTECTED]> wrote: > Hi, > > I want to draw a boxplot and I want to get rid of the surrounding box. > So that there are only the axis lines

Re: [R] ANOVA

2008-08-10 Thread stephen sefick
I don't know of a way to do it without the data On Sun, Aug 10, 2008 at 11:35 AM, Angelo Scozzarella <[EMAIL PROTECTED]> wrote: > Hi, > > How can I make an ANOVA if I haven't got all data set but I know the > numbers of subjects for each group, the mean and di standard deviation for > each group

  1   2   3   4   5   6   7   8   9   >