[R] Subsetting without partial matches

2009-01-31 Thread Jonathan Dushoff
David: Thank you for your very valuable response. In fact, I was trying to _avoid_ partial matching, not accomplish it. Subset is a _much_ better way of doing what I was trying to do. Humorously, however, your code also reproduces the mistake that brought me here, AFAICT. I think my code

Re: [R] Extracting coordinates for cluster::clusplot()

2009-01-31 Thread Mark Difford
Dear Friends, How can I extract the coordinates used in the plot? They are not made available. You can get them as follows: x - rbind(cbind(rnorm(10, 0, 0.5), rnorm(10, 0, 0.5)), cbind(rnorm(15, 5, 0.5), rnorm(15, 5, 0.5))) pp - pam(x, 2) ## The coordinates are in obj$scores

Re: [R] reshape with two time variables

2009-01-31 Thread Peter Dalgaard
hadley wickham wrote: On Fri, Jan 30, 2009 at 5:57 PM, Neil Stewart neil.stew...@warwick.ac.uk wrote: I have a data frame in wide format that I'd like to convert to long format. For example, in wide format I have: id A1B1A1B2A2B1A2B2 1 1 400 475 420

[R] Splitting a data frame with break points where factor changes value

2009-01-31 Thread Titus von der Malsburg
I have a data frame called s3. This data frame has a column called saccade which has two levels 1 and -1. head(s3$saccade, 100) [1] NA NA NA NA -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 [26] -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 1 1 1 [51]

Re: [R] Regression

2009-01-31 Thread Paul Hiemstra
Edwin Wibisono wrote: Hello, I have problem I have a, and b in regression then I can't plot x,y with (and) a, b lines can you help me ? thx [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Multidimensional scalling

2009-01-31 Thread Tomek Wlodarski
Dear Stephen, Thanks a lot! now I see that cmdscale is not the best option for my problem So I am wondering if you can advice me other method of MDS or different approach to my problem: I have matrix which describes distances between object and I would like to visualise this matrix onto 2D in

[R] Using ggplot2 I need to move the location of legend to on the plot

2009-01-31 Thread Jason Rupert
  Thanks again for the hints about adding the vertical line to the hist plot and in ggplot.  That worked great.    Based on that advice I've been flipping through the ggplot2 doc and ggplot-static\index.html webpage more looking for the answer to the next question.    Unfortunately, I haven't

Re: [R] Q about how to use Anova.mlm

2009-01-31 Thread pgseye
Thanks a lot John - appreciate your help. Regards, Paul -- View this message in context: http://www.nabble.com/Q-about-how-to-use-Anova.mlm-tp21739443p21757533.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

Re: [R] bug in R2WinBUGS

2009-01-31 Thread Mike Meredith
Problem with the change is that none of our old scripts work! If the model.file is in R's current working directory we either have to use working.directory=getwd() or specify a full path with model.file=file.path(getwd(), mymodel.bug). It would be really nice if 'bugs' looked in R's c.w.d. for

Re: [R] Multidimensional scalling

2009-01-31 Thread Titus von der Malsburg
Hi Tomek, have a look at R News, Volume 3/3, December 2003. There you find an article about different algorithms that are available in R. Titus On Sat, Jan 31, 2009 at 01:36:29AM +0100, Tomek Wlodarski wrote: now I see that cmdscale is not the best option for my problem So I am wondering if

Re: [R] Splitting a data frame with break points where factor changes value

2009-01-31 Thread jim holtman
Here is one approach. You might want to change NAs to 0 if you want them included in the split set.seed(1) x - sample(c(1, -1), 30, TRUE) x # determine where changes occur change - cumsum(c(0, diff(x) != 0)) change split(x, change) On Sat, Jan 31, 2009 at 5:25 AM, Titus von der Malsburg

[R] Test Driven Development in R

2009-01-31 Thread Jose Quesada
Hi, I wonder what kind of interest there is on Test Driven Development (TDD) in R. Test Driven Development consists of writing the test before the function, and iteratively build the function until it passes the test. Python and Ruby (specially Ruby) have very strong test-oriented cultures.

[R] Creating Sub-Matrix

2009-01-31 Thread Rofizah Mohammad
Hello, How do I get sub-matrix? Example, I would like to get matrix of size 3x7 from the matrix of size 4x7. Meaning that I try to exclude one row of the original matrix. Thanks -rofizah- [[alternative HTML version deleted]] __

Re: [R] Creating Sub-Matrix

2009-01-31 Thread Philipp Pagel
How do I get sub-matrix? Example, I would like to get matrix of size 3x7 from the matrix of size 4x7. Meaning that I try to exclude one row of the original matrix. Just exclude the row by negative indexing - e.g. foo[-2,] You may want to consider reading the Introduction to R where all these

[R] sas.get under Linux

