Re: [R] Problem with toString

dvkirankumar wrote: > > I got a problem with* "toString()"* function > its behaving odd in some times most of the times its giving correct > results > but in peculier cases its giving results like *"æ]"* > ... > especially while i am working with c# +R-project > This has to do w

[R] Problem with toString

Hi all, I got a problem with* "toString()"* function its behaving odd in some times most of the times its giving correct results but in peculier cases its giving results like *"æ]"* especially while i am working with c# +R-project thats also after getting some error like after trying

[R] drop1: warning message

Dear R collegues, I am new to R and do not understad the warning messages, which appear when I am trying to apply type III procedure to my glm: inkm$Gamie glm.inkmZ drop1(glm.inkmZ,test="Ch") Single term deletions Model: InkMb ~ IncStart + Vol + Gamie + PC1 + PC2 + PC3 + SpotsN + PC1:SpotsN

[R] Revolutions blog: December roundup

I write about R every weekday at the Revolutions blog: http://blog.revolution-computing.com , and every month I post a summary of articles from the previous month of particular interest to readers of r-help. You can find older summaries at http://blog.revolution-computing.com/roundups . (By the way

Re: [R] Solving graph theory problems with R ? (minimum vertex cover)

Tal Galili schrieb: I just realized (after many discussion with friends), that I might need to solve a (classical) graph theory problem with R. My specific problem is called: Minimum vertex cover for a hypergraph

Re: [R] How to get minimum value by group

Try this: > aggregate(myframe[4], myframe[1:3], min) factor1 factor2 factor3 y 1 1 a x 1 2 2 a x 9 3 1 b x 5 4 2 b x 13 5 1 a y 3 6 2 a y 11 7 1 b y 7 8 2

Re: [R] Problem about Box-Cox transformation (topic in html form)

Thank you, now I understand. If I plot the distribution of c888.dl.ma080+1200, then i will get a normally looked histogram. - -- Saji Ren from Shanghai China GoldenHeart Investment Group

Re: [R] how to track a number in a row

On Jan 11, 2010, at 10:16 PM, David Winsemius wrote: On Jan 11, 2010, at 9:38 PM, Márcio Resende wrote: Hi, I have a 100x15 matrix and in each row a set of 15 random numbers out of 25. for example: b <- c(1:25) a <-matrix(0,10,15) for (i in 1:10){ a[i,] <- sample(b,15,replace = FALSE) }

Re: [R] how to track a number in a row

On Jan 11, 2010, at 9:38 PM, Márcio Resende wrote: Hi, I have a 100x15 matrix and in each row a set of 15 random numbers out of 25. for example: b <- c(1:25) a <-matrix(0,10,15) for (i in 1:10){ a[i,] <- sample(b,15,replace = FALSE) } I would like to create another matrix (25x100), for ex

Re: [R] How to get minimum value by group

Tena koe Justin Try aggregate(): e.g., aggregate(myframe[,4], myframe[,1:3], min) HTH Peter Alspach > -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of JustinNabble > Sent: Tuesday, 12 January 2010 1:58 p.m. > To: r-help

Re: [R] How to get minimum value by group

On Jan 11, 2010, at 7:58 PM, JustinNabble wrote: I'd like to get a long data set of minimum values from groups in another data set. The following almost does what I want. (Note, I'm using the word factor differently from it's meaning in R; bad choice of words) myframe = data.frame(facto

Re: [R] how to track a number in a row

>Therefore I need to track, for example, if number 1 is present in the first row (d[1,1]) (which would give >me an probability of 1 out of 1). Just to make it clear, I need to track on matrix "a" if number 1 is present in the first row to fill the spot d[1,1] where the first collumn is related to

[R] how to track a number in a row

Hi, I have a 100x15 matrix and in each row a set of 15 random numbers out of 25. for example: b <- c(1:25) a <-matrix(0,10,15) for (i in 1:10){ a[i,] <- sample(b,15,replace = FALSE) } I would like to create another matrix (25x100), for example "d" with the probability of each number from the fir

