Re: [R] Mean(s) from values in different row?

Hi > > Hello: > > Is there a way to get a mean from values stored in different rows? > > The data looks like this: > YEAR-1, JAN, FEB, ..., DEC > YEAR-2, JAN, FEB, ..., DEC > YEAR-3, JAN, FEB, ..., DEC > > What I want is the mean(s) for just the consecutive winter months: > YEAR-1.DEC

Re: [R] reporting multiple objects out of a function

Thanks, Gabor! When a beginner (like myself) asks a question, it seems that the thing that we believe we are confused about, or want to learn, may not be the thing that would actually help us the most if it were clearly understood. Your response is what I consider ideal: Answer my question, the

[R] Error in Penalized Command

How do I fix this error in the penalized command? L<-penalized(y,x,lambda1=10) Error in .checkinput(match.call(), parent.frame()) : the length of "response" (20) does not match the row count of "penalized" (1000) Thanks! [[alternative HTML version deleted]] __

Re: [R] About stepwise regression problem

chris, I'm not using lmer, i just use gam mixed with smoothing function and linear function and summary of the model, it shows Family: gaussian Link function: log Formula: newNO2 ~ pressure + s(maxtemp, bs = "cr") + s(avetemp, bs = "cr") + s(mintemp, bs = "cr") + RH + s(solar, bs = "cr") +

[R] ANOVA/ANCOVA Repeated Measure Mixed Model

Hello, I am trying to test some results I have for significance. It has been recommended that I use R and I am completely new to this. Set-up: Groups: two groups of 8 subjects (16 total) Two conditions: alert and passive Measurements: responses for three different stimuli (A, B, and C)

Re: [R] factors in probit regression

I need to quote David Winsemius on this one again: "The advancement of science would be safer if you knew what you were doing." Note that the whole model screams at you that it is wrongly modeled. You are running a fully interacted model with factor variables. Thus, you have 19 regressors plus the

Re: [R] ANOVA/ANCOVA Repeated Measure Mixed Model

You are modeling Condition * Stimulus * Group as fully interacted fixed effects. I do not actually think that you would need a complex random effect structure for this. A simple random effect for the individual might suffice. You could try to model this with lmer (in the lme4 library) and inspect w

Re: [R] Best Package to Extract Numbers from this Image?

You can probably use the package jpg , It will be a pain. You will have to read it in as a jpg then isolate the numbers ( by row and column ) then figure out some way ( probably not that hard) to identify the numbers from the bit pattern. I tried something similar with the png package. it wa

Re: [R] sum of functions

On 6 October 2011 16:20, Dimitris.Kapetanakis wrote: > Dear all, > > I would like to create a code for semiparametric Klein and Spady's > estimator. For that I created a function that provides the log-likelihood > function for each observation (so it is a function of betas and i, where i > denotes

Re: [R] Image Enhancing in R

I am not really sure what you mean, but perhaps you want to convert the image into RGB mode, and zero out the other two (GB) vectors? Look up ?rgb and allied functions: may help. Ranjan On Thu, 6 Oct 2011 18:58:19 -0500 Lorenzo Isella wrote: > Dear All, > I am not sure that R is exactly the r

Re: [R] R shell line width