2009-01-31 Thread Adrian Dusa
Dear all, I am trying to import a SAS file into R (in fact I only need the value labels from the formats file), using Hmisc package, but I get this error: my.sas - sas.get(/home/adi/3, fis1_sgg) sh: sas: not found Error in sas.get(/home/adi/3, fis1_sgg) : SAS job failed with status 32512 I

Re: [R] Creating Sub-Matrix

2009-01-31 Thread David Winsemius
Several ways, two are shown: matrix(1:(4*7), 4,7) [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,]159 13 17 21 25 [2,]26 10 14 18 22 26 [3,]37 11 15 19 23 27 [4,]48 12 16 20 24 28 mtx - matrix(1:(4*7), 4,7) mtx[-1,]

Re: [R] Creating Sub-Matrix

2009-01-31 Thread Rofizah Mohammad
Many thanks David :-) I think I got that idea :-) On Sat, Jan 31, 2009 at 2:16 PM, David Winsemius dwinsem...@comcast.netwrote: Several ways, two are shown: matrix(1:(4*7), 4,7) [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,]159 13 17 21 25 [2,]26 10 14 18

Re: [R] sas.get under Linux

2009-01-31 Thread Frank E Harrell Jr
Adrian Dusa wrote: Dear all, I am trying to import a SAS file into R (in fact I only need the value labels from the formats file), using Hmisc package, but I get this error: my.sas - sas.get(/home/adi/3, fis1_sgg) sh: sas: not found Error in sas.get(/home/adi/3, fis1_sgg) : SAS job failed

[R] This site may harm your computer - Google warning about cran website

2009-01-31 Thread Christos Argyropoulos
This is extremely annoying (and it smells of rat). As of this time (1013 am ET), Google issues the following nasty warning: This site may harm your computer when I try to visit any of the CRAN package repository sites or CRAN task view. The funny thing is that Google's Safe Browsing

Re: [R] This site may harm your computer - Google warning about cran website

2009-01-31 Thread Gabor Grothendieck
I am having this same problem with a huge number of sites, not just CRAN. On Sat, Jan 31, 2009 at 10:19 AM, Christos Argyropoulos argch...@hotmail.com wrote: This is extremely annoying (and it smells of rat). As of this time (1013 am ET), Google issues the following nasty warning: This site

Re: [R] This site may harm your computer - Google warning about cran website

2009-01-31 Thread Erik Iverson
What OS/browser versions are you guys using? Christos Argyropoulos wrote: This is extremely annoying (and it smells of rat). As of this time (1013 am ET), Google issues the following nasty warning: This site may harm your computer when I try to visit any of the CRAN package repository sites or

Re: [R] This site may harm your computer - Google warning about cran website

2009-01-31 Thread Barry Rowlingson
2009/1/31 Gabor Grothendieck ggrothendi...@gmail.com: I am having this same problem with a huge number of sites, not just CRAN. Google seem to have fixed it now. Almost every site had this warning attached. Don't panic. The internet is safe (unless Jen drops it...[1]). Barry [1] The IT Crowd

Re: [R] This site may harm your computer - Google warning about cran website