Re: [R] Help me! using random Forest package, how to calculate Error Rates in the training set ?

Thank you, Andy I just read a paper, and they try to compare error rate among oob, test set, and training set and throung a figure showing random forest is not overfitting. when error rate in the training set come to zero, and oob and test set error rate do not increase. I am just a beginn

[R] How to get minimum value by group

I'd like to get a long data set of minimum values from groups in another data set. The following almost does what I want. (Note, I'm using the word factor differently from it's meaning in R; bad choice of words) myframe = data.frame(factor1 = rep(1:2,each=8), factor2 = rep(c("a","b"),each=4, time

[R] HoltWinters Forecasting

Hi R-users, I have a question relating to the HoltWinters() function. I am trying to forecast a series using the Holt Winters methodology but I am getting some unusual results. I had previously been using R for Windows version 2.7.2 and have just started using R 2.9.1. While using version 2.7.2

[R] Credit Migration Matrix

Hi user, I would like to know how can I compute credit rating migration matrix using R. I have 10 years data (monthly rates for each firm) and I would like to compute 12 (and more) months ahead migrations. Any hints? Best Rick [[alternative HTML version deleted]]

Re: [R] Problem with S4 generic function print

Hi Sven -- Sven Laur wrote: > First, if I have an S4 class boo should I define S3 method for printing > print.boo(..) or should I define S4 method print(..). > I mainly need to define print function only because some user might call > it and then I would like to dispatch the call to > right functi

Re: [R] 2 parameter exponential distribution

Ben, Say you have some variable (x) ~ exp f(x|lambda,theta) where f(x|lambda, theta) = lamba*(e^(minus lambda*(x-theta))) Lambda is the inverse scale parameter and theta is a location parameter (I erred in calling theta a shape parameter in the original question). John On Mon, Jan 11, 2010 at 4

Re: [R] 2 parameter exponential distribution