Thank you very much for your responses! This is exactly what I needed. On Fri, Sep 16, 2011 at 8:13 PM, Joshua Wiley wrote: > Hi Mike, > > Look at ?options  particularly something like: > > options(width = 120) > > 80 is the default, I believe.  On 1920 pixels I can comfortably get > around 220 (

[R] barplots

Hello, I have somewhat of a weird data set and am attempting to create a barplot with it. I have 8 columns with different variables and their percentages. I have 1 column with representations of 4 different treatments the variables undergo. I also have 1 column with year the data was recorded. I

[R] Sets of general objects in R

Hi, Is there a simple package to manage sets of general objects, and perform simple operations like: - union - intersection - difference - test for membership... but where the equality test between 2 objects is under the control of the programmer? In the base package, it seems only operations on

[R] Coefficients for lagged plm model variables not calculated

Hello, So I am afraid I am having a recurring problem that I just can't figure out. I am using the plm package to conduct a panel analysis - although I am not sure if the problem is arising as a result of the plm package or something more general. I am trying to run a fixed effects model with eff

[R] linear classifiers with sparse matrices

I've been trying to get some linear classifiers (LiblineaR, kernlab, e1071) to work with a sparse matrix of feature data. In the case of LiblineaR and kernlab, it seems I have to coerce my data into a dense matrix in order to train a model. I've done a number of searches, read through the manuals

[R] Best Package to Extract Numbers from this Image?

I've done a bit of searching around on the forums looking for an R Package to extract the numbers from an image like the following: http://www.whnt.com/extras/whnt/weather/7day_640.jpg I would like to be able to do this in an automated way, so thanks for any hints and tips. [[alternativ

Re: [R] Wide to long form conversion

Hi: > I have some data 'myData' in wide form (attached at the end), and > would like to convert it to long form. I wish to have five variables > in the result: > > 1) Subj: factor > 2) Group: between-subjects factor (2 levels: s / w) > 3) Reference: within-subject factor (2 levels: Me / She) > 4)

Re: [R] Coefficients for lagged plm model variables not calculated

You probably should not do panel data analysis but multiple time series analysis as your T is much larger than N. You only have seven units of observation but some 150 observations on each unit. Also, if the values on each unit of observation are very close to each other between t and t-1, then yo

[R] Image Enhancing in R

Dear All, I am not sure that R is exactly the right tool for doing what I have in mind, but a look at EBImage and rimage convinced me that this should be doable with R. I have a color image and I would like to do the following: convert the image into a grayscale image retaining only the red col

Re: [R] Question about ggplot2 and stat_smooth

Dennis & Hadley, This does exactly what I need — thank you so much! Regards, Tom On 10/4/11 5:34 PM, Dennis Murphy wrote: Hi Hadley: When I tried your function on the example data, I got the following: dd<- data.frame(year = rep(2000:2008, each = 500), y = rnorm(4500)) g<- function(df, qs

Re: [R] Wide to long form conversion

Thanks for the pointer! I still couldn't figure out how to convert my data because the example at stackoverflow seems to have only one variable (Year) while there are three within-subject variables plus one between-subjects variable. Any further help? Gang On Thu, Oct 6, 2011 at 4:34 PM, Andrew M

Re: [R] non-cumulative hazard in Cox model with time-dependent covariates

On Fri, Oct 7, 2011 at 9:27 AM, David Winsemius wrote: > > On Oct 6, 2011, at 7:00 AM, koshihaku wrote: > >> Dear all, >> Is there a way to calculate the non-cumulative hazard (instantaneous >> hazard), which is the product of baseline hazard and exp{beta*covariate} ? >> I knew in survfit, we can

Re: [R] Duplicate elements of a vector

thanks guys! I also found this solution but I think yours is better. as.vector(rbind(c(1, 5, 23))) Cheers syrvn -- View this message in context: http://r.789695.n4.nabble.com/Duplicate-elements-of-a-vector-tp3879561p3879907.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Wide to long form conversion

Take a look here. http://stackoverflow.com/questions/2185252/reshaping-data-frame-from-wide-to-long-format Andrew Miles Department of Sociology Duke University On Oct 6, 2011, at 4:28 PM, Gang Chen wrote: > I have some data 'myData' in wide form (attached at the end), and > would like to conver

[R] Wide to long form conversion

I have some data 'myData' in wide form (attached at the end), and would like to convert it to long form. I wish to have five variables in the result: 1) Subj: factor 2) Group: between-subjects factor (2 levels: s / w) 3) Reference: within-subject factor (2 levels: Me / She) 4) F: within-subject fa

Re: [R] non-cumulative hazard in Cox model with time-dependent covariates

On Oct 6, 2011, at 7:00 AM, koshihaku wrote: Dear all, Is there a way to calculate the non-cumulative hazard (instantaneous hazard), which is the product of baseline hazard and exp{beta*covariate} ? I knew in survfit, we can get the estimator of cumulative baseline hazard, but how can we g

