Re: [R] importing sas datasets

2007-04-27 Thread [EMAIL PROTECTED]
Hi Dan, thanks for the help, I'll try to resolve the sas.get problem in a second moment , actually I bypassed the problem creating a .csv file from my dataset and then using the read.csv function to import it into R. The only problem is that R imports only a 1000 observations, is that normal? Or

[R] Fwd: RE: importing sas datasets

2007-04-27 Thread [EMAIL PROTECTED]
Hi Dan, sorry for the precedent email, I created a csv file of only 1000 records by mistake, that's why R has only read 1000 records. regards Anna - Messaggio Inoltrato - Da : [EMAIL PROTECTED] [EMAIL PROTECTED] A : Daniel Nordlund [EMAIL PROTECTED],r-help@stat.math.ethz.ch Oggetto :

[R] Randomising matrices

2007-04-27 Thread Nick Cutler
I would like to be able to randomise presence-absence (i.e. binary) matrices whilst keeping both the row and column totals constant. Is there a function in R that would allow me to do this? I'm working with vegetation presence-absence matrices based on field observations. The matrices are

[R] Jarque-Bera and rnorm()

2007-04-27 Thread Murali Menon
Folks, I'm a bit puzzled by the fact that if I generate 100,000 standard normal variates using rnorm() and perform the Jarque-Bera on the resulting vector, I get p-values that vary drastically from run to run. Is this expected? Surely the p-val should be close to 1 for each test? Are

Re: [R] Jarque-Bera and rnorm()

2007-04-27 Thread Andrew Robinson
On Fri, April 27, 2007 7:02 pm, Murali Menon wrote: Folks, I'm a bit puzzled by the fact that if I generate 100,000 standard normal variates using rnorm() and perform the Jarque-Bera on the resulting vector, I get p-values that vary drastically from run to run. Is this expected? Surely the

[R] correlation tab

2007-04-27 Thread elyakhlifi mustapha
hello, I know how to use the function cor to do correlation frame but it remains a problem because I wanna calculate correlation between vector with differents dimensions it is possible? for exemple my data : V [[1]] [1] 97 104 103 98 99 92 99 99 98 87 95 101 91 104 98 96 84

Re: [R] Can a hyperlink be placed, thorugh R programming, on a graph?

2007-04-27 Thread Jim Lemon
Aldi Kraja wrote: Hi, If I use x-1:10 y-rnorm(10,0,1) ### pdf(file=c:\\aldi\\test.pdf) plot(x,y) segments(x,y,x+2,y+2) segments(x,y,x+0.5,y+0.5,col=3) ### dev.off() ### q() Is there a way that I can imbed in the graph plot for each point defined by x to x+0.5 and y to y+0.5

[R] Issues with Sweave and inclusion of graphics from nlme

