Re: [R] ON MAC, how to copy a plot on to Word document?

2008-09-07 Thread Nicky Chorley
2008/9/7 asdfjkl; [EMAIL PROTECTED]: Yes, I don't know how to copy the plot on Mac and paste on to Word because you can't right click on the graph and say copy as metafile. I'm so surprised I can't find any information about this anywhere on the Internet... The obvious way would be to save

Re: [R] Corrupted PDF files

2008-09-07 Thread Prof Brian Ripley
How are you executing the script? Quite possibly FAQ Q7.22 applies (it will it you use source(), for example). On Sat, 6 Sep 2008, Nathan Teuscher wrote: I have the following code that when executed from the command line works properly and produces a proper PDF. When the script is executed,

[R] XML - get node by name

2008-09-07 Thread Antje
Hi there, I try to rewrite some Java-code with R. It deals with reading XML files. I started with the XML package. In Java, I had a very useful method which gave me a node by using: name of the node index of appearance start point: global (false) / local (true) So, I could do something like

Re: [R] Test for equality of complicatedly related average correlations

2008-09-07 Thread Ralph79
Thank you very much, Adam. I have to get a bit more familiar with the model you propose in order to understand if it applies to my problem as well. My question is not really does time show a different effect but which one of two measures is more reliable: My respondents have completed exactly

Re: [R] loop

2008-09-07 Thread Davide Crapis
It's exactely what I was looking for. Thanks a lot -- View this message in context: http://www.nabble.com/loop-tp19346683p19356409.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] Mode value

2008-09-07 Thread Jim Lemon
Carlos Morales wrote: Hello everyone, I would like to know if there is any function to calculate the mode value, or I have to build one to do it. Hi Carlos, If you mean the mode of a sample from a discrete distribution, try Mode in the prettyR package. Jim

[R] R_USER - in which file should I include it?

2008-09-07 Thread Eduardo M. A. M.Mendes
Hello I am a newbie. I had my R upgraded from 2.7.1 to 2.7.2 and in doing so I decided to install all 2.7 versions under c:\program files\R\2.7 from now on (2.7.1 is located under .\2.7.1) Although I don't like the idea (I am running Vista), I have edited etc\Renviron.site to contain:

[R] Problem with starting and using R

