similar to: Packages

Dear R Users, Suppose we are interested for generating a new vector ( x ) from a current vector (y) of length 1000 so that x includes the sum of every 5 values in y respectively from the first to the end of length y. The same length of y for x is desired, so that other 4 positions (indices) in x are filled out with NA. For generating such a new vector, I have no idea. I tried in some

Dear R Users, For running wavelet functions using dwt( ), modwt( ), and mra( ), a wavelet filter algorithm is applied. For all these functions, default is "la8" and other possibility is "haar". In related documents, another possibilities like as symlet and coiflet ... are not cited. Besides "la8" and "haar", which wavelet filters can be used? Thank

Hi R Users, I don't know how much is difficult my problem and even it is possible to solve in R or not. Given a vector with 2000 observations. I want to creat a new vector from that vector so that new vector be the sum of every 5 observations sequently. That is , each new observation is sum of every 5 sequent observations, from the first observation of previous vector to the end. So

I use smooth.spline(x, y) in package modreg and I would like to get value of penalized log likelihood and preferable also its two parts. To make clear what I am asking for (and make sure that I am asking for the right thing) I clarify my problem trying to use the same notation as in help(smooth.spline): I want to find the natural cubic spline f(x) such that L(f) = \sum_{k=1}{n} w[k](y[k] -

Dear R Users, Suppose we want to creat a new vector ( x ) from a current vector (y) of length 1000. The current vector y includes negative, zero and positive values. We want our new vector x includes the negative values in y, otherwise NA with the same length as y. For this, we have x=y[y<0] . Now x includes a subset of y with shorter length than y. With x=match(y,x) we would

Dear R Friends , I want to name my variables( more than 2 variables in a single plot) within a plot to distinct them from each other, but I cann't. How it is possible? I don't mean x and y axis using xlab or ylab. At the below , it follows some lines, only as an example that you could try please, if it is possible. I really thanks for your attention. Amir library(graphics) y<-

Hi dear R users, Is it basically correct that a problem is ( linearly on nonlinearly ) modeled so that the slack variable is bounded by an upper bound ? If so, how it can be handled and coded practically ? for example: x1+ x2 =< b so ----> x1 + x2 + s=b s=b- x1 - x2 b- x1 - x2 =< upper value But algorithms can not calculate b- x1 - x2 , because

Dear R users, Sometimes it is needed to read only some columns from a table, in particulare for high frequency data. How it is possible to read just some certain columns using read.table ( ). The reason could be keeping space in R and in particular accelerating in reading data when the number of rows are huge and some of them are not needed. Thank you very much, Amir

Hi Every Body, I don't know why some regression functions have no related prediction function. For example lowess, loess, lokerns, lpridge, lpepa, and ksmooth. What could help? Is there any global or wrapper function so that can help? Regards, Amir Safari --------------------------------- [[alternative HTML version deleted]]

Dear All , Since dissimilarity is one of quality measures in clustering , I'm trying to access to the sum of dissimilarity as a whole measure. But after running my data using CLARA I obtain : 1128 dissimilarities, summarized : Min. 1st Qu. Median Mean 3rd Qu. Max. 0.033155 0.934630 2.257000 2.941600 4.876600 8.943700 But I can not find the sum of dissimilarity.How can i

Dear R Users, Could somebody please compare the packages of unit root test ( Uroot, Ucra, tseries and fseries ) regarding the type of test ( without constant and trend, with constant , and with constant and trend ) ? Regards, Amir Safari --------------------------------- Rekindle the Rivalries. Sign up for Fantasy Football [[alternative HTML version deleted]]

Hi every body , I simply want to calculate the 2th power of a vector without changing the sign of values. How it is possible in R ? Thanks a lot for any idea. Amir --------------------------------- [[alternative HTML version deleted]]

Dear R Users, Before running ksmooth( ), a suitable bandwidth selection is needed. I use some functions for this task and receive these results for my data: width.SJ(y,nb=100,method="ste") : 40.25 bcv(y,nb=100) : 40.53 ucv(y) : 41.26 bandwidth.nrd(y) : 45.43 After implementing the function ksmooth(x,y, bandwidth= each of abovementioned bandwidths), I have some NAs

Dear R Users, When using locpoly function, number of output values is smaller than the number of input values. How is it possible to get number of output component $y equal to the number of inputs. Thanks a lot, Amir --------------------------------- [[alternative HTML version deleted]]

Dear R Users, My question is overall and not necessarily related to R. Suppose we face to a situation in which MSE( Mean Squared Error) shows desired results but Bias shows undesired ones, Or in advers. How can we evaluate the results. And suppose, Both MSE and Bias are important for us. The ecact question is that, whether there is any combined measure of two above metrics. Thank you so

Hi Dear R users, For a pair plotting, usaully we use par( ) function. Apparently it does not work anywhere. I want to have 3 plots in a single figure, like this: par(mfrow=c(3,1)) densityplot( a) densityplot(b) densityplot(c) But it does not work. How is it possible to have such a figure with densityplot( ) in a single figure? So many thanks for any help. Amir Safari

I want to report the progress of my attempts to get R running on HP-UX with gcc: Firstly, I could work around the PIC problem by leaving out all the lib stuff (I rebuilt gcc with --enable-share, but this (alone) did not help), e.g. the line gcc -shared -fPIC -o ctest.sl ansari.o fexact.o kendall.o ks.o prho.o swilk.o -L/usr/local/lib -lg2c -lm -L/lib/pa1.1 -L/usr/lib/pa1.1 -u main

Hello, I had to set CPPFLAGS and FC=fort77 but configure saw atlas and used it instead of blas. However atlas for hp seems to missing some libs. I had to add -lblas in Makeconf. Everything builds ok. But modreg and tcltk doesn't load. R : Copyright 2000, The R Development Core Team Version 1.2.0 Under development (unstable) (2000-August-9) R is free software and comes with ABSOLUTELY NO

Full_Name: George N. White III Version: R-1.3.0 OS: Irix 6.5 Submission from: (NULL) (142.176.61.253) R-1.3.0 was built on Irix 6.5 using gcc and MIPSPro f90. This combination has worked for previous versions, in particular R-1.2.3. $ cd test/Examples $ ../../bin/R --vanilla < modreg-Ex.R > modreg-Ex.Rout This worked in R-1.2.3. It looks like there were a number of cosmetic changes

Hi, I'm not sure what the recommended thing to do here would be: In the lattice library, panel.loess needs to use loess.smooth(), which is in the modreg library. Now should I (1) make lattice depend on modreg, and call require(modreg) in zzz.R OR (2) use autoload("loess", "modreg") I like the second option better, as a particular session might not need panel.loess.