2007-04-27 Thread caroline.truntzer
Dear all, I'm using the nlme package. I would like to introduce the plot of my groupedData object in my latex report using Sweave but it doesn't work. I have no error message in R (Sweave(my.file_Rnw)) and the plot is generated in R, but when I want to see the pdf (I'm using TexnicCenter) of my

[R] Inflate/Expand/Resize an array

2007-04-27 Thread Mark Payne
Gudday, I've had a good look everywhere trying to figure out how to do this, but I'm afraid I can seem to find an answer anywhere - maybe its because I'm not using the right terms, or maybe its because I'm a dummy. But unfortunately, I am not completely and utterly stuck. Here's the problem: I

Re: [R] correlation tab

2007-04-27 Thread Ted Harding
On 27-Apr-07 09:38:55, elyakhlifi mustapha wrote: hello, I know how to use the function cor to do correlation frame but it remains a problem because I wanna calculate correlation between vector with differents dimensions it is possible? for exemple my data : V [[1]] [1] 97 104 103

[R] bootstrap

2007-04-27 Thread Marc Bernard
Dear All, I would like to use a nonparametric bootstrap to calculate the confidence intervals for the 5% and 95% quantiles using boot.ci. As you know, boot.ci requires the use of boot to generate bootstrap replicates for my statistic. However this last function doesn't work in my case

Re: [R] Help with saptial analysis (cluster)

2007-04-27 Thread Roger Bivand
ONKELINX, Thierry Thierry.ONKELINX at inbo.be writes: Dear Fransico, The distance matrix would be 102000 x 102000. So it would contain 1040400 values. If you need one bit for each value, this would requier 9,7 GB. So the distance matrix won't fit in the RAM of your computer. Perhaps

Re: [R] Inflate/Expand/Resize an array

2007-04-27 Thread Robin Hankin
Hi f.dims - c(10,25,1,1,1,14) m.dims - c(10, 1,1,1,1,14) f - array(1:prod(f.dims),f.dims) m - array(1:prod(m.dims),m.dims) jj - do.call(abind,c(rep(list(m),25),along=2)) Then f + jj will give you what you want. The more general problem is a bit harder, as you say HTH rksh

Re: [R] Unwanted white borders on semi-transparent polygons?

2007-04-27 Thread Matthew Neilson
Thanks for your fast response. I'm using R version 2.1.1 on OS X 10.3.9 to create the pdfs. I have tried viewing the pdf output in both Acrobat 6 and 7 (both display a white border around each polygon) as well as Preview (displays fine). I have emailed the pdf file to some correspondents

Re: [R] Issues with Sweave and inclusion of graphics from nlme

2007-04-27 Thread Duncan Murdoch
On 4/27/2007 5:39 AM, [EMAIL PROTECTED] wrote: Dear all, I'm using the nlme package. I would like to introduce the plot of my groupedData object in my latex report using Sweave but it doesn't work. I have no error message in R (Sweave(my.file_Rnw)) and the plot is generated in R, but when I

[R] add arrows to barchart with groups

2007-04-27 Thread GOUACHE David
Hello Rhelpers, I am trying to represent the following data (hereafter named donnees) in a barchart with a grouping variable : sitetraitement dates res 1 NT 17/10/2005 normal 76.2 1 T 17/10/2005 normal 103.2 1 NT 23/11/2005

Re: [R] Inflate/Expand/Resize an array

2007-04-27 Thread Robin Hankin
[replying to myself] it might be better to use sweep(f,(1:6)[-2],m,+) instead. On 27 Apr 2007, at 11:56, Robin Hankin wrote: Hi f.dims - c(10,25,1,1,1,14) m.dims - c(10, 1,1,1,1,14) f - array(1:prod(f.dims),f.dims) m - array(1:prod(m.dims),m.dims) jj -

Re: [R] Unwanted white borders on semi-transparent polygons?

2007-04-27 Thread Matthew Neilson
Hi Roger, Yes, that's the FAQ to which I was referring in my original email. Unfortunately, turning off smoothing for line art (as suggested in the FAQ) did not work for me. Many thanks, -Matt On Fri Apr 27 12:32 , Roger Bivand [EMAIL PROTECTED] sent: Matt: Is this by any chance a

Re: [R] Unwanted white borders on semi-transparent polygons?

2007-04-27 Thread Prof Brian Ripley
The posting guide says For questions about unexpected behavior or a possible bug, you should, at a minimum, copy and paste the output from sessionInfo() into your message. If you are using an old version of R and think it does not work properly, upgrade to the latest version and try

Re: [R] Randomising matrices

2007-04-27 Thread Finny Kuruvilla
Hi Nick, This way isn't the most elegant but works well, especially if the matrices aren't too large: # This function works on 2x2 arrays, randomizing them, but # preserving row and column totals shuffle_matrix - function(x) { nrow = dim(x)[1] ncol = dim(x)[2] rmargins - apply(x,1,sum)

Re: [R] Jarque-Bera and rnorm()

2007-04-27 Thread Martin Becker
Murali Menon schrieb: Folks, I'm a bit puzzled by the fact that if I generate 100,000 standard normal variates using rnorm() and perform the Jarque-Bera on the resulting vector, I get p-values that vary drastically from run to run. Is this expected? Yes. Surely the p-val should be close

[R] how to be clever with princomp?

2007-04-27 Thread Simon Pickett
Hi all, I have been using princomp() recently, its very useful indeed, but I have a question about how to specify the rows of data you want it to choose. I have a set of variables relating to bird characteristics and I have been using princomp to produce PC scores from these. However since I

[R] proportional sampling

2007-04-27 Thread tibi . codilean
Dear All, I wonder if you could help me. I have a table with a series of sites that I would like to sample from. The table has 5 columns: ID X Coordinate Y Coordinate Value Factor The conditions are that each site can be selected more than once and the probability of it being selected (or

Re: [R] Randomising matrices

2007-04-27 Thread Stéphane Dray
Your solution is not suitable for quite large matrices. There are several papers on the subject : e.g. : @ARTICLE{SD751, author = {Gotelli, N.J. and Entsminger, G.L.}, title = {Swap and fill algorithms in null model analysis:rethinking the knight's tour}, journal = {Oecologia}, year

Re: [R] proportional sampling

2007-04-27 Thread jim holtman
sample(df$ID, numOfTime, replace=TRUE, prob=df$Factor) On 4/27/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Dear All, I wonder if you could help me. I have a table with a series of sites that I would like to sample from. The table has 5 columns: ID X Coordinate Y Coordinate Value

Re: [R] Unwanted white borders on semi-transparent polygons?

2007-04-27 Thread Matthew Neilson
Hi Brian, Terribly sorry if I accidentally broke a rule. sessionInfo() produces the following: sessionInfo() R version 2.2.1, 2005-12-20, powerpc-apple-darwin7.9.0 attached base packages: [1] methods stats graphics grDevices utils datasets [7] base I cannot upgrade to the

Re: [R] Randomising matrices

2007-04-27 Thread Finny Kuruvilla
Stephanie Dray's comments are absolutely correct. This method will only work with small matrics. The below modification is much faster than the first version that I wrote, but still slow on large matrices. shuffle_matrix - function(x) { nrow = dim(x)[1] ncol = dim(x)[2] cmargins -

[R] partitioning variation using the Vegan CCA routine?

2007-04-27 Thread Matthew McArthur
Hello I am using Jari Oksanen's CCA routine from the Vegan package on some estuary data, following a technique applied in (Anderson, M.J. Gribble, N.A., 1998, Partitioning the variation among spatial, temporal and environmental components in a multivariate data set, Australian Journal of Ecology

Re: [R] proportional sampling

2007-04-27 Thread Ted Harding
On 27-Apr-07 12:15:29, [EMAIL PROTECTED] wrote: Dear All, I wonder if you could help me. I have a table with a series of sites that I would like to sample from. The table has 5 columns: ID X Coordinate Y Coordinate Value Factor The conditions are that each site can be selected

Re: [R] Randomising matrices

2007-04-27 Thread Gabor Grothendieck
Try ?r2dtable On 4/27/07, Nick Cutler [EMAIL PROTECTED] wrote: I would like to be able to randomise presence-absence (i.e. binary) matrices whilst keeping both the row and column totals constant. Is there a function in R that would allow me to do this? I'm working with vegetation

Re: [R] Randomising matrices

2007-04-27 Thread Gavin Simpson
On Fri, 2007-04-27 at 09:45 +0100, Nick Cutler wrote: Hi Nick (Been meaning to reply to your private email to me but I've been away on vacation and at meeting for several weeks) I would like to be able to randomise presence-absence (i.e. binary) matrices whilst keeping both the row and

Re: [R] partitioning variation using the Vegan CCA routine?

2007-04-27 Thread Gavin Simpson
On Fri, 2007-04-27 at 16:03 +1000, Matthew McArthur wrote: Hello I am using Jari Oksanen's CCA routine from the Vegan package on some estuary data, following a technique applied in (Anderson, M.J. Gribble, N.A., 1998, Partitioning the variation among spatial, temporal and environmental

[R] error returned by make check in R-2.5.0

2007-04-27 Thread Eric Thompson
Today I tried to install the R-2.5.0 (currently running R-2.4.1) on Mandriva Linux. The ./configure and make commands seem to run fine, but make check gives the following messages: running regression tests make[3]: Entering directory `/home/ethomp04/R-2.5.0/tests' running code in 'reg-tests-1.R'

[R] Update of xlsReadWrite package

2007-04-27 Thread Hans-Peter
The (*windows only*) xlsReadWrite package has been updated and the new version is available on CRAN. xlsReadWrite natively reads and writes Excel files (v97 - 2003). o Changes from v1.1.1 to v1.3.3 (versions between never really published) - ROWNAMES-support also for matrices (not only for

Re: [R] how to be clever with princomp?

2007-04-27 Thread Gavin Simpson
On Fri, 2007-04-27 at 12:58 +0100, Simon Pickett wrote: Hi all, I have been using princomp() recently, its very useful indeed, but I have a question about how to specify the rows of data you want it to choose. I have a set of variables relating to bird characteristics and I have been

[R] (no subject)

2007-04-27 Thread elyakhlifi mustapha
hi, I don't understand my error when I submit this program for(k in LR) { + donRep[[k]] - subset(don2, Id_Essai == 1006961 Id_Cara == 193 Id_Rep == k , select = Val_O) + M[[193]] - matrix(rep(1, 3*min(C)), nrow = min(C), ncol = 3) +

Re: [R] (no subject)

2007-04-27 Thread Matthew Neilson
Hey, The syntax is: if(condition){ #COMMANDS_A } else{ #COMMANDS_B } Hope this helps, -Matt On Fri Apr 27 15:07 , elyakhlifi mustapha [EMAIL PROTECTED] sent: hi, I don't understand my error when I submit this program for(k in LR) { + donRep[[k]] +

[R] [R-pkgs] Hmisc Version 3.3-2 is now in CRAN repository

2007-04-27 Thread Charles Dupont
Hmisc Version 3.3-2 has been posted to the CRAN repository. Change Log 3.3-2 4/26/2007 Fixed bug with combine- function not handling NULL values. 3.3-1 3/29/2007 Created functions trunc.POSIXt, ceil.POSIXt, round.POSIXt to do those ops on POSIXt objects. Removed chron

[R] Write text in the device region of a graphic device

2007-04-27 Thread Felix Wave
Hello, I started a graphic device: par(oma = c(2,0,0,0), mfrow=c(2,2) ) in the cols and rows are three images. Now I want to write a text in the device region, it's the main title of the graphic device. But when I use mtext() I can only write in the figure region of my four plots. Has anybody an

Re: [R] Write text in the

2007-04-27 Thread Matthew Neilson
Hey Felix, So basically what you want is a figure containing a block of four plots, with a main title for the figure? If that's the case then something like this should work: # BEGIN CODE # par(oma=c(0,0,1,0), mfrow=c(2,2)) for(i in 1:4){ plot(NA,xlim=range(0,10),ylim=range(-5,5))

[R] SweaveInput and absolute paths

2007-04-27 Thread Kevin R. Coombes
The path is correct, and the file exists. Here is a transcript of the R session: - sessionInfo() R version 2.4.0 (2006-10-03) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United

[R] R-2.5.0:BioCv2install

2007-04-27 Thread Anna Lobley
I am receiving the following errors when trying to install gcrma and a number of other Bioconductor packages when the installWithVers flag is set to be TRUE. source(http://www.bioconductor.org/getBioC.R;); getBioC(gcrma, installWithVers=T); Running biocinstall version 2.0.7 with R version 2.5.0

Re: [R] select if + other questions

2007-04-27 Thread Natalie O'Toole
Hi, Does anyone know how to skip variables (or columns) in R. Say, for example i had PUMFID position1 and Y_Q10A position 33 and i do not want to include all the variables in between. Is there a way to do this in R when you are extracting variables from a large .txt file with many, many

[R] [ANN:] new package xlsReadWrite*Pro*: announce + update.

2007-04-27 Thread Hans-Peter
I am happy to announce that a new version of xlsReadWritePro is available. xlsReadWritePro natively reads and writes Excel files (v97 - 2003) on Windows. The pro version is more powerful than the free version but, beeing a small company, we ask you to support our development and maintenance

Re: [R] A coding question involving variable assignments in ifelse()

2007-04-27 Thread xpRt.wannabe
Is this an ad hominem comment or a comment of brevity? Unless my eyes are playing tricks on me, I can't seem to find any language in the Posting Guide on what is considered a reasonable vs. unreasonable request from an anonymous poster. Kindly point me to it if it exists. In any case, thanks

Re: [R] A coding question involving variable assignments in ifelse()

2007-04-27 Thread Gabor Grothendieck
I don't think there is any requirement to identify yourself in any way nor should their be. Many people on the list are in academia and in those cases they probably want their name in lights but others may wish to have a lower profile and its common to use an alias on the net for privacy. On

Re: [R] select if + other questions

2007-04-27 Thread John Kane
You need to have a look at Chapter 5 of the Intro to R. I would recommend downloading the pdf and printing it out. It is not an easy read but it should help. newdata - mydata[, c(PUMDID, Y_Q10A)] or newdata - mydata[, c(1,33)] should do the trick --- Natalie O'Toole [EMAIL PROTECTED]

[R] Scatter Plot in R - Help

2007-04-27 Thread A Ezhil
Dear All, I am using the following commands to do the scatter plot of two vectors, say X and Y. plot(X,Y, col=blue) abline(a=1,b=1, col=red) abline(a=-1,b=1, col=green) I would like to split the scatter plot into 3 part with 3 different colors: (i) points lies between 2 lines, (ii) points above

Re: [R] A coding question involving variable assignments in ifelse()

2007-04-27 Thread Duncan Murdoch
On 4/27/2007 11:41 AM, Gabor Grothendieck wrote: I don't think there is any requirement to identify yourself in any way nor should their be. Many people on the list are in academia and in those cases they probably want their name in lights but others may wish to have a lower profile and its

Re: [R] A coding question involving variable assignments in ifelse()

2007-04-27 Thread rolf
I think Duncan's suggestion that the poster work it out for him/her self is perfectly reasonable. People can be anonymous all they like. Those of us who are made suspicious by this anonymity (what have they got to hide?) can decline to provide assistance. cheers,

Re: [R] select if + other questions

2007-04-27 Thread Natalie O'Toole
Hi John, I figured out the skipping of variables using a first subset then making a second subset it worked! Thanks, Natalie __ You need to have a look at Chapter 5 of the Intro to R. I would recommend downloading the pdf and printing it out. It is not an easy read but

Re: [R] Issues with Sweave and inclusion of graphics from nlme

2007-04-27 Thread Friedrich Leisch
On Fri, 27 Apr 2007 07:15:22 -0400, Duncan Murdoch (DM) wrote: On 4/27/2007 5:39 AM, [EMAIL PROTECTED] wrote: Dear all, I'm using the nlme package. I would like to introduce the plot of my groupedData object in my latex report using Sweave but it doesn't work. I have no error

[R] Logistic Regression Question: Risk Threshold

2007-04-27 Thread diana
Hi, I am working on problem 2 of Chapter 8 in Data Analysis and Graphics Using R and don't know how to approach the second half of the question: In the data set (an artificial one of 3121 patients, that is similar to a subset of the data analyzed in Stiell et al., 2001) head.injury, obtain a

[R] Example of mcmcsamp() failing with lmer() output

2007-04-27 Thread Michael Kubovy
Hi, I would appreciate help with the following model 1= gunload - read.table(hh('datasets/gunload.dat'), header = T) gunload$method - factor(gunload$method, labels = c('new', 'old')) gunload$physique - factor(gunload$group, labels = c('slight', 'average', 'heavy')) gunload$team9 -

Re: [R] Scatter Plot in R - Help

2007-04-27 Thread J . delasHeras
Quoting A Ezhil [EMAIL PROTECTED]: Dear All, I am using the following commands to do the scatter plot of two vectors, say X and Y. plot(X,Y, col=blue) abline(a=1,b=1, col=red) abline(a=-1,b=1, col=green) I would like to split the scatter plot into 3 part with 3 different colors: (i)

[R] Siegel-Tukey Test

2007-04-27 Thread Ralf Finne
Hi everybody! Can anybody help me to find an R-function for the Siegel-Tukey test. Than you in advance. Ralf Finne SYH, University of Applied Sciences Vasa Finland __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] Integrating R-programs into larger systems

2007-04-27 Thread Ralf Finne
Hi experts! Have anybody experience in including an R-program as part of a larger system? In Matlab there is a toolbox that converts a m-script into C-code. One application in mind is that I do the model building in R, for estimating the risk for cancer based on clinical measurements. When

Re: [R] Unwanted white borders on semi-transparent polygons?

2007-04-27 Thread Matthew Neilson
Hey again, Just to let you know that I managed to test out R-2.5.0, and the transparency glitch has thankfully disappeared. I apologise for wasting your time with this. One final question, though. Did you mean to say that the known transparency bug was fixed in version 2.5.0? Or was it fixed

[R] weight

2007-04-27 Thread Natalie O'Toole
Hi, I have the file below called happyguys. It is a subset of data. How do I apply the weight variable (WTPP) to this file? Can i just multiply each column (except the first column because it is a record id) by WTPP? If the answer is yes, how do I multiply one variable name by another?

Re: [R] xyplot() and controlling panel.polygon()

2007-04-27 Thread Deepayan Sarkar
On 4/26/07, Michael Kubovy [EMAIL PROTECTED] wrote: I'm not sure I understand. I'm trying to add a filled band to show CIs around the data (which is why I tried to use xYplot, where 'filled bands' have bug). So with these data: est - c(1:4, 3:6, 7, 9, 11, 13, 12, 15, 18, 21) cond -

Re: [R] Unwanted white borders on semi-transparent polygons?

2007-04-27 Thread Prof Brian Ripley
On Fri, 27 Apr 2007, Matthew Neilson wrote: Hey again, Just to let you know that I managed to test out R-2.5.0, and the transparency glitch has thankfully disappeared. I apologise for wasting your time with this. One final question, though. Did you mean to say that the known transparency

Re: [R] add arrows to barchart with groups

2007-04-27 Thread Deepayan Sarkar
On 4/27/07, GOUACHE David [EMAIL PROTECTED] wrote: Hello Rhelpers, I am trying to represent the following data (hereafter named donnees) in a barchart with a grouping variable : sitetraitement dates res 1 NT 17/10/2005 normal 76.2 1 T

[R] [R-pkgs] New packages: contrast and desirability

2007-04-27 Thread Kuhn, Max
The contrast and desirability packages are available for all platforms at cran.r-project.org (and coming soon to a mirror near you). The contrast package extends Frank Harrell's contrast.Design function for one degree of freedom contrasts of model parameters to other types of models, such as lm,

Re: [R] Integrating R-programs into larger systems

2007-04-27 Thread Frank E Harrell Jr
Ralf Finne wrote: Hi experts! Have anybody experience in including an R-program as part of a larger system? In Matlab there is a toolbox that converts a m-script into C-code. One application in mind is that I do the model building in R, for estimating the risk for cancer based on

Re: [R] xyplot() and controlling panel.polygon()

2007-04-27 Thread Deepayan Sarkar
On 4/27/07, Michael Kubovy [EMAIL PROTECTED] wrote: Hi Deepayan, Thanks for your advice. This is moving along, however: The following is drawing the same polygons in each panel. I'm trying to get a different polygon (confidence band) for each group in each panel. That's why I thought I would

Re: [R] Integrating R-programs into larger systems

2007-04-27 Thread Rajarshi Guha
On Fri, 2007-04-27 at 20:29 +0300, Ralf Finne wrote: Hi experts! Have anybody experience in including an R-program as part of a larger system? In Matlab there is a toolbox that converts a m-script into C-code. One application in mind is that I do the model building in R, for estimating

[R] Protocol for data inclusion in new packages

2007-04-27 Thread Doran, Harold
In the near future I will release MiscPsycho, a package that contains various functions useful for applied psychometricians. I would like to include some data sets for distribution in the package, but have not created any of these on my own, but have used data distributed in other packages such as

Re: [R] Scatter Plot in R - Help

2007-04-27 Thread John Kane
Here is a recent posting by Petr Klasterecky that does not seem to be on the archive yet that may help. -- What do you mean by background? Maybe this is enough: plot(seq(-3,3,.01), dnorm(seq(-3,3,.01)), type=n, xlab=x, ylab=f(x), main=Normal density)

Re: [R] xyplot() and controlling panel.polygon()

2007-04-27 Thread Michael Kubovy
Hi Deepayan, Thanks for your advice. This is moving along, however: The following is drawing the same polygons in each panel. I'm trying to get a different polygon (confidence band) for each group in each panel. That's why I thought I would need to pass groups and subscripts to the

[R] Quadratcount() plotting in R spatstat

2007-04-27 Thread beamer
Hello, I am trying to plot a quadratcount object over a ppp object in the spatstat package. I ultimately want to get something like this http://bg9.imslab.co.jp/Rhelp/R-2.4.0/src/library/spatstat/man/images/big_quadratcount_001.png

[R] digest readability/ thunderbird email client/ R-help

2007-04-27 Thread Sam McClatchie
System: Linux kernel 2.6.15 Ubuntu dapper R version 2.4.1 ESS on Emacs 21.4.1 - Colleagues Has anyone figured out how to make the R-help digest more easily readable in the Thunderbird mail client? I think the digest used to be more readable than is currently the case

Re: [R] xyplot() and controlling panel.polygon()

2007-04-27 Thread Michael Kubovy
Hi Deepayan, Your solution works, anf the polygon are drawn where I wanted them to go. I thought that I could figure out how to gain control over the colors of the four ensuing polygons (I'm trying to get two lighter shades of the lines). I've tried, for example, to see if I could control

[R] Unsubscription Confirmation

2007-04-27 Thread Egoldsystem
Thank you for subscribing. You have now unsubscribed and no more messages will be sent. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] weight

2007-04-27 Thread Natalie O'Toole
__ Hi Dr. Kubovy, Here is my code so far: My question is: how do I then get a frequency count of Y_Q10A with the WTPP applied to it? myfile-(c:/test2.txt) mysubset-myfile mysubset$Y_Q02 -mysubset$DVSELF -NULL mysubset2-mysubset mysubset2$Y_Q10B -mysubset2$GP2_07 -NULL

Re: [R] bootstrap

2007-04-27 Thread Cody_Hamilton
Marc, I believe the problem is with the function you are passing to the statistic option. According to the documentation for boot, the statistic option provides: A function which when applied to data returns a vector

Re: [R] xyplot() and controlling panel.polygon()

2007-04-27 Thread Deepayan Sarkar
On 4/27/07, Michael Kubovy [EMAIL PROTECTED] wrote: Hi Deepayan, Your solution works, anf the polygon are drawn where I wanted them to go. I thought that I could figure out how to gain control over the colors of the four ensuing polygons (I'm trying to get two lighter shades of the lines).

[R] acf and pacf plot

2007-04-27 Thread tom soyer
Hi, I noticed that whenever I ran acf or pacf, the plot generated by R always includes two horizontal blue doted lines. Furthermore, these two lines are not documented in the acf documentation. I don't know what they are for, but it seems that they are important. Could someone tell me what they

[R] how to evaluate a dynamic variable?

2007-04-27 Thread Bobby Prill
Please help me evaluate a dynamic variable as show in the code snippet below. # regular variable a = c(1,2,3,4) # dynamic variable set at runtime DV = a eval(DV) -- a eval(eval(DV)) -- a # what I want something_goes_here(DV) -- 1,2,3,4 Can someone teach me how to do this? Thanks very

Re: [R] acf and pacf plot

2007-04-27 Thread Eric Thompson
The lines indicate the confidence interval (95% by default). I think you mean that it is not documented in help(acf), but it directs you to plot.acf in the See Also secion. From ?plot.acf: Note: The confidence interval plotted in 'plot.acf' is based on an _uncorrelated_ series and

Re: [R] xyplot() and controlling panel.polygon()

2007-04-27 Thread Michael Kubovy
I can't thank you enough! (I am also grateful to Sundar Dorai-Raj and Bert Gunter for their initial responses.) On Apr 27, 2007, at 4:14 PM, Deepayan Sarkar wrote: On 4/27/07, Michael Kubovy [EMAIL PROTECTED] wrote: Hi Deepayan, Your solution works, anf the polygon are drawn where I wanted

Re: [R] how to evaluate a dynamic variable?

2007-04-27 Thread Finny Kuruvilla
I'm not sure if I'm understanding what you're asking for but this may be it: DV=b assign(DV,1:4) DV [1] b b [1] 1 2 3 4 On Fri, 27 Apr 2007, Bobby Prill wrote: Please help me evaluate a dynamic variable as show in the code snippet below. # regular variable a = c(1,2,3,4) # dynamic

Re: [R] how to evaluate a dynamic variable?

2007-04-27 Thread Greg Snow
What are you trying to do with that? Most times when someone asks this question, they are trying to accomplish something that can be done better a different way, so if you tell us what you are trying to accomplish, we can suggest better approaches. The short answer to your question is ?get.

Re: [R] Scatter Plot in R - Help

2007-04-27 Thread Petr Klasterecky
If I understand correctly what A. Ezhil asked for, polygons won't help here. For coloring the individual points check ?points and use subsets, e.g. plot(x,y) points(x[xy-1],y[xy-1], col=red) or something similar. Petr John Kane napsal(a): Here is a recent posting by Petr Klasterecky that

Re: [R] how to evaluate a dynamic variable?

2007-04-27 Thread Bobby Prill
Mark, You're right. Thank you. Look, it works: a = c(1,2,3) a [1] 1 2 3 b = a b [1] a get(b) [1] 1 2 3 On Apr 27, 2007, at 5:05 PM, Leeds, Mark (IED) wrote: check out get using ?get. I'm not an expert but that might help or work. -Original Message- From: [EMAIL

Re: [R] weight

2007-04-27 Thread Natalie O'Toole
Does anyone know why it is giving me this error? Any help would be greatly appreciated!! Thanks, Nat myfile-(c:/test2.txt) mysubset-myfile mysubset$Y_Q02 -mysubset$DVSELF -NULL mysubset2-mysubset mysubset2$Y_Q10B -mysubset2$GP2_07 -NULL

Re: [R] how to evaluate a dynamic variable?

2007-04-27 Thread Don MacQueen
Does get(DV) do what you want? -Don At 4:35 PM -0400 4/27/07, Bobby Prill wrote: Please help me evaluate a dynamic variable as show in the code snippet below. # regular variable a = c(1,2,3,4) # dynamic variable set at runtime DV = a eval(DV) -- a eval(eval(DV)) -- a # what I want

[R] Restructuring Hadley CET data

2007-04-27 Thread Ted Harding
Hi Folks, I have a nasty data restructuring problem! I can think of one or two really clumsy ways of doing it with 'for' loops and the like, but I can't think of a *neat* way of doing it in R. The data are the Hadley Centre Central England Temperature series, daily from 01/01/1772 to 31/03/2007,

Re: [R] acf and pacf plot

2007-04-27 Thread tom soyer
Thanks Eric! I also noticed that in R, acf returns ac at lag 0, while pacf does not (pac for pacf starts at lag 1). Do you know if there is a good reason for that? Shouldn't ac at lag 0 always be 1? On 4/27/07, Eric Thompson [EMAIL PROTECTED] wrote: The lines indicate the confidence interval

Re: [R] Scatter Plot in R - Help

2007-04-27 Thread John Kane
When I reread it, I think you are right. I read split the data plot and assumed a backround rather than just the points . --- Petr Klasterecky [EMAIL PROTECTED] wrote: If I understand correctly what A. Ezhil asked for, polygons won't help here. For coloring the individual points check

Re: [R] Restructuring Hadley CET data

2007-04-27 Thread hadley wickham
Hi Ted, melt(df, id=c(year,DoM)) should get you started. Hadley On 4/27/07, Ted Harding [EMAIL PROTECTED] wrote: Hi Folks, I have a nasty data restructuring problem! I can think of one or two really clumsy ways of doing it with 'for' loops and the like, but I can't think of a *neat* way

Re: [R] Restructuring Hadley CET data

2007-04-27 Thread hadley wickham
Ooops, add library(reshape) before that, and see http://had.co.nz/reshape for more info Hadley On 4/27/07, hadley wickham [EMAIL PROTECTED] wrote: Hi Ted, melt(df, id=c(year,DoM)) should get you started. Hadley On 4/27/07, Ted Harding [EMAIL PROTECTED] wrote: Hi Folks, I have a nasty

[R] xy plot key colors don't match plot

2007-04-27 Thread Ken Nussear
Hi Trying to finish a simple xyplot, but the key colors don't match those specified on the plot. I have only been able to change the color of the key text, but not the points. Hours - c (25.88,32.31,24.04,23.27,16.64,19.75,27.25,17.9,20.45,29.65,20.83,27.17,

[R] Limit on vector evaluation

2007-04-27 Thread Robert Barber
Dear R Experts, Why I try to run this expression: x-sapply(rnorm(rep(10,10),mean=9,sd=1.5),mean) it evaluates the first 1 values and then stops, but does not return to the command prompt. My cpu keeps running at 100%. When I exit the expression with CTL-C, I then see that x holds

Re: [R] Restructuring Hadley CET data

2007-04-27 Thread Ted Harding
On 27-Apr-07 22:43:58, hadley wickham wrote: Ooops, add library(reshape) before that, and see http://had.co.nz/reshape for more info Hadley On 4/27/07, hadley wickham [EMAIL PROTECTED] wrote: Hi Ted, melt(df, id=c(year,DoM)) should get you started. Hadley On 4/27/07, Ted Harding

Re: [R] xy plot key colors don't match plot

2007-04-27 Thread Deepayan Sarkar
On 4/27/07, Ken Nussear [EMAIL PROTECTED] wrote: Hi Trying to finish a simple xyplot, but the key colors don't match those specified on the plot. I have only been able to change the color of the key text, but not the points. Hours - c

Re: [R] Restructuring Hadley CET data

2007-04-27 Thread Gabor Grothendieck
Here is a bit more. It creates a zoo series. URL - http://hadobs.metoffice.com/hadcet/cetdl1772on.dat; DF - read.table(URL, col.names = c(Year, DoM, 1:12), check = FALSE) library(reshape) long - melt(DF, c(Year, DoM)) long - subset(long, value != -999) library(zoo) dd - as.Date(with(long,

Re: [R] xy plot key colors don't match plot

2007-04-27 Thread Ken Nussear
Deepayan fixed it in the end with this command. Posting for future questions. Thanks! Ken xyplot(DistanceMeters~Hours, groups=Surveyor, xlab=Time (h), ylab=Distace Moved(m), par.settings = list(superpose.symbol = list ( pch = c(1, 16), col=c(black,black))), auto.key=list(space =

[R] like SPSS

2007-04-27 Thread Natalie O'Toole
Hi, I've written code to extact a pumf file in R, subset it, and weight it like you would do in SPSS. My code is below it works great. My question is: how do i then calculate the frequencies of smokers (1) versus non-smokers (2) after having weighted my file? or even the process that SPSS is

Re: [R] Limit on vector evaluation

2007-04-27 Thread Duncan Murdoch
On 4/27/2007 7:13 PM, Robert Barber wrote: Dear R Experts, Why I try to run this expression: x-sapply(rnorm(rep(10,10),mean=9,sd=1.5),mean) it evaluates the first 1 values and then stops, but does not return to the command prompt. My cpu keeps running at 100%. When I exit the

[R] Configuring R-2.5.0 for X11 on Red Hat WS 4 x86-64

2007-04-27 Thread Bruce Foster
I have been unable to configure R-2.5.0 for the X11 interface on Red Hat WS 4 x86-64 (AMD Opteron). The configure script looks for XtToolkitInitialize in libXt, and does not find it because it's not in libXt (which is present, however). I tried installing openmotif and openmotif-devel, and

  1   2   >