2008-09-07 Thread Thomas Lo
Dear all, I encountered a problem on starting and using the R v 2.7.2 installation on my PC running Windows Vista and would appreciate your help. When R was first started, the Rgui returned several error messages: Error in structure(.Internal(Sys.getenv(as.character(x), as.character(unset)$

Re: [R] Label 2 groups in PCA different colours

2008-09-07 Thread Pedro Mardones
here is a simple approach to, for instance, plot scores for PC1 and PC2 using diff colors: scores - prcomp(yourdata)$x plot(scores[1:100,1], scores[1:100,2], pch = 20, col = blue) points(scores[101:200,1], scores[101:200,2], pch = 20, col = red) PM On Sat, Sep 6, 2008 at 11:44 PM, pgseye [EMAIL

[R] cohen's kappa

2008-09-07 Thread Weiwei Shi
Dear all, I have a question on Cohen's kappa: Assume I have two datasets, one has 500 objects, 10 methods and the other, 1000 different objects, 20 different methods. Could I compare between the two datasets to conclude the 10 methods are more concordant than the 20 ones by looking at some

Re: [R] ON MAC, how to copy a plot on to Word document?

2008-09-07 Thread John Kane
I think you need to save the plot and import it into Word. AFAIK you can only copy and paste a plot in Windows. Have a look at ?png (There are other formats available) --- On Sun, 9/7/08, asdfjkl; [EMAIL PROTECTED] wrote: From: asdfjkl; [EMAIL PROTECTED] Subject: [R] ON MAC, how to copy

Re: [R] XML - get node by name

2008-09-07 Thread Ajay ohri
well not sure how its done in R , but heres a way to do it in simple Excel. http://decisionstats.com/2008/parsing-xml-files-easily/ Parsing XML files easily To parse a XML (or KML or PMML) file easily without using any complicated softwares, here is a piece of code that fits right in your excel

Re: [R] XML - get node by name

2008-09-07 Thread Dirk Eddelbuettel
On 7 September 2008 at 10:22, Antje wrote: | I try to rewrite some Java-code with R. It deals with reading XML files. I [...] | Now, I'd like to do something like this in R. Most important would be to | retrieve a node just by its name, not by the whole path. How is it possible? | | Can

[R] Fwd: request: most repeated sequnce

2008-09-07 Thread jim holtman
-- Forwarded message -- From: jim holtman [EMAIL PROTECTED] Date: Sun, Sep 7, 2008 at 11:42 AM Subject: Re: [R] request: most repeated sequnce To: Muhammad Azam [EMAIL PROTECTED] This should do it for you:

Re: [R] XML - get node by name

2008-09-07 Thread Gabor Grothendieck
In particular try this: Lines - ' + root + data loc=1 +val i=t1 22 /val +val i=t2 45 /val + /data + data loc=2 +val i=t1 44 /val +val i=t2 11 /val + /data + /root + ' library(XML) doc - xmlTreeParse(Lines, asText = TRUE, trim = TRUE, useInternalNodes = TRUE) root -

Re: [R] ON MAC, how to copy a plot on to Word document?

2008-09-07 Thread Dr Eberhard W Lisse
I just CMD-C'd it and pasted it into OpenOffice with CMD-V. el On 07 Sep 2008, at 16:57 , John Kane wrote: I think you need to save the plot and import it into Word. AFAIK you can only copy and paste a plot in Windows. Have a look at ?png (There are other formats available)

Re: [R] Axis tick label format and rotation

2008-09-07 Thread hadley wickham
Hi Kurt, Please tell me how to format data in a data frame so when currency amount is displayed in a chart the axis tick labels contain leading $ signs. The easiest way is add a custom scale: vals - seq(0, 100, by = 10) qplot(...) + scale_x_continuous(breaks = vals, labels = paste($, vals,

Re: [R] using nls to fit a curve to data

2008-09-07 Thread jpl
Thanks Ben! Switching over to the gamma pdf and using the algorithm=plinear did the trick. jpl -- View this message in context: http://www.nabble.com/using-nls-to-fit-a-curve-to-data-tp19332210p19360761.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] XML - get node by name

2008-09-07 Thread Duncan Temple Lang
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi Antje Well, the XML package gives you a variety of ways to parse an XML document and manipulate it in R. Perhaps the approach that best matches the Java-style you outline is to use XPath to access nodes. To do this, you use doc =

[R] an error to call 'gee' function in R

2008-09-07 Thread Qinglin Wu
Dear List: I found an error when I called the 'gee' function. I cannot solve and explain it. There are no errors when I used the 'geeglm' function. Both functions fit the gee model. The project supervisor recommends me to use the 'gee' function. But I cannot explain to him why this error

Re: [R] re ferring to a group of vectors without explicit enumeration

2008-09-07 Thread shalu
Thanks very much. This is what I implemented. I found a similar example elsewhere too. It works fine now. jholtman wrote: I would suggest that you use a list to store the values since it is easier to create and reference: output - list() for (i in 1:10) output[[i]] - seq(i) output

[R] Regression with nominal data

2008-09-07 Thread soeren . vogel
Hi, y is nominal (3 categories), x1 to 3 is scale. What I want is a regression, showing the probability to fall in one of the three categories of y according to the x. How can I perform such a regression in R? Thanks for your help Sören __

Re: [R] Regression with nominal data

2008-09-07 Thread Dimitris Rizopoulos
check: help(multinom, package = nnet) I hope it helps. Best, Dimitris [EMAIL PROTECTED] wrote: Hi, y is nominal (3 categories), x1 to 3 is scale. What I want is a regression, showing the probability to fall in one of the three categories of y according to the x. How can I perform such a

Re: [R] XML - get node by name

2008-09-07 Thread Gabor Grothendieck
On Sun, Sep 7, 2008 at 12:10 PM, Gabor Grothendieck [EMAIL PROTECTED] wrote: In particular try this: Lines - ' + root + data loc=1 +val i=t1 22 /val +val i=t2 45 /val + /data + data loc=2 +val i=t1 44 /val +val i=t2 11 /val + /data + /root + ' library(XML) doc

Re: [R] XML - get node by name

2008-09-07 Thread Antje
Thanks a lot to Gabor and Duncan! I didn't know that XPath is a standard. I'll give it a deeper look to better understand it. Oh, I guess I understand a bit more xpathApply(doc, //val, function(n) xmlValue(n)) would search globally for all nodes named val and return its values :-) So that's

[R] Averaging 'blocks' of data

2008-09-07 Thread Steve Murray
Dear all, I have a large dataset which I hope to reduce in size, to make it more useable. I hope to do this by taking an average of each 60 x 60 blockof values and forming a new data frame out of the averaged values. How would I go about taking averages of 60 x 60 'blocks' in R, and cycling

Re: [R] XML - get node by name

2008-09-07 Thread Gabor Grothendieck
On Sun, Sep 7, 2008 at 2:56 PM, Antje [EMAIL PROTECTED] wrote: Thanks a lot to Gabor and Duncan! I didn't know that XPath is a standard. I'll give it a deeper look to better understand it. Oh, I guess I understand a bit more xpathApply(doc, //val, function(n) xmlValue(n)) or just

Re: [R] Averaging 'blocks' of data

2008-09-07 Thread Gabor Grothendieck
This was answered last month: http://tolstoy.newcastle.edu.au/R/e4/help/08/08/19091.html On Sun, Sep 7, 2008 at 3:32 PM, Steve Murray [EMAIL PROTECTED] wrote: Dear all, I have a large dataset which I hope to reduce in size, to make it more useable. I hope to do this by taking an average

Re: [R] Averaging 'blocks' of data

2008-09-07 Thread Dylan Beaudette
On Sun, Sep 7, 2008 at 12:32 PM, Steve Murray [EMAIL PROTECTED] wrote: Dear all, I have a large dataset which I hope to reduce in size, to make it more useable. I hope to do this by taking an average of each 60 x 60 blockof values and forming a new data frame out of the averaged values.

Re: [R] Averaging 'blocks' of data

2008-09-07 Thread Steve Murray
Gabor - thanks for your suggestion... I had checked the previous post, but I found (as a new user of R) this approach to be too complicated and I had problems gaining the correct output values. If there is a simpler way of doing this, then please feel free to let me know. Dylan - thanks, your

Re: [R] Averaging 'blocks' of data

2008-09-07 Thread Steve Murray
Gabor - thanks for your suggestion... I had checked the previous post, but I found (as a new user of R) this approach to be too complicated and I had problems gaining the correct output values. If there is a simpler way of doing this, then please feel free to let me know. Dylan - thanks, your

[R] how to draw a vertical line from points to x-axis

2008-09-07 Thread Anny Huang
Hello, I want to know how to draw a line connecting each point to the x-axis perpendicularly (i.e. a vertical line). abline(v=...) seems not to work for my purpose, because it runs over the data point. Can anyone help? Thanks. Anny [[alternative HTML version deleted]]

[R] creating a vignette

2008-09-07 Thread Edna Bell
Dear R Gurus: How would I create a vignette, please? Why would a vignette be better than examples, please? Thanks, Edna Bell __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] restructuring datset problem

2008-09-07 Thread jim holtman
This should do it for you: CODE NAME 13 aaa 23 aab 33 aac 44 bba 54 bbb 64 bbc 74 bbd 85 cca 95 ccb x.s - split(x$NAME, x$CODE) maxLine - max(table(x$CODE)) # pad out the lines x.pad - lapply(x.s, function(line){ + # convert to character +

Re: [R] how to draw a vertical line from points to x-axis

2008-09-07 Thread Peter Alspach
Anny Here's one way: plot(0:10, 0:10, pch=16) lines(rep(0:10, each=3), t(matrix(c(0:10, rep(c(0,NA), each=11)), ncol=3))) HTH Peter Alspach -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Anny Huang Sent: Monday, 8 September 2008 8:49 a.m.

Re: [R] Averaging 'blocks' of data

2008-09-07 Thread jim holtman
Here is a way to do it by reading in 60 lines at a time and computing the means: # create some test data n - 360 x - matrix(runif(360*16800), nrow=16800) cat(x, file=/tempxx.txt) # now process the data 60 lines at a time, averaging each 60x60 block result - matrix(0, nrow=6, ncol=280) nextLine

[R] R_USER - in which file should I include it?

2008-09-07 Thread Eduardo M. A. M.Mendes
Hello I am a newbie. I had my R upgraded from 2.7.1 to 2.7.2 and in doing so I decided to install all 2.7 versions under c:\program files\R\2.7 from now on (2.7.1 is located under .\2.7.1) Although I don't like the idea (I am running Vista), I have edited etc\Renviron.site to contain:

Re: [R] how to draw a vertical line from points to x-axis

2008-09-07 Thread Barry Rowlingson
2008/9/7 Anny Huang [EMAIL PROTECTED]: Hello, I want to know how to draw a line connecting each point to the x-axis perpendicularly (i.e. a vertical line). abline(v=...) seems not to work for my purpose, because it runs over the data point. Can anyone help? Thanks. If your x-axis is at

Re: [R] restructuring datset problem

2008-09-07 Thread Gabor Grothendieck
Try this: # read in data ensuring NAME is character, not factor Lines - CODE NAME + 13 aaa + 23 aab + 33 aac + 44 bba + 54 bbb + 64 bbc + 74 bbd + 85 cca + 95 ccb + DF - read.table(textConnection(Lines), header = TRUE, as.is = TRUE) DF$seq =

[R] run optim() on a list

2008-09-07 Thread Weidong Gu
Hi, I am at the end of my wit to figure out how to run the optim function on a list. Basically, I have a data set of three columns as Site, Pool and Positivity ( the full data set is copied at the end). I want to run the maximal likelihood estimation separately on subsets split by Site

[R] Request for advice on character set conversions (those damn Excel files, again ...)

2008-09-07 Thread Emmanuel Charpentier
Dear list, I have to read a not-so-small bunch of not-so-small Excel files, which seem to have traversed Window 3.1, Windows95 and Windows NT versions of the thing (with maybe a Mac or two thrown in for good measure...). The problem is that 1) I need to read strings, and 2) those

Re: [R] an error to call 'gee' function in R

2008-09-07 Thread Thomas Lumley
From the gee() help page: Data are assumed to be sorted so that observations on a cluster are contiguous rows for all entities in the formula. -thomas On Sun, 7 Sep 2008, Qinglin Wu wrote: Dear List: I found an error when I called the 'gee' function. I cannot solve and explain

Re: [R] creating a vignette

2008-09-07 Thread Duncan Murdoch
On 07/09/2008 4:51 PM, Edna Bell wrote: Dear R Gurus: How would I create a vignette, please? Why would a vignette be better than examples, please? You can create a vignette any way you like, but the most common way is with Sweave: write a document that is mainly LaTeX, but with R code

[R] Help parametric boot

2008-09-07 Thread ctu
Hi R users Is there any example for nonlinear parametric boot? I google it but I can't find it. I am interested in the parameter estimators of a nonlinear model. But I really don't know how to code it in the ran.gen statement (data set from ?nls) fm1 - nls(weight ~

Re: [R] an error to call 'gee' function in R

2008-09-07 Thread Qinglin Wu
Dear Thomas Lumley: According to your suggestion I have sorted the 'newdata' by id. When I print the results the error was still there. So I tried to sort the 'newdata' by y, x although I don't think it make sense. The error was still there. Here it was the part of the code:

Re: [R] extracting max row from data matrix

2008-09-07 Thread Jorge Ivan Velez
Dear Srini, Here is one way: # Data set x=read.table(textConnection(fruit weight 1 apple1.3 2 apple1.5 3 apple1.6 4 orange1.4 5 orange1.6),header=TRUE) x[tapply(x$weight,x$fruit,which.max),] apple orange 1.61.6 or Try also

Re: [R] Averaging 'blocks' of data

2008-09-07 Thread Robert A LaBudde
I'm not sure I exactly understand your problem, but if you are looking for a recursive algorithm for calculating the average by addition of one record only at a time, consider: y[k] = y[k-1] + (x[k] - y[k-1])/k, where y(0) = 0, k = 1, 2, ... At each stage, y[k] = (x[1]+...+x[k])/k. At

Re: [R] cohen's kappa

2008-09-07 Thread Weiwei Shi
one more question, The third value (kappa (2*PA-1)) is adjusted for prevalence using the method proposed by Byrt, Bishop and Carlin (1993) --- from ?cohen.kappa What does the prevalence refer to? On Sun, Sep 7, 2008 at 10:43 PM, Weiwei Shi [EMAIL PROTECTED] wrote: Dear all, I have a

Re: [R] Regression with nominal data

2008-09-07 Thread paulandpen
Soren, It sounds like you are new to R so I will refer you to some packages that I think some people would find more user friendly as beginners. Zelig is excellent. You could run a series of logistic regressions coding your dependent variables as follows (a versus b, a versus c, b versus c)

[R] Poisson Distribution - Chi Square Test for Goodness of Fit

2008-09-07 Thread saggak
Dear R-help,   Chi Square Test for Goodness of Fit     Problem Faced :   I have got a discrete data as given below (R script)   No_of_Frauds -c