2009-01-31 Thread Christos Argyropoulos
I am currently logged in my Vista OS. The problem seems to be Firefox specific (I'm using Firefox 3.05), since the warning does not show up when I browse with IE. At this time I cannot log in to my Linux partition, so cannot tell whether Firefox has the same problem under Linux. Date: Sat,

Re: [R] sas.get under Linux

2009-01-31 Thread Ajay ohri
Hi, have you looked at the third party SAS language compilers WPS ( 600 dollars per desktop version http://www.teamwpc.co.uk/home/ ) and Carolina ( http://dullesopen.com/) http://dullesopen.com/ http://dullesopen.com/ if you need just base SAS. I think SAS institute existing products have

Re: [R] Test Driven Development in R

2009-01-31 Thread Tobias Verbeke
Jose Quesada wrote: Hi, I wonder what kind of interest there is on Test Driven Development (TDD) in R. Test Driven Development consists of writing the test before the function, and iteratively build the function until it passes the test. Python and Ruby (specially Ruby) have very strong

Re: [R] reshape with two time variables

2009-01-31 Thread Neil Stewart
Thank you all so much for your help. I've gone with repeated_measures_data.csv: sub,A1B1,A1B2,A2B1,A2B2 s1,400,475,420,510 s2,390,500,470,472 s3,428,512,555,610 s4,703,787,801,822 s5,611,634,721,705 s6,543,522,612,788 s7,411,488,506,623 s8,654,644,711,795 library(reshape)

Re: [R] This site may harm your computer - Google warning about cran website

2009-01-31 Thread Duncan Murdoch
On 31/01/2009 10:19 AM, Christos Argyropoulos wrote: This is extremely annoying (and it smells of rat). As of this time (1013 am ET), Google issues the following nasty warning: This site may harm your computer when I try to visit any of the CRAN package repository sites or CRAN task view. The

Re: [R] Using ggplot2 I need to move the location of legend to on the plot

2009-01-31 Thread Juliet Hannah
Try these options: p + opts(legend.position=top) p + opts(legend.position=c(0.5,0.5)) On Fri, Jan 30, 2009 at 5:53 PM, Jason Rupert jasonkrup...@yahoo.com wrote: Thanks again for the hints about adding the vertical line to the hist plot and in ggplot. That worked great. Based on that

Re: [R] problem using identify() after plot()

2009-01-31 Thread Uwe Ligges
Christopher W. Ryan wrote: I can't seem to use the point-and-click identify() function properly. I'm running R 2.5.1 (I know, I need to get around to upgrading) under Win XP. The problem is, when I click on a point on the graph, I get an error, no point within 0.25 inches. But in some areas, I

Re: [R] Test Driven Development in R

2009-01-31 Thread Spencer Graves
Hi, Jose: How does the following I have not studied Test Driven Development, but I routinely write documentation *.Rd files before I write code, and the documentation files include test cases in the \examples section. With \dontshow, I hide tests that I don't think would benefit

[R] Fwd: The Internet is an Error by Google

2009-01-31 Thread Ajay ohri
Dear List, In continuation of the Website Error Messages. Ajay www.decisionstats.com Looks like the boys from Mountain View did some testing for anti spam or denial of service ( depends on if you like/trust/dislike G) , and went live instead of sand boxing the tests... Article from

[R] interaction plot in R for factorial experiment

2009-01-31 Thread Swanton0822
hi, i did a 2^4 factorial experiment, and i got the following result: Effect t p (computer) p -484.52494328125 -5.64590926071629 0.0001 d -450.67095078125 -5.25142684568607 0.0001 pd 438.80508046875

Re: [R] bootstrapping in regression

2009-01-31 Thread Charles C. Berry
On Fri, 30 Jan 2009, Stephan Kolassa wrote: Hi Thomas, Thomas Mang schrieb: I have a question here: I am not sure if I understand your 'fit the full model ... to the permuted data set'. Am I correct to suppose that once the residuals of the reduced-model fit have been permuted and added

Re: [R] interaction plot in R for factorial experiment

2009-01-31 Thread John Fox
Dear Swanton0822, You could refit the model with the statistically significant high-order terms and try plot(allEffects(mod)) from the effects package, where mod is the model object. A caution, however: If you used the default treatment contrasts for the factors, then a lower-order term such as

Re: [R] bootstrapping in regression

2009-01-31 Thread David Winsemius
On Jan 31, 2009, at 1:27 PM, Charles C. Berry wrote: On Fri, 30 Jan 2009, Stephan Kolassa wrote: Hi Thomas, Thomas Mang schrieb: I have a question here: I am not sure if I understand your 'fit the full model ... to the permuted data set'. Am I correct to suppose that once the

[R] Model Based Bootstrap of an MA(1)-Modell (with R-code)!

2009-01-31 Thread Andreas Klein
Hello. I want to calculate percentile intervals for the coefficient of an MA(1)-Model, but it doesn't work. Code for model based bootstrap based upon a MA(1)-Modell and building a bootseries recursivley (takes around 4 minutes to compute): y -

[R] can't get package boot to load

2009-01-31 Thread Gary Smith
Hi, I am new to R and I'm totally confused by this problem. I'm trying to load data and run a simple correlation using corr() in the boot package. Yesterday my code worked. Today it can't find the function corr(). I've tried searching the web and haven't found the root of my problem. Apologies

[R] FW: can't get package boot to load

2009-01-31 Thread Gary Smith
Changing library(RODBC,boot) to library(RODBC) library(boot) seems to have solved the problem. _ _ From: Gary Smith [mailto:gary.smit...@comcast.net] Sent: Saturday, January 31, 2009 12:55 PM To: 'r-help@R-project.org'

[R] Question on Sweave-Latex and examples in the Sweave Manual

2009-01-31 Thread eugen pircalabelu
Hi List, I have a problem with using Latex and Sweave for creating a document. So I downloaded the Sweave manual from http://www.statistik.lmu.de/~leisch/Sweave/Sweave-manual.pdf and i have tried to replicate the example on pages 4-5, but i encounter the following problem: IT DOES NOT WORK. I

Re: [R] sub question

2009-01-31 Thread David Hajage
Thank you, it's perfect. david 2009/1/30 Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no David Hajage wrote: Hello R users, I have a string, for example: x - \t\tabc\t def This string can contain any number of tabulations. I want to replace each tabulation of the begining

Re: [R] sub question

2009-01-31 Thread Wacek Kusnierczyk
David Hajage wrote: Thank you, it's perfect. to extend the context, if you were to solve the problem in perl, the regex below would work in perl 5.10, but not in earlier versions of perl; another approach is to replace the unwanted leading characters with equally many replacement characters

Re: [R] sub question

2009-01-31 Thread Gabor Grothendieck
On Sat, Jan 31, 2009 at 4:46 PM, Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no wrote: David Hajage wrote: Thank you, it's perfect. to extend the context, if you were to solve the problem in perl, the regex below would work in perl 5.10, but not in earlier versions of perl;

Re: [R] sub question

2009-01-31 Thread Wacek Kusnierczyk
Gabor Grothendieck wrote: On Sat, Jan 31, 2009 at 4:46 PM, Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no wrote: David Hajage wrote: Thank you, it's perfect. to extend the context, if you were to solve the problem in perl, the regex below would work in perl 5.10,

Re: [R] sub question

2009-01-31 Thread Gabor Grothendieck
On Sat, Jan 31, 2009 at 6:01 PM, Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no wrote: Gabor Grothendieck wrote: On Sat, Jan 31, 2009 at 4:46 PM, Wacek Kusnierczyk waclaw.marcin.kusnierc...@idi.ntnu.no wrote: David Hajage wrote: Thank you, it's perfect. to extend the context, if

[R] Tunnelling X for R graphics

2009-01-31 Thread Adam D. I. Kramer
Dear colleagues, I run R on a few different machines, and view graphs and the like by tunnelling X through SSH to my local machine. This is useful for me because my local machine can't easily handle some of the data sets I work with. However, when an ssh connection dies, the

Re: [R] FW: can't get package boot to load

2009-01-31 Thread Adam D. I. Kramer
If you check ?library, you will see that the package argument only takes one packages. So your code was like saying library(package=RODBC, help=boot) which does not make sense...but which would indeed load RODBC and not boot. That is why your problem occurred and why your solution worked. --Adam

[R] embed in webpage

2009-01-31 Thread Ista Zahn
Hi all, I want to use R to provide participants in psychology experiments with their scores on various measures. I looked around and found several programs designed to do this. I settled on Rpad because it seemed easy to set up and use. I have it working (see

Re: [R] Tunnelling X for R graphics

2009-01-31 Thread Dylan Beaudette
Try starting your R session after starting a 'screen' session. Like this: $ screen $ R # do stuff, when taking a break do CTRL-A D to disconnect # use as normal See the man page for screen, it is basically a terminal multiplexer that can gracefully accommodate connection failures. If you get

Re: [R] Tunnelling X for R graphics

2009-01-31 Thread Erik Iverson
Dylan, Can you confirm that you've done that before while tunneling R graphics over ssh? I thought I'd tried and screen and failed with a similar message when using graphics as Adam had. I could be wrong though. Dylan Beaudette wrote: Try starting your R session after starting a 'screen'

[R] question on statistical test

2009-01-31 Thread ANJAN PURKAYASTHA
hi, this is more of a statistical methodology question than an R question. however, since there are quite a few expert statisticians in this forum i think i can expect to get some useful feedback. here is my problem: i am studying a set of genes. for each gene i have a set of probes (there are 2-5

[R] Simulation of linear lists

2009-01-31 Thread Christian Hoffmann
Hi there has anybody experimented with the implementation in R of linear lists, as described e.g. in books on Pascal (Jensen and Wirth), or Wirth, Algorithms and data structures ? See \url{www.dbnet.ece.ntua.gr/~adamo/csbooksonline/AD.pdf } Pointers would be welcome. Christian -- --

Re: [R] This site may harm your computer - Google warning about cran website

2009-01-31 Thread Hesen Peng
Hey all, Just now I ran into this update on Google Blog explaining the whole incident: http://googleblog.blogspot.com/2009/01/this-site-may-harm-your-computer-on.html On Sat, Jan 31, 2009 at 11:03 AM, Duncan Murdoch murd...@stats.uwo.ca wrote: On 31/01/2009 10:19 AM, Christos Argyropoulos

[R] Extracting Coefficients and Such from mle2 Output

2009-01-31 Thread Tom La Bone
The mle2 function (bbmle library) gives an example something like the following in its help page. How do I access the coefficients, standard errors, etc in the summary of a? x - 0:10 y - c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8) LL - function(ymax=15, xhalf=6) + -sum(stats::dpois(y,

Re: [R] Extracting Coefficients and Such from mle2 Output

2009-01-31 Thread David Winsemius
On Jan 31, 2009, at 9:13 PM, Tom La Bone wrote: The mle2 function (bbmle library) gives an example something like the following in its help page. How do I access the coefficients, standard errors, etc in the summary of a? ?coef ?vcov eeep. Further comment on etc not possible at this time.