John Westbury wrote: > Ben, > > Say you have some variable (x) ~ exp f(x|lambda,theta) where f(x|lambda, > theta) = lamba*(e^(minus lambda*(x-theta))) > > Lambda is the inverse scale parameter and theta is a location parameter > (I erred in calling theta a shape parameter in the original questi

Re: [R] warning inside loop

What type of error do you get if it has a false converge? Does it raise an error (if so, ?try). Is there some value you can test for? You need to provide more details as to what happens. On Mon, Jan 11, 2010 at 6:06 PM, Rense wrote: > > Hi, > > I'm running some data simulations using (mixed e

Re: [R] warning inside loop

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Rense > Sent: Monday, January 11, 2010 3:07 PM > To: r-help@r-project.org > Subject: [R] warning inside loop > > > Hi, > > I'm running some data simulations using (mixed effect

[R] warning inside loop

Hi, I'm running some data simulations using (mixed effects)* regression models that show difficulty to converge. Therefore, I seek a way of capturing warnings (of false convergence) inside a loop. Inside that loop, I modify data and estimate a model. I do so many times with slightly different mo

Re: [R] Getting a date out of an indice in a time series

You cannot faithfully map year and week to ts since years do not have the same number of weeks and ts can only represent regular series. If you wish to use ts for this and you want it to be faithful then use ts(x) and 1 will represent the first week, 2 the second, etc. Then if o is a Date class va

Re: [R] 2 parameter exponential distribution

John Westbury gmail.com> writes: > > Hello, > > I am new to R and am trying to figure out how to specify an exponential > distribution with scale and shape parameters. I can specify an exponential > distribution with an inverse scale but would like to know how to specify the > distribution wit

Re: [R] Illustrating kernel distribution in wheat ears

On 01/12/2010 12:41 AM, Carl-Göran CG. Pettersson wrote: Thanks a lot for the quick response! The suggested code worked fine up to a certain point: the actual plotting from the datasets... With the sample dataset in "samp" my code looks like this at the end: for(cultivar in 1:6) +

Re: [R] help needed to find zero areas in a vector

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of boshao zhang > Sent: Monday, January 11, 2010 1:20 PM > To: r-help@r-project.org > Subject: [R] help needed to find zero areas in a vector > > Dear Helpers: >   > I spend more th

[R] Solving graph theory problems with R ? (minimum vertex cover)

I just realized (after many discussion with friends), that I might need to solve a (classical) graph theory problem with R. My specific problem is called: Minimum vertex cover for a hypergraph (Please s

Re: [R] Eigenvectors and values in R and SAS

On Jan 11, 2010, at 2:31 PM, jso717 wrote: Hi, I was wondering if function eigen() does something different from the function call eigen() in SAS. I'm in the process of translating a SAS code into a R code and the values of the eigenvectors and eigenvalues of a square matrix came out to be

Re: [R] help needed to find zero areas in a vector

Try this: with(rle(x), lapply(lengths[values == 0], rep, x = 0)) On Mon, Jan 11, 2010 at 7:19 PM, boshao zhang wrote: > Dear Helpers: > > I spend more than half a day to solve this problem in R: > > Let x be a vector of a string of 0s and 1s, such as > x<-c(0,0,1,1,0,0,0,0,1,1,0,0,0,0). It can b

Re: [R] Eigenvectors and values in R and SAS

On 12/01/2010, at 8:31 AM, jso717 wrote: Hi, I was wondering if function eigen() does something different from the function call eigen() in SAS. I'm in the process of translating a SAS code into a R code and the values of the eigenvectors and eigenvalues of a square matrix came out to be

[R] help needed to find zero areas in a vector

Dear Helpers:   I spend more than half a day to solve this problem in R:   Let x be a vector of a string of 0s and 1s, such as x<-c(0,0,1,1,0,0,0,0,1,1,0,0,0,0). It can be a very long vector. How to sub vectors of 0s? In the above example, I would like get the vectors (0,0), (0,0,0,0), (0,0,0,0).

Re: [R] Calling FING.EXE under RGui.EXE for windows.

I just wanted to add that it is not only caused by the FIND.EXE, but I also noticed the same problem with DOSKEY.EXE. It too causes the same fault, but the redirection ( To: Duncan Murdoch Cc: Uwe Ligges , John Schexnayder/San Jose/i...@ibmus, r-h...@stat.math.ethz.ch Date: 01/10/2010 10:42 AM

Re: [R] Finally, the first R spam!

I got that same email from "Steve Blum": Hi folks, Which is the best statistical package out there. I am out of work and looking to get the best buck for my money. Tx - ed --j Max Kuhn wrote: ... and this one I received this morning: Hi folks, Which is the best statistical package

[R] Eigenvectors and values in R and SAS

Hi, I was wondering if function eigen() does something different from the function call eigen() in SAS. I'm in the process of translating a SAS code into a R code and the values of the eigenvectors and eigenvalues of a square matrix came out to be different from the values in SAS. I would also

Re: [R] Calling FING.EXE under RGui.EXE for windows.

On 11/01/2010 2:56 PM, John Schexnayder wrote: I just wanted to add that it is not only caused by the FIND.EXE, but I also noticed the same problem with DOSKEY.EXE. It too causes the same fault, but the redirection (pointed out with FIND. I imagine there are number of other executables that w

Re: [R] Limiting number of tickmarks in lattice bwplot

On Jan 11, 2010, at 2:42 PM, GL wrote: Have a simple bwplot with 24 ordered factors across the x axis. I would like to only label every 4th tick mark so that the labels fit. I tried scales=list(x=list(tick.number=6)), but I still seem to get 24 tickmarks and 24 labels. Full code is below

[R] 2 parameter exponential distribution

Hello, I am new to R and am trying to figure out how to specify an exponential distribution with scale and shape parameters. I can specify an exponential distribution with an inverse scale but would like to know how to specify the distribution with a shape parameter as well. thanks much, John

[R] Limiting number of tickmarks in lattice bwplot

Have a simple bwplot with 24 ordered factors across the x axis. I would like to only label every 4th tick mark so that the labels fit. I tried scales=list(x=list(tick.number=6)), but I still seem to get 24 tickmarks and 24 labels. Full code is below: bwplot(SumOfIn.Use ~ Hour | Period, scales=li

[R] trouble with installing pbatR

Hello folks, I am having a problem installing pbatR package. (My R version is 2.8.1.) As I execute a command > install.packages(lib='/usr/local/apps/R-2.8.1/lib64/R/library', "pbatR") it detects lack of rootSolve package and attempts to install it. But cannot seem to compile it, and the res

[R] Getting a date out of an indice in a time series

I have a weekly data set imported via: tsSource=ts(sh1$I000,start=c(2004,1),freq=52) I am now getting to some 'spit and polish' but I realize something I can't wrap my head around. Given an outlier I find at say tsSource[54] ... how can get translate index 54 into the date\week. I mean I can figu

Re: [R] Illustrating kernel distribution in wheat ears

Here is how I tried to do a simple side-by-side boxplot: # Create new dataframe (cultivar, values, idx) # cbind goes by column so the cn column is repeated for us cultivars <- cbind(cultivar=wheat$cn, stack(wheat[2:15])) boxplot(values ~ cultivar, data=cultivars, main="Kernel distribution"

Re: [R] find the corresponding mean y

Annie, If I understand what you are asking, you want to sample some number of observations, with replacement, then find the mean of x and y for those observations. In this case, you would want to sample the row indices instead of the values. For example, #test data set.seed(1) df <- da

Re: [R] Directory operations

On 1/10/10 5:46 AM, anupam sinha wrote: Hi Uwe, Thanks for your suggestion . Here's my code. I am confused as to how to initialize an empty list . Here I have used pairlist() *list.files()->org_xml_dirs ## the following is much preferred: org_xml_dirs <- list.files() for

Re: [R] find the corresponding mean y

On Jan 11, 2010, at 12:42 PM, luciferyan wrote: Hello, I have 49 paired data, x, y. I have sampled x (where replacement is true), and find its mean. If you show how you did that, then we can show you the next steps. How can I find the corresponding mean y, which is the paired data of abov

[R] find the corresponding mean y

Hello, I have 49 paired data, x, y. I have sampled x (where replacement is true), and find its mean. How can I find the corresponding mean y, which is the paired data of above sample x? Thank you very much, Annie -- View this message in context: http://n4.nabble.com/find-the-corresponding-mean-y

Re: [R] Standard errors from a randomization test?

On Jan 8, 2010, at 6:45 AM, jjh wrote: Hello- Is it possible to estimate standard errors for a multiple regression model using a randomization test approach? I have seen a lot on using the procedure to get a test statistic, but nothing that talks about getting actual standard errors. Is

Re: [R] interpolation

On Jan 11, 2010, at 11:49 AM, René Mayer wrote: My problem is that x values increas with y x is mostly decreasing in the order you presented: plot(x, type="l") until some point then the pattern reverses. The whole line which line? is a kind of U-shape with a right-buttom to middel-top d

Re: [R] R for windows 64 bit

On 11.01.2010 17:56, alessia matano wrote: Many thanks for all your suggestions. And do you think then there is any particular issue to increase the amount of memory to be used? It allocates up to 8Gb on my machine (that does not have more RAM). But note: it is experimental. Uwe Ligges

Re: [R] R for windows 64 bit

Many thanks for all your suggestions. And do you think then there is any particular issue to increase the amount of memory to be used? 2010/1/11 Henrique Dallazuanna : > Try this version (beta of development version): > > http://www.stats.ox.ac.uk/pub/RWin/Win64/R-2.11.0dev-win64.exe > > On Mon,

Re: [R] Help me! using random Forest package, how to calculate Error Rates in the training set ?

From: bbslover > > now I am learining random forest and using random forest > package, I can get > the OOB error rates, and test set rate, now I want to get the > training set > error rate, how can I do? > > pgp.rf<-randomForest(x.tr,y.tr,x.ts,y.ts,ntree=1e3,keep.forest > =FALSE,do.trace=1e2)

Re: [R] interpolation

My problem is that x values increas with y until some point then the pattern reverses. The whole line is a kind of U-shape with a right-buttom to middel-top diagonal at the end of it (a look at the plot makes it clearer). The interpolation (approx, spline) makes a zick-zack aut of it. What I

Re: [R] R for windows 64 bit

Try this version (beta of development version): http://www.stats.ox.ac.uk/pub/RWin/Win64/R-2.11.0dev-win64.exe On Mon, Jan 11, 2010 at 2:29 PM, alessia matano wrote: > Dear all, > > do you know if there is any particular version of R to implement with > windows 64 bit, in such a way to increase

Re: [R] R for windows 64 bit

On 11/01/2010 11:29 AM, alessia matano wrote: Dear all, do you know if there is any particular version of R to implement with windows 64 bit, in such a way to increase the amount of memory it can use? There are two that I know of. Revolution Computing sells one, built with a commercial co

[R] R for windows 64 bit

Dear all, do you know if there is any particular version of R to implement with windows 64 bit, in such a way to increase the amount of memory it can use? How should I increase the memory, and more importantly to set a higher max vector size? It still stops me saying "Could not allocate vector of

Re: [R] How to control spaces between axis, tick and label in xyplot or xYplot?

Thanks a lot for your help! I have tried Peter's code: it works perfectly. Hopefully, this method will also be helpful to others with similar problems. Best regards, Jianghua -- View this message in context: http://n4.nabble.com/How-to-control-spaces-between-axis-tick-and-label-in-xyplot-or-xYpl

Re: [R] How to control number of significant digits (figures) in y-axis?

Dear Peter, Many thanks! Your code works perfectly. I hope this method will also be helpful to other users with similar problems. Best regards, Jianghua Peter Ehlers wrote: > > Use the 'scales=' argument together with formatC: > > x <- 1:10 > y <- sample(10) > xyplot(y ~ x, >scales = list

Re: [R] Run R in server web

I haven't used any of these yet, and don't know how they would fit in with Tomcat, but a few R web options are rapache (http://biostat.mc.vanderbilt.edu/rapache/), R.rsp (http://cran.case.edu/web/packages/R.rsp/index.html) and CGIwithR (http://www.omegahat.org/CGIwithR/). -- John On 1/3/10, Tal G

Re: [R] sparseM and kronecker product_R latest version

I can't reproduce this on my mac: > sessionInfo() R version 2.10.1 (2009-12-14) x86_64-apple-darwin9.8.0 locale: [1] C attached base packages: [1] graphics grDevices datasets utils stats methods base other attached packages: [1] fortunes_1.3-7 quantreg_4.44 SparseM_0.83

[R] Plot uncertainty in lmer predictions

Hi everyone. I have just started using R a couple of months ago, so I have a lot to learn yet. Anyway, I have fitted a linear mixed model for my data, with variable intercept and slope. I would like to plot the median estimate (wich I already know how to do) and uncertainty in the regression by pl

[R] Exponential regression

There are a couple of points to keep in mind when doing a log-transform of an exponential model, such as -- y = a*exp(b*x) 1. The implicit statistical model is multiplicative in the error. The implied statistical model of the log transform is -- log(y) = log(a) + b*x +

Re: [R] Strange behaviour of as.integer()

On 1/11/2010 10:37 AM, Martin Maechler wrote: Magnus Torfason Of course there is still the problem that: > 1+1 == sqrt(2)*sqrt(2) [1] FALSE and gmp will not solve this . I don't know if there is an R-package for arbitrary-precision reals floating around, but probably not. Yes, there's pac

Re: [R] sparseM and kronecker product_R latest version

I had an earlier problem where cutting and pasting from a help window caused obscure problems. Prof Ripley suggested trying example() and that did show that the example on the help page "worked". I just ran example(`%x%`), as well as example(matrix.csr) and encountered no errors. I also tr

Re: [R] Functions for QUAIDS and nonlinear SUR?

Werner, I know that S-Plus package Finmetrics has a NLSUR function. This is a commercial package, but maybe if you write the authors asking for code only, or some hints... Rick -- From: "Werner W." Sent: Friday, January 08, 2010 10:21 PM To: Sub

Re: [R] Plotting numeric values against non numeric items

Check out the interaction.plot function. On Fri, 8 Jan 2010, lse1986 wrote: l > l > Hi i want do a line graph. l > l > My y axis contains numeric values. My x axis contains non numeric l > statements. l > l > This is what i want the graph to look like. l > l > When i try to plot this graph

[R] Forming Portfolios for Fama / French Regression

Hi mates, I have a problem. I am new to R and want to conduct the Fama/French asset pricing test. As I am from Germany, I cannot use the already computed factors from French's website, but need to compute them myself. So I have to sort a number of stocks into different portfolios using one facto

Re: [R] Directory operations

To initialize an list: foo <- list() ## and then extend it as in your code foo <- c(foo, list(x=1)) foo $x [1] 1 foo <- c(foo, list(y=1)) foo $x [1] 1 $y [1] 1 At 7:16 PM +0530 1/10/10, anupam sinha wrote: Hi Uwe, Thanks for your suggestion . Here's my code. I am conf

Re: [R] Strange behaviour of as.integer()

> "MT" == Magnus Torfason > on Thu, 07 Jan 2010 10:01:28 -0500 writes: MT> There have been some really great responses to this question. I would MT> like to make a minor contribution by throwing the 'gmp' package into the MT> mix. MT> On 1/7/2010 8:12 AM, Duncan Mu

Re: [R] LD50 and SE in GLMM (lmer)

Sorry for the delay in response. I had a somewhat similar need recently with the difference that I used a logit link for a bioassay. The design had different dose-response "replicates" that I modeled as "blocks". It looks like you are concentrating on estimation of fixed effects and thus the popula

Re: [R] sparseM and kronecker product_R latest version

Many thanks for it. However it is strange that when I put the numbers rather than ncol(R) (a matrix with ncol=36698) it worked. Look below > dim(res2) [1] 170471 25822 > D<- as.matrix.csr(0,nrow(tmpb),25822) > D<- as.matrix.csr(0,nrow(tmpb),ncol(res2)) Error in if (length(x) == nrow * ncol) x <-

Re: [R] example() leaves console prompting to press Enter for new graphs

On Jan 11, 2010, at 9:28 AM, Michael Friendly wrote: Env: Win XP, R 2.9.2] It's a minor annoyance, but I find that after I run example(), the console state of prompting to press Enter for each new graph remains as it is under example. It's more annoying that it seems difficult to turn thi

Re: [R] sparseM and kronecker product_R latest version

> "am" == alessia matano > on Mon, 11 Jan 2010 15:38:39 +0100 writes: am> Many thanks for your suggestions, am> fortunately when I closed R and opeedn it again, it did not give again am> that error. However, now there is something more strange happening, am> related ag

Re: [R] example() leaves console prompting to press Enter for new graphs

On 11/01/2010 9:28 AM, Michael Friendly wrote: Env: Win XP, R 2.9.2] It's a minor annoyance, but I find that after I run example(), the console state of prompting to press Enter for each new graph remains as it is under example. It's more annoying that it seems difficult to turn this behavior

Re: [R] apply a function down each column

On Jan 11, 2010, at 9:18 AM, Steve Lianoglou wrote: Hi, On Mon, Jan 11, 2010 at 8:41 AM, Laetitia Schmid wrote: Hello World, I have a function that makes pairwise comparisons between two strings. I would like to apply this function to my data (which consists of columns with different s

Re: [R] sparseM and kronecker product_R latest version

Many thanks for your suggestions, fortunately when I closed R and opeedn it again, it did not give again that error. However, now there is something more strange happening, related again to sparse matrix, and I am afraid it could concern memory problems (I put 4000 of memory limit). I am defining

[R] example() leaves console prompting to press Enter for new graphs

Env: Win XP, R 2.9.2] It's a minor annoyance, but I find that after I run example(), the console state of prompting to press Enter for each new graph remains as it is under example. It's more annoying that it seems difficult to turn this behavior off! In a fresh session: > getOption("exampl

Re: [R] apply a function down each column

Hi, On Mon, Jan 11, 2010 at 8:41 AM, Laetitia Schmid wrote: > Hello World, > I have a function that makes pairwise comparisons between two strings. I > would like to apply this function to my data (which consists of columns with > different strings) in the way that it compares the first with th

Re: [R] Help with Order

David , Duncan Thanks for the swift response. You guys hit the nail on the head. That's exactly what the problem was. All the best Steve - Original Message - From: "David Winsemius" To: "Duncan Murdoch" Cc: "Steve Sidney" ; Sent: Monday, January 11, 2010 3:49 PM Subject: Re: [R] He

Re: [R] Standard errors from a randomization test?

Standard errors for the coefficients. Thanks! Tal Galili wrote: > > Hi jjh, > I wasn't able to understand: for what statistic do you want the standard > error for? (The coefficients, the R^2 ?) > > Tal > > > Contact > Details:--

[R] Problem with S4 generic function print

First, if I have an S4 class boo should I define S3 method for printing print.boo(..) or should I define S4 method print(..). I mainly need to define print function only because some user might call it and then I would like to dispatch the call to right function instead. Second, it seems that

[R] sparseM and kronecker product_R latest version

Dear all, I just installed the new version of R, 2.10.1, and I am currently using the package sparseM. (I also use a 64 bit windows version) I got a problem that I never had: when I try to multiply with a kronecker product (%x%) two sparse matrixes I get the following message: Error in dim(x) <-

Re: [R] Illustrating kernel distribution in wheat ears

Hi: It wasn't clear to me precisely what you wanted, but here are a couple of ideas in the hope that it will help. I used ggplot2 for the graphics, so it requires some manipulation of your dataset from 'wide' format to 'long'. I also add an indicator for side of the ear (odd is side one (L?), even

Re: [R] sparseM and kronecker product_R latest version

On Jan 11, 2010, at 7:55 AM, Peter Ehlers wrote: Do you have the same problem with the example on the help page? ?'%x%-methods' Works for me on Windows Vista (32-bit OS) and R version 2.10.1 Patched (2010-01-05 r50896). -Peter Ehlers alessia matano wrote: Dear all, I just installed the new

Re: [R] Problem about Box-Cox transformation (topic in html form)

Dear Saji Ren, Dieter Menne has already pointed out that you lost the negative values in the transformation. Another point is that since you selected the transformation based on the "started" data c888.dl.ma080 + 1200, then you should transform c888.dl.ma080 + 1200 and not c888.dl.ma080. But as Di

Re: [R] Help with Order

On Jan 11, 2010, at 7:49 AM, Duncan Murdoch wrote: On 11/01/2010 7:37 AM, Steve Sidney wrote: Dear List As a fairly new R programmer I seem to have run into a strange problem - probably my inexperience with R After reading and merging successive files into a single data frame, I find that

Re: [R] interpolation

On Jan 11, 2010, at 7:44 AM, René Mayer wrote: Dear R-users, I have a complex line by xy-values (ordered by z). And I would like to get interpolated y-values on the positions of x = 0:600. How do I get the correct points? x = c (790,790,790,790,790,786,783,778,778,766,763,761,761,761,715

Re: [R] Illustrating kernel distribution in wheat ears

Thanks a lot for the quick response! The suggested code worked fine up to a certain point: the actual plotting from the datasets... With the sample dataset in "samp" my code looks like this at the end: > for(cultivar in 1:6) + boxplot(unlist(samp[samp$cn==cultivar[cultivar],]), +add=TRUE,a

[R] apply a function down each column

Hello World, I have a function that makes pairwise comparisons between two strings. I would like to apply this function to my data (which consists of columns with different strings) in the way that it compares the first with the second entry, and then the third with the fourth, and then the fift

Re: [R] Finally, the first R spam!

... and this one I received this morning: > Hi folks, Which is the best statistical package out there. I am out of work > and looking to get the best buck for my money. > > Thanks, > Mike. > I guess this is still better than those AOL disks. On Thu, Jan 7, 2010 at 3:49 AM, Liviu Andronic wrote

Re: [R] sparseM and kronecker product_R latest version

Do you have the same problem with the example on the help page? ?'%x%-methods' Works for me on Windows Vista (32-bit OS) and R version 2.10.1 Patched (2010-01-05 r50896). -Peter Ehlers alessia matano wrote: Dear all, I just installed the new version of R, 2.10.1, and I am currently using th

Re: [R] Help with Order

On 11/01/2010 7:37 AM, Steve Sidney wrote: Dear List As a fairly new R programmer I seem to have run into a strange problem - probably my inexperience with R After reading and merging successive files into a single data frame, I find that order does not sort the data as expected. I have mu

Re: [R] K-means recluster data with given cluster centers

That kmeans returns an error if there is an empty cluster is a bit of a nuisance. It should not be too difficult to get rid off the kmeans function for what you call "reclustering". You could write your own function that assigns every point of the new data to the closest initial center. That s

[R] interpolation

Dear R-users, I have a complex line by xy-values (ordered by z). And I would like to get interpolated y-values on the positions of x = 0:600. How do I get the correct points? x=c(790,790,790,790,790,786,783,778,778,766,763,761,761,761,715,628,521,350,160,134,134,129,108,101,93,111,161,249,288,243

[R] Help with Order

Dear List As a fairly new R programmer I seem to have run into a strange problem - probably my inexperience with R After reading and merging successive files into a single data frame, I find that order does not sort the data as expected. I have multiple references in each file but each file

[R] Non-inferiority log-rank test

Hi Group, I'm looking for a method to calculate sample size for a non-inferiority survival analysis. Initially I've used cpower() of Hmisc, switching the roles of type-I and type-II errors, as cpower() assumes no treatment difference as the null hypothesis. Now, I would like to introduce event ra

[R] K-means recluster data with given cluster centers

K-means recluster data with given cluster centers Dear R user, I have several large data sets. Over time additional new data sets will be created. I want to cluster all the data in a similar/ identical way with the k-means algorithm. With the first data set I will find my cluster centers and s

[R] sparseM and kronecker product_R latest version

Dear all, I just installed the new version of R, 2.10.1, and I am currently using the package sparseM. (I also use a 64 bit windows version) I got a problem that I never had: when I try to multiply with a kronecker product (%x%) two sparse matrixes I get the following message: Error in dim(x) <-

Re: [R] Simple proportion barchart

Thanks Dieter. I had played with it before but was having trouble with it. I seem to be doing better now I have discovered prop.table! Will From: Dieter Menne [via R] [mailto:ml-node+1010645-1539655...@n4.nabble.com] Sent: 10 January 2010 12:00 pm To: Will Parry Subject: Re: [R] Simple prop

Re: [R] Standard errors from a randomization test?

Hi jjh, I wasn't able to understand: for what statistic do you want the standard error for? (The coefficients, the R^2 ?) Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Heb

Re: [R] Point plot comparisions

Are you asking how to add points to a plot with points in it ? (for that there is ?points) If not, please try to elaborate the question, since I didn't fully get it. Contact Details:--- Contact me: tal.gal...@gmail.com | 972-5

  1   2   >