Re: [R] Urgent help needed for honours project - breaks between races in one year

This might help: if x is a vector of the race days then max(diff(sort(x))) finds the biggest gap between consecutive values. >>> From: "Jana.K" To: Date: 7/Oct/2011 4:47a Subject: [R] Urgent help needed for honours project - breaks between races in one year Hi to anyone who is willing to he

Re: [R] Weird behaviour of tab characters in a string in R (vs Python)

Things like hard-tabs are usually going to vary by text-editor / GUI. Python is pretty peculiar in its use of tabs, so I wouldn't expect R to replicate. My Matlab license is buggy right now, but I think you'd see similar behavior there, while interactive Ruby gives an autocomplete. I think that hit

Re: [R] Titles changing when a plot is redrawn

You can work around the problem by making each call to bquote() in a different environment, each containing its own value of 'i'. E.g., par(mfrow=c(2,1)) for(i in 1:2) { x <- 1:100 rmse <- sin(x/5) # fake data plot(x,rmse) str1 <- local({ i <- i ; bquote( paste("local RMSE(",thet

Re: [R] creating a loop for a function

To see what sort of things are available from Box.test() (and most everything else), hit ?Box.test and scroll down to "Value". There's actually more but its not called so frequently in my experience Michael On Thu, Oct 6, 2011 at 3:04 PM, upananda pani wrote: > Respected Sir, > I am grateful to

Re: [R] creating a loop for a function

Respected Sir, I am grateful to you for helping me out. Earlier i used to directly the formula without calling for statstic, parameter or p-value. As i am learning R first time, i can see how deep i have to go to learn it. With regards, Upananda On Fri, Oct 7, 2011 at 12:00 AM, R. Michael Weyla

Re: [R] Duplicate elements of a vector

Check function rep(), e.g., rep(c(1, 5, 23), each = 2) I hope it helps. Best, Dimitris On 10/6/2011 8:58 PM, syrvn wrote: Hi, let's assume I have the following vector a: 1 5 23 How can I use R to duplicate the elements so that my new vector looks like: 1 1 5 5 23 23 Many thanks, Syrv

Re: [R] Duplicate elements of a vector

rep(c(1,5,23), each=2) On Thu, Oct 6, 2011 at 2:58 PM, syrvn wrote: > Hi, > > let's assume I have the following vector a: > > 1 5 23 > > How can I use R to duplicate the elements so that my new vector looks like: > > 1 1 5 5 23 23 > > > Many thanks, > Syrvn > > -- > View this message in context:

Re: [R] counts in quantiles in and from a matrix

Excellent! Thank you! ben On Wed, Oct 5, 2011 at 9:18 PM, Dennis Murphy wrote: > Hi: > > Here's one way: > > m <- matrix(rpois(100, 8), nrow = 5) > f <- function(x) { >q <- quantile(x, c(0.1, 0.9), na.rm = TRUE) >c(sum(x < q[1]), sum(x > q[2])) >} > > t(apply(m, 1, f)) > > HTH, > De

[R] Duplicate elements of a vector

Hi, let's assume I have the following vector a: 1 5 23 How can I use R to duplicate the elements so that my new vector looks like: 1 1 5 5 23 23 Many thanks, Syrvn -- View this message in context: http://r.789695.n4.nabble.com/Duplicate-elements-of-a-vector-tp3879561p3879561.html Sent from

Re: [R] plot titles and labels XXXX

On Oct 6, 2011, at 2:51 PM, Dan Abner wrote: Hi everyone, Is it possible to to specify titles and axis labels on 2 separate lines when using the plot fn? What about tick mark labels? Example of main title: ?Quotes try: title(main="Figure 1: Side-by-side Boxplots\n of Y by X") in

[R] TCGA expression data: plotting ....

Hi, I am new to R. I am trying to figure out how to graph expression data from the TCGA database. If I understand correctly the expression data I have downloaded is from a microarray using the AgilentG4502A. I've had trouble reading into R in the level I, level II, and the gene expression analysi

[R] plot titles and labels XXXX

Hi everyone, Is it possible to to specify titles and axis labels on 2 separate lines when using the plot fn? What about tick mark labels? Example of main title: Figure 1: Side-by-side Boxplots of Y by X instead of Figure 1: Side-by-side Boxplots of Y by X [[alternativ

Re: [R] Titles changing when a plot is redrawn

On Oct 6, 2011, at 2:29 PM, David Winsemius wrote: On Oct 6, 2011, at 8:25 AM, John Nolan wrote: Thank you for telling me a fix. But I still don't know if this behavior is what is intended. I used bquote(...) because the plotmath(...) help page refers to bquote and gives an example lik

[R] Ignore my last message

It was an editor interaction -- the changes weren't being saved. __ 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.html and provide commented, minim

[R] Sweave and bibliographies

I added a bibliographic reference to one of my .Rnw documents, to wit: \begin{thebibliography}{9} \bibitem{Jose} Jos{\'{e}} C. Pinheiro and Douglas M. Bates, \emph{Mixed-Effects Models in S and S-PLUS}, Springer, 2000. \end{thebibliography} just before the \end{document}. When I run Sw

[R] factors in probit regression

Hi to all of you, I'm fitting an full factorial probit model from an experiment, and I've the independent variables as factors. The model is as follows: fit16<-glm(Sube ~ as.factor(CE)*as.factor(CEBO)*as.factor(Luz), family=binomial(link="probit"), data=experimento) but, when I took a look to t

Re: [R] gefp() boundaries?

Thank you very much for the answer. If I take Poisson model and follow "Generalized M-fluctuation tests for parameter instability", A. Zeileis and K. Hornik, Statistica Neerlandica (2007) Vol. 61, N. 4, p. 500-501 (section 4.3): data("Boston") n <- 506; my.X <- as.matrix(cbind(1, Boston["crim"], B

Re: [R] creating a loop for a function

Well you said you wanted statistics from the test, but you didn't say which statistics you wanted: any of the following would work: sapply(1:10, function(i) Box.test (lfut, lag = i, type="Ljung")$statistic) sapply(1:10, function(i) Box.test (lfut, lag = i, type="Ljung")$parameter) sapply(1:10, fun

Re: [R] About stepwise regression problem

I hadn't seen that page Dennis, that makes the case much more succinctly than my anti stepwise ramblings! Furthermore, "pigpigmeow" if you are using a random effects model i.e lmer - where are you getting your p-values from? And what do they mean in this context? I would strongly advise using i

Re: [R] Advice in model construction

Hi Chris, Linear regression model of categorical variables is equivalent to anova. If one of estimates of coefficients is significant in lm, it is interpreted as samples of Family did not come from the same population and the signficant factor level is compared to the reference level. If you want

Re: [R] help with regexp

Thanks to all who replied! With all these possible solutions it will be hard to find the best one :-). --- Gabor Grothendieck schrieb am Mi, 5.10.2011: > Von: Gabor Grothendieck > Betreff: Re: [R] help with regexp > An: "Jannis" > CC: r-h...@stat.math.ethz.ch > Datum: Mittwoch, 5. Oktober, 20

Re: [R] About stepwise regression problem

Hi: Please read this: http://www.stata.com/support/faqs/stat/stepwise.html Dennis On Thu, Oct 6, 2011 at 6:41 AM, pigpigmeow wrote: > using AIC/BIC, I'm not know too much about this. I just know using p-value to > perform stepwise regression > > if I used p-value to perform multimodel stepwise

Re: [R] SPlus to R

On Oct 6, 2011, at 11:17 AM, Scott Raynaud wrote: I'm re-posting this since it did not appear at the end of the thread. Sorry for the inconvenience. Not sure why it's giving the message: An embedded and charset-unspecified text was scrubbed... As far as I know my replies are set up as t

Re: [R] any way to convert back to DateTime class when "accidental" conversion to numeric?

Hi Dr. Ripley, All, Thanks for the succinct advice! Perfectly what I needed! Jeff, Absolutely I agree that this is a dangerous path, and I would never consider doing it for something that needs to be robust. But in 'R' type casting is a bit messed up, so I've come to accept that someti

Re: [R] Running *slow*

Patrick is right, most of the time is probably taken up for the reasons documented in the (masterful) R Inferno, namely the rbind() calls. There is another problem though and it gets at the very core of R, and for that matter, all interpreted languages that I'm familiar with. I'll give a fairly el

Re: [R] Fitting parabolic function to data

Hi Duncan, Thanks for your reply. You saved my Day :) Thanks, Henri On Thu, Oct 6, 2011 at 12:07 PM, Duncan Murdoch wrote: > On 11-10-06 4:59 AM, Henri Mone wrote: >> >> Dear R users and experts, >> >> I want to fit a shifted parabolic function with the following >> functional form to my data

Re: [R] Running *slow*

Probably most of the time you're waiting for this you are in Circle 2 of 'The R Inferno'. If the values are numbers, you might also be in Circle 1. On 06/10/2011 13:37, Thomas wrote: Anyone got any hints on how to make this code more efficient? An early version (which to be fair did more than t

Re: [R] setting up the RPostgreSQL_0.2-0.tar

I got this message in a private message before and asked you to read the posting guide before posting here! I also asked you to specify R version and operating system assuming install.packages() does not work. And I still do not know what the error was you got. Uwe Ligges On 06.10.2011 17:13

Re: [R] converting commas to points

Hi Anna, have a look at ?write.csv2, which deals with the Excel conventions for CSV in German locale. cheers Am 06.10.2011 17:39, schrieb Anna Lee: > Hello everyone! > > I work with a german excell version which uses commas instead of > points for seperating decimal places. R work with points s

Re: [R] Limitations of R

On 10/06/2011 08:33 AM, Uwe Ligges wrote: On 06.10.2011 16:49, Alaios wrote: Dear all, I have a few binary files like 9Gb or even of 50Gb. Hi Alaios -- Maybe you have a particular domain you are interested in (e.g., high-throughput sequence analysis) and there are packages (e.g., at http:

[R] "undefined slot classes" in package code

Dear all, I am facing a problem (warning message) when building a package that I am unable to fix: * checking whether package ‘myRcppTest’ can be installed ... WARNING Found the following significant warnings: Warning: undefined slot classes in definition of "S4test": Rcpptest(class "Rcpp_

Re: [R] converting commas to points

Sorry guys, I allready found the solution. Excell showed some of the numbers in the format: 1,90053-E05 and R interpreted it as 1.9... I changed that in Excel Cheers, Anna Am 6. Oktober 2011 17:48 schrieb Uwe Ligges : > > > On 06.10.2011 17:39, Anna Lee wrote: >> >> Hello everyone! >> >> I work w

Re: [R] converting commas to points

On 06.10.2011 17:39, Anna Lee wrote: Hello everyone! I work with a german excell version which uses commas instead of points for seperating decimal places. R work with points so in order to be able to save my excell tables without changing the commas to points, whenever I load a table I type i

Re: [R] unique possible bug

The error I am referring to is in unique.c in Base R, it cannot accomodate greater than 2^29 values, even though it appears the overflow protection should be 2^30. The only relevance of the arules package is I was using it while I discovered this issue. Thanks, Patrick 2011/10/6 Uwe Ligges : >

Re: [R] converting commas to points

>From ?read.csv: read.csv2( file, header = TRUE, sep = ";", quote="\"", dec=",", fill = TRUE, comment.char="", ...) I think this is specifically set up for German decimal commas. Rgds, Rainer On Thursday 06 October 2011 17:39:46 Anna Lee wrote: > Hello everyone! > > I work with a ger

[R] initial value in ComboBox tkwidget

Hello, I took me very long to find out how to set the initial value in the combobox widget in tk (see the example below). My question is now : Why does "textvariable" has to be a tclVar while "values" can be a normal vector? My next question is: How could I have known this much earlier? Is there a

[R] setting up the RPostgreSQL_0.2-0.tar

I am a new R user, and am looking to call Postgres through R. I am not so technical and need help in setting up. I have downloaded the tar file but cannot get passed this step. Is there documentation for dummies to get working with this? I don't even know where I should save this file? [

Re: [R] About stepwise regression problem

using AIC/BIC, I'm not know too much about this. I just know using p-value to perform stepwise regression if I used p-value to perform multimodel stepwise regression, is it correct in the first message box? -- View this message in context: http://r.789695.n4.nabble.com/About-stepwise-regression-

[R] sum of functions

Dear all, I would like to create a code for semiparametric Klein and Spady's estimator. For that I created a function that provides the log-likelihood function for each observation (so it is a function of betas and i, where i denotes the observation). Now, in order to maximize the log-likelihood f

[R] Urgent help needed for honours project - breaks between races in one year

Hi to anyone who is willing to help, I have a csv. file which has 1999 horses as the rows and the age(in years) of the horse at each race as columns. Ive read this file into R and called it 'horses'. Im trying to find the longest break between each race in the horse's first year of racing. I alrea

[R] converting commas to points

Hello everyone! I work with a german excell version which uses commas instead of points for seperating decimal places. R work with points so in order to be able to save my excell tables without changing the commas to points, whenever I load a table I type in: read.table(..., dec = ",") only R puts

Re: [R] Limitations of R

On 06.10.2011 16:49, Alaios wrote: Dear all, I have a few binary files like 9Gb or even of 50Gb. I would like to ask you what are the known limits of the R for the data processing part, I have a system with a lot of RAM (500Gb) Really accessible from one core? Amazing. but I am not sure ab

Re: [R] Problem with .C

Quoting Uwe Ligges : I don't agree that it's overkill -- you get to sidestep the whole `R CMD SHLIB ...` and `dyn.load` dance this way while you experiment with C(++) code 'live" using the inline package. You need two additional packages now where you have to rely on the fact those are avai

Re: [R] KS test and theoretical distribution

On Thu, 6 Oct 2011, francogrex wrote: x <- runif(100) y <- runif(100) ks.test(x,y) Two-sample Kolmogorov-Smirnov test data: x and y D = 0.11, p-value = 0.5806 alternative hypothesis: two-sided ok I expected that, but: ks.test(runif(100), "runif") One-sample Kolmogorov-Smirn

Re: [R] SPlus to R

I'm re-posting this since it did not appear at the end of the thread.  Sorry for the inconvenience.  Not sure why it's giving the message: An embedded and charset-unspecified text was scrubbed...  As far as I know my replies are set up as text and not html.  I did some testing and I believe the

Re: [R] KS test and theoretical distribution

On Thu, 6 Oct 2011, francogrex wrote: x <- runif(100) y <- runif(100) ks.test(x,y) Two-sample Kolmogorov-Smirnov test data: x and y D = 0.11, p-value = 0.5806 alternative hypothesis: two-sided ok I expected that, but: ks.test(runif(100), "runif") One-sample Kolmogorov-Smirn

[R] KS test and theoretical distribution

> x <- runif(100) > y <- runif(100) > ks.test(x,y) Two-sample Kolmogorov-Smirnov test data: x and y D = 0.11, p-value = 0.5806 alternative hypothesis: two-sided ok I expected that, but: > ks.test(runif(100), "runif") One-sample Kolmogorov-Smirnov test data: runif(100) D =

Re: [R] unique possible bug

I see. For now: Yes, you need to change and recompile. I will take a look what was actually changed and will run some test cases. Best, Uwe On 06.10.2011 16:50, Patrick McCann wrote: The error I am referring to is in unique.c in Base R, it cannot accomodate greater than 2^29 values, even thou

Re: [R] Dealing with proportions

On Oct 6, 2011, at 10:33 AM, David Winsemius wrote: On Oct 5, 2011, at 4:08 PM, Sam wrote: Dear list, I have very little experience in dealing with proportions, i am sure this is a very simple question Sometimes making that claim in a group of statisticians can provoke an extended dis

[R] Limitations of R

Dear all, I have a few binary files like 9Gb or even of 50Gb. I would like to ask you what are the known limits of the R for the data processing part, I have a system with a lot of RAM (500Gb) but I am not sure about the "internal" limitations of the R. How long for example a vector can be? C

Re: [R] Dealing with proportions

On Oct 5, 2011, at 4:08 PM, Sam wrote: Dear list, I have very little experience in dealing with proportions, i am sure this is a very simple question Sometimes making that claim in a group of statisticians can provoke an extended discussion to which there will be no eventual agreement.

Re: [R] distance coefficient for amatrix with ngative valus

Did you read any of the comments I made regarding working examples, meaningful question asking, or replying to the entire list? If you look at the code, you'll see pco is just a very elementary wrapper for cmdscale, the author of which is active on this list and could have seen your question and r

Re: [R] SPlus to R

I did some testing and I believe the program is operating properly.  It takes some time to finish, especially as sample sizes get larger, but I seem to be able to reproduce the results from the original paper.  Right now I'm most interested in method 3.  I set nc=40 and d=.2 as in the paper.  Th

Re: [R] Help with wireframe graphics problem (newbie)

On 06.10.2011 00:50, David Wiley wrote: All, I've read several tutorials re: generating wireframes, but am clearly missing something. I have data along the lines of: tbl [1:10,] Visits Activity Course.Grade 1 17218.31 2 7 1120.67 3 9

Re: [R] Wilcox Test / Mann Whitney U Test

And I figured it out, sorry to bother the list. The normal approximation I was using is not accurate in the presence of ties. Sam On Thu, Oct 6, 2011 at 10:56 AM, Sam Stewart wrote: > So I checked it with the wilcox_test in the coin library, and got the > same result.  That makes me more confid

Re: [R] unique possible bug

On 05.10.2011 22:15, Patrick McCann wrote: Hi, I am trying to read in a rather large list of transactions using the arules library. You mean the arules package? It seems in the coerce method into the dgCmatrix, it somewhere calls unique. Unique.c throws an error when n> 536870912; howev

[R] Fwd: break.axis all range of data

But in this way for moving also the point have to add 2 at all the x point after 8. axis(1,at=c(1:8,12:39),labels=1:36) plot(1:36,1:36,xaxt="n") I increased the gap, and if I don increase the x for the points at X 9, 10 they will appear in the gap. HC On Thu, Oct 6, 2011 at 8:20 AM, Heverkuhn

Re: [R] Wilcox Test / Mann Whitney U Test

So I checked it with the wilcox_test in the coin library, and got the same result. That makes me more confident that I made a mistake, but still doesn't help me find it d = data.frame(value=c(dropouts,remain),group=c(rep("dropout",length(dropouts)),rep("remain",length(remain wilcox_test(valu

Re: [R] Problem with .C

On 06.10.2011 15:41, Steve Lianoglou wrote: Hi, 2011/10/6 Uwe Ligges: On 06.10.2011 14:51, Jan van der Laan wrote: An obvious reason might be that your second argument should be a pointer to int. As others have mentioned, you might want to have a look at Rccp and/or inline. The documentati

Re: [R] how to install a locally built package

On 06.10.2011 15:45, Bond, Stephen wrote: Uwe, Are u saying 1) I can add the new func in one of the existing .R files in the R dir?? Yes, sure, or add an additional file. Or 2) add a new .R file in the same dir and ignore the lack of a matching .Rd file? Yes. You can also generate a sing

Re: [R] how to install a locally built package

Uwe, Are u saying 1) I can add the new func in one of the existing .R files in the R dir?? Or 2) add a new .R file in the same dir and ignore the lack of a matching .Rd file? Thank you. Stephen Bond | Senior Analyst | Treasury Analytics | 416-956-3092 -Original Message- From: Uwe Li

Re: [R] Problem with .C

Hi, 2011/10/6 Uwe Ligges : > On 06.10.2011 14:51, Jan van der Laan wrote: >> >> An obvious reason might be that your second argument should be a pointer >> to int. >> >> As others have mentioned, you might want to have a look at Rccp and/or >> inline. The documentation is good and I find it much e

[R] Wilcox Test / Mann Whitney U Test

Hello List, I'm trying to prepare some lecture notes on non parametric methods, and I can't manually reproduce the results of the wilcox.test function for ordinal data. The data I'm using are from David Howell's website, available here http://www.uvm.edu/~dhowell/StatPages/More_Stuff/OrdinalChis

Re: [R] [R-SIG-Finance] AsOf join in R

On Thu, Oct 6, 2011 at 2:15 AM, Ulrich Staudinger wrote: > A bit late,  but here is what I always do: > > m = merge(bid, ask, tick) > m<- interpNA(m, method="before") > > intrepNA can also interpolate NAs in different ways, for example linearly. > Its not clear precisely what sort of objects bid,

Re: [R] gamm: problems with corCAR1()

Karel V ugent.be> writes: > > Dear all, > > I’m analyzing this dataset containing biodiversity indices, measured over > time (Week), and at various contaminant concentrations (Treatment). We have > two replicates (Replicate) per treatment. > I’m looking for the effects of time (Week) and conta

Re: [R] Usng MCMCpack, error is "initial value in vmmin is not finite"

On 05.10.2011 20:38, yiy83102 wrote: __ 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.html and provide commented, minimal, self-contained, reprod

Re: [R] break.axis all range of data

Thank you Jim It seems exactly what I was looking for :) Claudio On Thu, Oct 6, 2011 at 5:11 AM, Jim Lemon wrote: > On 10/06/2011 12:31 AM, Heverkuhn Heverkuhn wrote: > >> ..all the point from 8 to 13. >> >> On Wed, Oct 5, 2011 at 8:28 AM, Heverkuhn Heverkuhn* >> *wrote: >> >> The problem with

Re: [R] Running *slow*

?unique x <- matrix(c(1:6, 6:1),ncol=2) x.temp <- x dim(x.temp) <- NULL unique(x.temp) Michael On Thu, Oct 6, 2011 at 8:37 AM, Thomas wrote: > Anyone got any hints on how to make this code more efficient? An early > version (which to be fair did more than this one is) ran for 330 hours and >

Re: [R] how to install a locally built package

On 06.10.2011 15:10, Bond, Stephen wrote: Is there a way to put all R code in a single file? I have too many small files now, and it is inconvenient to edit (I still have to put everything in one buffer) and when I add just one new func I have to go through the process of manually editing al

Re: [R] how to install a locally built package

Is there a way to put all R code in a single file? I have too many small files now, and it is inconvenient to edit (I still have to put everything in one buffer) and when I add just one new func I have to go through the process of manually editing all help files one by one. When I put all the

Re: [R] Problem with .C

On 06.10.2011 14:51, Jan van der Laan wrote: An obvious reason might be that your second argument should be a pointer to int. As others have mentioned, you might want to have a look at Rccp and/or inline. The documentation is good and I find it much easier to work with. For example, your exampl

Re: [R] error messages

Joyce Flynt yahoo.com> writes: > > Im doing a project using Rcommander. I have a dataset with 5 variables, > BMI being one of them. I would like to > find the variance, standard deviation, range, median, etc of BMI. > However, everytime I type a command > into the script window, (such as var(B

Re: [R] Fwd: Re: Party extract BinaryTree from cforest?

-- Forwarded message -- Date: Wed, 5 Oct 2011 21:09:41 + From: Chris To: r-h...@stat.math.ethz.ch Subject: Re: [R] Party extract BinaryTree from cforest? I found an internal workaround to this to support printing and plot type simple, tt<-party:::prettytree(cf at ensembl

Re: [R] How to do operations on zoo/xts objects with Monthly and Daily periodicities

Thanks for your suggestion Gabor, it has helped me in my developments during the last few months. I'm now trying to go further by being able to do the following: I have a collection of xts time series that are spaced in different ways (weekdays only, monthly, quarterly etc.). Is there an efficie

  1   2   >