similar to: OT: Statistics question

Can someone please explain what the $minimum result of the optimize function actually is? I'm trying to optimize the function: fitIT<-function(ampFac,ts_wave1,ts_template){ template<-stretchWaveTime(ts_template,ampFac) fit<-calcFit(ts_wave1,template) return(fit) } with >optimize(f=fitIT,interval=c(0.5,4),ts_wave1=test.data[,1],ts_template=test.data[,1]) $minimum [1]

Hi everyone, I'm hoping you can give me some pointers. I have a requirement to draw multiple (103) xy line plots onto one output device. Ideally the plots should be displayed in a hexagonal grid (example at www.maladmin.com/example.jpg). I can calculate the locations for each waveform but am wondering how to create multiple plotting areas. I have come accross references to a package grid

Can someone please give me a pointer here. I have two matrices matA A B C 1 5 2 4 2 2 4 3 3 1 2 4 matB A B C 1 TRUE FALSE TRUE 2 FALSE TRUE TRUE 3 FALSE FALSE FALSE how do I extract all the values from matA where the coresponding entry in matB == TRUE (or FALSE), perferably in vector form. Many thanks tom

I was able to install the curl library with no problems. However, when I tried to install ggplot (install.packages("ggplot2")), I got the same message as earlier, that R can't load internet.so. Thanks for your help! On Wed, Apr 20, 2016 at 1:24 PM, Tom Wright <tom at maladmin.com> wrote: > apt-get install curl libcurl4-openssl-dev > > On Wed, 2016-04-20 at 09:36

Dear all, I am trying to use the rstable(n, alpha, beta, gamma = 1, delta = 0, pm = c(0, 1, 2))) function to generate positive stable random numbers. For positive stable distribution, beta==1 and alpha is in (0,1), which defines random variables with support (0, infinity). So, I used rstable(100, 0.5, 1) for an example. I found that this gives me some negative numbers. For example, >

Dear R users, This is most likely very basic question but I am new to R and would really appreciate some tips on those two problems. 1) I need to plot variables from a data frame. Because of some few high numbers my graph is really strange looking. How could I plot a fraction of the samples (like 0.1 (10%), 0.2 up to for example 0.6) on x axis and values 'boundaries' (like any value

Can someone tell me how I create a vector of numbers where the step isn't 1? i.e. x<-(0.0,0.5,1.0,1.5....) Thanks tom

I'm just starting to learn to use R, and although I'm seeing lots of functions aimed at doing orthodox statistical analyses, I don't see the same for Bayesian analyses. What support does R have for Bayesian statistics?

I'm not an engineer so I hope I'm using the correct terminology here. I have a recorded waveform that I want to apply low and high pass filters too, are tehre already R functions existing to do this or am I going to have to program my own? thanks for any pointers tom

Can someone please show me what I need to get something like this to work for(a in c(1:5)){ data$a<-c(a:10) } so that i end up with a structure data$1<-[1,2,3,4,5,6,7,8,9,10] data$2<-[2,3,4,5,67,8,9,10] data$3<-[3,4,5,67,8,9,10] data$4<-[4,5,67,8,9,10] data$5<-[5,67,8,9,10] thanks loads Tom

Can somebody please take a look at this and tell me whats going wrong? It seems to be parsing wronly around the 'if' statement and gives me a directory listing. Thanks in advance Tom N.B. datan is an invented dataset xvals<-c(1,0.4,0.2) datan<-data.frame(s1=c(3,4,5),s2=c(5,5,5),s3=c(21,55,34),s4=c(5,3,2)) datan$sint<-NA datan$sgrad<-NA for(icount in 1:dim(datan)[1]) {

Sorry to keep posting but I want to do this right and I'm hoping for some pointers I now have two time series objects which I need to subtract. Unfortunatly the two series dont have the same sample rates. When I try to subtract them avgSub<-avg1-avg2 The time series object is clever enough to object. So I guess I need to write a function for subtraction of the time series objects which

> On Apr 21, 2016, at 1:23 PM, Paul Tremblay <paulhtremblay at gmail.com> wrote: > > I was able to install the curl library with no problems. However, when I > tried to install ggplot (install.packages("ggplot2")), I got the same > message as earlier, that R can't load internet.so. Is the libcurl directory in your search path? David. > > Thanks for your

Hello, I have been googling for 2 days and I cannot find the answer in previous posts. I have a set of d-dimensional data elements (d=11 .. 14), each data point can be censored at different values both Lower-limit and upper limit. N = 2000 sets of vectors of D=11 data points per vector. Each of the N*D points can have different upper and lower limits. I "simply" want to fit a

Dear All, What is the preferred way to get Bayesian analysis results (such as those from MCMCpacki, MCMCglmm, and DPpackage) into LaTeX table automatically? I have been using the "apsrtable" package and similar functions in "memisc" package, but neither seems to handle MCMC output directly. Many thanks. Shige

Dear all i use R only a few days and don't understand the difference between fivenum(x) und summary(x). > x [1] 20.77 22.56 22.71 22.99 26.39 27.08 27.32 27.33 27.57 27.81 28.69 29.36 [13] 30.25 31.89 32.88 33.23 33.28 33.40 33.52 33.83 33.95 34.82 > fivenum(x) [1] 20.770 27.080 29.025 33.280 34.820 > summary(x) Min. 1st Qu. Median Mean 3rd Qu. Max. 20.77 27.14

LD_LIBRARY_PATH=/usr/lib/x86_64-linux-gnu/:$LD_LIBRARY_PATH source ~/.bash_profile I still get the same error. On Thu, Apr 21, 2016 at 7:45 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > > On Apr 21, 2016, at 1:23 PM, Paul Tremblay <paulhtremblay at gmail.com> > wrote: > > > > I was able to install the curl library with no problems. However, when

Dear list members, I have a statistical question, that doesn't belong to this list, and I apologise for that in advance but I would appreciate your help very much. Is there some convention for selecting the a level for significance testing in scientific (e.g. chemical processes) studies? Most people use the 0.05 level but I could not find a reference to justify this. Why not 0.01 or 0.1?

I'm trying to create an app using TclTk and R Can someone please explain how I bind a click event to the tree widget (http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/TreeWidget.html) Ideally I'd like to bind to particular elements in the tree but tkbind doesnt seem to work. thanks tom

I'm stuck on what I feel should be a minor problem. I have a dataseries obtained from a MS Access database that consists of a series of numbers seperated by carridge returns (\r) Currently this data is in R as mode numeric??? I want to separate this into a vector of the componant numbers. Perhaps a little code will help describe what I've got here!! library(rodbc) s_sql<-'SELECT