Re: [R] na.omit - Is it working properly?

2011-05-04 Thread peter dalgaard
On May 3, 2011, at 21:18 , Kalicin, Sarah wrote: I have a work around for this, but can someone explain why the first example does not work properly? I believed it worked in the previous version of R, by selecting just the rows=200525 and omitting the na's. I just upgraded to 2.13. I am

[R] problem with package adapt for R in Mac

2011-05-04 Thread Matías Ramírez Salgado
Hi, How i can install the package adapt in some version of R for mac? i try in 2.13, 2.9,2.7 and other previous versions... and nothing happens. and another question: There are some packages that do the same but that it is implemented for mac? (calculate integrals in 2 or more dimmensions).

[R] nls problem with R

2011-05-04 Thread sterlesser
the original data are V2 =c(371000,285000 ,156000, 20600, 4420, 3870, 5500 ) T2=c( 0.3403 ,0.4181 ,0.4986 ,0.7451 ,1.0069 ,1.553) nls2=nls(V2~v0*(1-epi+epi*exp(-cl*(T2-t0))),start=list(v0=10^7,epi=0.9,cl=6.2,t0=8.7)) after execution error occurs as below

Re: [R] How to fit a random data into Beta distribution?

2011-05-04 Thread Shekhar
Hi Steven, Thanks for the quick reply. i have tried but its giving me error---Error in optim(x = c(38.1815173696765, -12.7988197976440, -3.88212459045077, : initial value in 'vmmin' is not finite i have tried something like this: library(MASS) x-rnorm(n=100,mean=10,sd=20);

Re: [R] RStudio -manipulate command

2011-05-04 Thread veepsirtt
Why the mean value h is not changing as the slider moves from 0 to 25 ?. It remains always constant. library(manipulate) example - function(x.max){ plot(cars, xlim=c(0,x.max)) abline(h=mean(cars$dist),col=blue,lty=2) } manipulate( example(x.max), x.max=slider(0,25, step=5) ) veepsirtt --

[R] HELP

2011-05-04 Thread Nekeisha Spencer
Accelerated Failure Time Model to Proportional Hazard Form Greetings R users: I have been working on a problem for a while and can't seem to get any result. I am trying to convert accelerated failure time estimates to proportional form. I keep getting an error that I can't understand and

Re: [R] help with the maxBHHH routine

2011-05-04 Thread Arne Henningsen
Dear Rohit On 3 May 2011 22:53, Rohit Pandey rohitpandey...@gmail.com wrote: Hello R community, I have been using R's inbuilt maximum likelihood functions, for the different methods (NR, BFGS, etc). I have figured out how to use all of them except the maxBHHH function. This one is

Re: [R] nls problem with R

2011-05-04 Thread Andrew Robinson
The fact that T2 and V2 are of different lengths seems like a likely culprit. Other than that, you need to find start points that do not lead to a singular gradient. There are several books that provide advice on obtaining initial parameter estimates for non-linear models. Google Books might

Re: [R] problem with package adapt for R in Mac

2011-05-04 Thread Andrew Robinson
Hi, Is there such a package? I can't find it on CRAN. Can you let us know exactly how you tried to install it, and what the error message was (if any)? Cheers Andrew On Wed, May 04, 2011 at 01:29:37AM -0300, Mat?as Ram?rez Salgado wrote: Hi, How i can install the package adapt in some

Re: [R] Constructing a histogram with words as labels as height as frequency?

2011-05-04 Thread Jim Lemon
On 05/04/2011 03:35 AM, Caitlin wrote: Hi all. I need to construct a plot showing words on the x-axis and how many times each word was given as a verbal response on the y-axis as solid bar (frequency). Is there a convenient function to do this in R? I considered hist(), but I'm not sure how to

[R] adding columns to dataframes contained in a list

2011-05-04 Thread E Hofstadler
hi there, I have a list of 5 identical dataframes: mydf - data.frame(x=c(1:5), y=c(21:25)) mylist - rep(list(mydf),5) and a factor variable with 5 levels: foo - c(letters[1:5]) foo - as.factor(foo) Question: I'd like to add a new variable to each dataframe in the list, each containing only

Re: [R] nls problem with R

2011-05-04 Thread Rubén Roa
In addition to Andrew's advice, you should get more familiar with your nonlinear model. From what you wrote, as T2 tends to infinity, V2 tends to v0*(1-epi). There you have a baseline on the Y-axis towards which your model tends, and this will give you sensible starting values for v0 and epi.

Re: [R] Simple loop

2011-05-04 Thread Petr Savicky
On Tue, May 03, 2011 at 12:04:47PM -0700, William Dunlap wrote: [...] ave() can deal that problem: cbind(x, newCol2 = with(x, ave(H, Site, Prof, FUN=function(y)y-min(y Site Prof H newCol2 111 24 8 211 16 0 311 67 51 412

Re: [R] Simple General Statistics and R question (with 3 line example) - get z value from pairwise.wilcox.test

2011-05-04 Thread JP
On 3 May 2011 20:50, peter dalgaard pda...@gmail.com wrote: On Apr 28, 2011, at 15:18 , JP wrote: I have found that when doing a wilcoxon signed ranked test you should report: - The median value (and not the mean or sd, presumably because of the underlying potential non normal

[R] Read last line of a file

2011-05-04 Thread Joel
Hi dear R users. I got a file that I need to extract the third column (or word) of the last line of files that has a diffrent amounts of rows. It works with x-read.tables(file) x[1,3] This returns the proper result but as the files is large this takes time and uses memory that is just

Re: [R] problem with package adapt for R in Mac

2011-05-04 Thread Uwe Ligges
On 04.05.2011 09:19, Andrew Robinson wrote: Hi, Is there such a package? There was such a package that is archived now. I can't find it on CRAN. Can you let us know exactly how you tried to install it, and what the error message was (if any)? The OP is probably looking for packages

Re: [R] Read last line of a file

2011-05-04 Thread Petr Savicky
On Wed, May 04, 2011 at 02:17:29AM -0700, Joel wrote: Hi dear R users. I got a file that I need to extract the third column (or word) of the last line of files that has a diffrent amounts of rows. It works with x-read.tables(file) x[1,3] This returns the proper result but as the

Re: [R] How to fit a random data into Beta distribution?

2011-05-04 Thread Steven Kennedy
Hi Shekhar, It looks from your error that you have values outside the range 0 to 1. The beta distribution is only defined between 0 and 1. Can you please post your data set (or at least a portion of it)? On Wed, May 4, 2011 at 3:24 PM, Shekhar shekhar2...@gmail.com wrote: Hi Steven, Thanks

[R] issue with strange characters (locale settings)

2011-05-04 Thread R.T.A.J.Leenders
WinXP-x32, R-21.13.0 Dear list, I have a problem that (I think) relates to the interaction between Windows and R. I am trying to scrape a table with data on the Hawai'ian Islands, This is my code: library(XML) u - http://en.wikipedia.org/wiki/Hawaii; tables -

Re: [R] adding columns to dataframes contained in a list

2011-05-04 Thread Patrick Connolly
On Wed, 04-May-2011 at 10:25AM +0300, E Hofstadler wrote: | hi there, | | I have a list of 5 identical dataframes: | | mydf - data.frame(x=c(1:5), y=c(21:25)) | mylist - rep(list(mydf),5) | | and a factor variable with 5 levels: | | foo - c(letters[1:5]) | foo - as.factor(foo) | | |

Re: [R] bootstrap vignette construction and package installation

2011-05-04 Thread Uwe Ligges
On 03.05.2011 22:20, Murat Tasan wrote: hi all - i'm trying to 'R CMD build' a package, but i have what appears to be a bootstrapping problem: i've included a vignette in my package, with R code interwoven (and built using Sweave), but in this documentation i have a code line:

Re: [R] issue with strange characters (locale settings)

2011-05-04 Thread Prof Brian Ripley
Oh, please! This is about the contributed package XML, not R and not Windows. Some of us have worked very hard to provide reasonable font support in R, including on Windows. We are given exceedingly little credit, just the brickbats for things for which we are not responsible. (We even work

Re: [R] Overlapping x axes using Lattice

2011-05-04 Thread Patrick Breheny
The scales argument in lattice controls the appearance of the axes. It consists of two lists, one for the x axis and one for the y axis. For example: histogram(~pesti[,1]|pesti[,2]+ pesti[,3],scales=list(x=list(at=c(2,6),labels=c(First,Second This allows you to place the labels

[R] what happens when I store linear models in an array?

2011-05-04 Thread Andrew D. Steen
I've got a bunch of similar datasets, all of which I've fit to linear models. I'd like to easily create arrays of a specific parameter from each linear model (e.g., all of the intercepts in one array). I figured I'd put the model objects into an array, and then (somehow) I could easily create

Re: [R] Bigining with a Program of SVR

2011-05-04 Thread Max Kuhn
train() uses vectors, matrices and data frames as input. I really think you need to read materials on basic R before proceeding. Go to the R web page. There are introductory materials there. On Tue, May 3, 2011 at 11:19 AM, ypriverol yprive...@gmail.com wrote: I saw the format of the caret data

[R] first occurrence of a value?

2011-05-04 Thread Albert-Jan Roskam
Hello, A simple question perhaps, but how do I, within each row, find the first occurence of the number 1 in the df below? I want to use this position to programmatically create the variable 'year'. I'v come up with a solution, but I find it downright ugly. Is there a simpler way? I was hoping

Re: [R] first occurrence of a value?

2011-05-04 Thread Dimitris Rizopoulos
another approach is: df - data.frame(j1999 = c(0,0,0,0,1,0), j2000 = c(NA,1,1,1,0,0), j2001 = c(1,0,1,0,0,0)) years - as.numeric(gsub(^[^0-9]+, , names(df))) ind - apply(sapply(df, ==, 1), 1, function (x) which(x)[1]) df$year - years[ind] I hope it helps. Best, Dimitris On 5/4/2011

Re: [R] first occurrence of a value?

2011-05-04 Thread Patrick Breheny
You may want to look into the function 'match', which finds the first occurrence of a value. In your example, df - data.frame(j1999=c(0,0,0,0,1,0), j2000=c(NA, 1, 1, 1, 0, 0), j2001=c(1, 0, 1, 0, 0, 0), year=c(2001, 2000, 2000, 2000, 1999, NA)) apply(df,1,match,x=1) [1] 3 2 2 2 1 NA

Re: [R] Simple General Statistics and R question (with 3 line example) - get z value from pairwise.wilcox.test

2011-05-04 Thread peter dalgaard
On May 4, 2011, at 11:03 , JP wrote: On 3 May 2011 20:50, peter dalgaard pda...@gmail.com wrote: On Apr 28, 2011, at 15:18 , JP wrote: I have found that when doing a wilcoxon signed ranked test you should report: - The median value (and not the mean or sd, presumably because of

Re: [R] first occurrence of a value?

2011-05-04 Thread Albert-Jan Roskam
Hi Patrick, Dimitri, Thank you! Yes, 'match' was exactly what I was looking for. I like it as it doesn't require too many functions to be nested. Cheers!! Albert-Jan ~~ All right, but apart from the sanitation, the medicine,

Re: [R] compute coefficient of determination (R-squared) for GLM (maximum likelihood)

2011-05-04 Thread KateR
Dear mr Joris Meys I would like to know, where to find this formulas in books or articles. # possibility 1 R2 - cor(y,predict(mod))^2 # possibility 2 R2 - 1 - ( sum( (y-predict(mod))^2 ) / sum( (y-mean(y))^2 ) ) The first calculation seems OK, it gives the logicala values in models

Re: [R] How to fit a random data into Beta distribution?

2011-05-04 Thread Ali Akbar Bromideh
Dear Shekar, For me I had such a problem, when I repeated with many dist. So, I drafted my own procedure for example for Weibull, Gamma and Log-Normal without any warnings. But you can ignore the Warning, since the result was somehow the same for my case. Regards, Ali Date: Wed, 4 May 2011

[R] select value from a column depending on a value in another column

2011-05-04 Thread tornanddesperate
Hi everybody I couldn't find the solution to what must be quite a simple problem. Maybe you can help? treatment session period stage wage_accepted market 1 11 11 25 public 2 11 11 19

[R] Format ddmmYYYY in date

2011-05-04 Thread smoff
Hello everybody, I'm quite new in using R so please do not kill me if I ask stupid questions. My problem is that I have a table containing dates in the first column of 10 years. These dates have the format ddmm at least in the csv-file. After importing the file using read.table() R deletes

[R] Is this confict of different versions of R or something else?

2011-05-04 Thread TheSavageSam
Hello! I have had some problems lately with use of R at home and school. At my home laptop (Ubuntu linux 64bit R-2.12.0) R works just fine. But when I take my codes to school(Windows XP 32bit R-2.10.1 I think) and run those there those codes probably won't work. Functions as combinations()

Re: [R] ID parameter in model

2011-05-04 Thread Mike Harwood
Thank you, Goran. Please see the package details below: packageDescription('eha') Encoding: UTF-8 Package: eha Version: 1.3-2 Date: 2011-03-01 Title: Event History Analysis Description: A package for survival and event history analysis License: GPL (= 3) Author: Göran Broström Depends: R (=

[R] bivariate linear interpolation

2011-05-04 Thread Halldór Björnsson
Hi, I have three matrices (X,Y,P) with the same dimension. The X,Y grid is regular and I want to perform linear interpolation to pick out certain points. In matlab appropriate call is something like Pout=interp2(X,Y,P,Xout,Yout, method=linear) where Xout and Yout are the locations where I want

[R] scatter plot with Z value

2011-05-04 Thread Matevž Pavlič
Hi all, I would like to create a scatter plot of two variables (y, x) whith third value (z) written on the plot? After that i would like to add a line (Y=0.7*(x-20)) to the graph. I tried plot(x~y) but there is no command for the third vairable to be shown on the graph also i

Re: [R] Format ddmmYYYY in date

2011-05-04 Thread Ben Bolker
smoff simon.frey at boku.ac.at writes: My problem is that I have a table containing dates in the first column of 10 years. These dates have the format ddmm at least in the csv-file. After importing the file using read.table() R deletes the first character if it is a zero. [snip] How

Re: [R] select value from a column depending on a value in another column

2011-05-04 Thread Andrew Robinson
Try subset(). Andrew On Wed, May 04, 2011 at 02:01:44AM -0700, tornanddesperate wrote: Hi everybody I couldn't find the solution to what must be quite a simple problem. Maybe you can help? treatment session period stage wage_accepted market 1 11 1

Re: [R] what happens when I store linear models in an array?

2011-05-04 Thread Andrew Robinson
Hi Andrew, try fitted(lms.ASP[1,1][[1]]) Cheers Andrew On Wed, May 04, 2011 at 01:49:45PM +0200, Andrew D. Steen wrote: I've got a bunch of similar datasets, all of which I've fit to linear models. I'd like to easily create arrays of a specific parameter from each linear model (e.g., all

Re: [R] Simple General Statistics and R question (with 3 line example) - get z value from pairwise.wilcox.test

2011-05-04 Thread JP
Peter thanks for the fantastically simple and understandable explanation... To sum it up... to find the z values of a number of pairwise wilcox tests do the following: # pairwise tests with bonferroni correction x - pairwise.wilcox.test(a, b, alternative=two.sided, p.adj=bonferroni, exact=F,

Re: [R] scatter plot with Z value

2011-05-04 Thread Jannis
I am sure you would have found answers to your questions if you would have searched the mailing list archive (http://r.789695.n4.nabble.com/R-help-f789696.html)! To get you started, have a look at: ?text (for the z values) ?abline (for the line) Jannis --- Matevž Pavlič

[R] Problems saving ff objects

2011-05-04 Thread Jannis
Dear list, I am trying to understand and use the ff package. As I had some problems saving some ff objects, and as I did not fully manage to understand the whole concept of *.ff, *.ffData and *.RData with the help of the documentation, I tried to reproduce the examples from the help of

Re: [R] Watts Strogatz game

2011-05-04 Thread Robert Baer
I have a erdos-renyi game with 6000 nodes and probability 0.003. g1 = erdos.renyi.game(6000, 0.003) How to create a Watts Strogatz game with the same probability. g1 = watts.strogatz.game(1, 6000, ?, ?) What should be the third and fourth parameter to this argument. According to

Re: [R] Problems saving ff objects

2011-05-04 Thread Jannis
Just did some more testing.May the problem be due to the fact that I am using a windows machine? I just ran the same code on a Linux machine and everything worked fine. If windows (or the file system of the disk) caused the problem, is there any way to resolve it? I know that using Linux

Re: [R] what happens when I store linear models in an array?

2011-05-04 Thread Jonathan Daily
It looks like your call: lms.ASP[1,1] - list(surf105.lm.ASP) makes a list of length 1 containing the lm object and puts that list into element [1,1] of your array. That is why you will need the extra indexing call of [[1]] Andrew Robinson suggested. On Wed, May 4, 2011 at 9:02 AM, Andrew

Re: [R] select value from a column depending on a value in another column

2011-05-04 Thread Scott Chamberlain
This for example could work df - data.frame(x = seq(1:10), y = rep(c(a,b),5)) df df$x[df$y %in% a] Scott On Wednesday, May 4, 2011 at 4:01 AM, tornanddesperate wrote: Hi everybody I couldn't find the solution to what must be quite a simple problem. Maybe you can help? treatment

[R] Storing data from a test as a vector or matrix

2011-05-04 Thread wwreith
I just finished a MANOVA test and got the following output: summary(M, test=Pillai) Df Pillai approx F num Df den DfPr(F) as.factor(X) 3 1.1922 6.5948 36360 2.2e-16 *** Residuals 129 --- Signif. codes: 0 '***'

[R] Superscript number before letter

2011-05-04 Thread Janhal
Salut, I have been struggling to superscript the 18 before the O without ^ visible and found only help to superscript numbers after the letter. Thanks to anyone who can help. xlab=expression(delta*18O VSMOW [‰]) Cheers, Janine -- View this message in context:

Re: [R] Unexp. behavior from boot with multiple statistics

2011-05-04 Thread algorimancer
Thanks, that clears things up quite a bit. Now I'm left wondering why there is so much bias, but that's a separate issue. -- View this message in context: http://r.789695.n4.nabble.com/Unexp-behavior-from-boot-with-multiple-statistics-tp3493300p3495590.html Sent from the R help mailing list

[R] create a folder with mode '0777'

2011-05-04 Thread Xian Zhang
Dear list, I am trying to create a folder structure, say 'test/sub', and set the folder and sub folder to be writable to everyone. By default dir.create('test/sub', recursive=TRUE, mode='0777') creates folders with mode: drwxr-xr-x After Sys.chmod('test/sub',mode='0777') The folder 'test'

Re: [R] nls problem with R

2011-05-04 Thread sterlesser
Thanks Ruben. Your suggestion about more deeper analysis about the model itself is really helpful. I am trying out some new initial values based on the analysis of the special T2 in the model. -- View this message in context:

Re: [R] ID parameter in model

2011-05-04 Thread Göran Broström
On Wed, May 4, 2011 at 2:18 PM, Mike Harwood maharw...@hotmail.com wrote: Thank you, Goran.  Please see the package details below: Thanks, I have uploaded a corrected version of eha to CRAN. Should be available soon. Göran [...] __

Re: [R] nls problem with R

2011-05-04 Thread sterlesser
Thanks Andrew. I am sorry for some typos that I omit some numbers of T2. Based on your suggestion,I think the problem is in the initial values. And I will read more theory about the non-linear regression. -- View this message in context:

Re: [R] Watts Strogatz game

2011-05-04 Thread Gábor Csárdi
On Wed, May 4, 2011 at 9:28 AM, Robert Baer rb...@atsu.edu wrote: I have a erdos-renyi game with 6000 nodes and probability 0.003. g1 = erdos.renyi.game(6000, 0.003) How to create a Watts Strogatz game with the same probability. g1 = watts.strogatz.game(1, 6000, ?, ?) What should be the

Re: [R] nls problem with R

2011-05-04 Thread Ravi Varadhan
In addition to the suggestion about finding a good initial value, you should also scale your response V2 (and, of course, V0). Divide V2 by 10^4, for example. Now your V0 should also be scaled by this factor. This would likely help with convergence. Ravi.

Re: [R] Superscript number before letter

2011-05-04 Thread Richard M. Heiberger
Dos this do what you want? plot(1:10, xlab=expression(delta*{}^18*O * VSMOW [‰])) The specific is to put an empty item there to hold the superscript. On Wed, May 4, 2011 at 9:37 AM, Janhal janine.hal...@unil.ch wrote: Salut, I have been struggling to superscript the 18 before the O

[R] Odp: Is this confict of different versions of R or something else?

2011-05-04 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 04.05.2011 12:44:04: TheSavageSam jujok...@jyu.fi Odeslal: r-help-boun...@r-project.org Hello! I have had some problems lately with use of R at home and school. At my home laptop (Ubuntu linux 64bit R-2.12.0) R works just fine. But when I

Re: [R] Simple General Statistics and R question (with 3 line example) - get z value from pairwise.wilcox.test

2011-05-04 Thread peter dalgaard
On May 4, 2011, at 15:11 , JP wrote: Peter thanks for the fantastically simple and understandable explanation... To sum it up... to find the z values of a number of pairwise wilcox tests do the following: # pairwise tests with bonferroni correction x - pairwise.wilcox.test(a, b,

Re: [R] create a folder with mode '0777'

2011-05-04 Thread Patrick Breheny
Linux systems have a user mask that limits the file mode creation possibilities of any processes launched from that shell. If you check your /etc/profile file, you will see the line umask 022 This prevents you by default from creating files with write access for everyone except the user.

Re: [R] Format ddmmYYYY in date

2011-05-04 Thread David Winsemius
On May 4, 2011, at 5:57 AM, Ben Bolker wrote: smoff simon.frey at boku.ac.at writes: My problem is that I have a table containing dates in the first column of 10 years. These dates have the format ddmm at least in the csv- file. After importing the file using read.table() R deletes the

Re: [R] what happens when I store linear models in an array?

2011-05-04 Thread David Winsemius
On May 4, 2011, at 4:49 AM, Andrew D. Steen wrote: I've got a bunch of similar datasets, all of which I've fit to linear models. I'd like to easily create arrays of a specific parameter from each linear model (e.g., all of the intercepts in one array). I figured I'd put the model objects

Re: [R] Format ddmmYYYY in date

2011-05-04 Thread David Winsemius
On May 4, 2011, at 5:57 AM, Ben Bolker wrote: smoff simon.frey at boku.ac.at writes: My problem is that I have a table containing dates in the first column of 10 years. These dates have the format ddmm at least in the csv- file. After importing the file using read.table() R deletes the

Re: [R] Storing data from a test as a vector or matrix

2011-05-04 Thread Ivan Calandra
Hi, I would suggest you to check the structure of your summary object with str(), like this: S - summary(M, test=Pillai) str(S) You will then see how to access each element of it. If you cannot manage to do it yourself, then provide an example, or at least the output from str(s). By the way,

Re: [R] Superscript number before letter

2011-05-04 Thread David Winsemius
On May 4, 2011, at 7:28 AM, Richard M. Heiberger wrote: Dos this do what you want? plot(1:10, xlab=expression(delta*{}^18*O * VSMOW [‰])) The specific is to put an empty item there to hold the superscript. I do not think that is necessary: plot(1:10, xlab=expression(delta^18*O~VSMOW[‰]))

Re: [R] Is this confict of different versions of R or something else?

2011-05-04 Thread David Winsemius
On May 4, 2011, at 3:44 AM, TheSavageSam wrote: Hello! I have had some problems lately with use of R at home and school. At my home laptop (Ubuntu linux 64bit R-2.12.0) R works just fine. But when I take my codes to school(Windows XP 32bit R-2.10.1 I think) and run those there those

Re: [R] create a folder with mode '0777'

2011-05-04 Thread Prof Brian Ripley
Please read the comments in the help about umask (and in the posting guide about the 'at a minimum' information required in postings, for the details do depend on the version of R and it seems yours is not current). In R 2.13.0: ‘dir.create’ creates the last element of the path, unless

[R] hurdle, simulated power

2011-05-04 Thread David Atkins
Hi all-- We are planning an intervention study for adolescent alcohol use, and I am planning to use simulations based on a hurdle model (using the hurdle() function in package pscl) for sample size estimation. The simulation code and power code are below -- note that at the moment the

[R] natural cubic splines

2011-05-04 Thread Ondřej Mikula
Dear R-helpers, I need to fit natural cubic spline with specified number of knots. I expected 'splines' package will be helpful, but I am confused by its help. Is more detailed documentation available for it or could you recommend another R function? Best regards Ondrej Mikula

Re: [R] natural cubic splines

2011-05-04 Thread Joshua Wiley
Hi Ondrej, What documentation have you looked at? Does this help at all? require(splines) ?ns ## one example summary(lm(y ~ ns(x, df = 3), data = data.frame(y = runif(100), x = rbinom(100, 9, .25)^2))) ## built in examples example(ns) Also, I am very fond of the book, Modern Applied

[R] Error Rscript: No such file or directory

2011-05-04 Thread xavier abulker
Hello, I'm trying to build a simple cpp file using the R CMD SHLIB command and I always receive the same error message: cygwin warning: MS-DOS style path detected: C:/PROGRA~1/R/R-212~1.1/etc/i386/Makeconf Preferred POSIX equivalent is: /cygdrive/c/PROGRA~1/R/R-212~1.1/etc/i386/Makeconf

Re: [R] Superscript number before letter

2011-05-04 Thread Richard M. Heiberger
David, That is not clear from the original request. The request was for {}^18*O It wasn't for delta^18 Therefore I put the space there to be sure that the 18 was seen as pre-superscript of O, not as a post-superscript of delta. I probably should also have used ~ as plot(1:10,

Re: [R] Simple loop

2011-05-04 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Petr Savicky Sent: Wednesday, May 04, 2011 12:51 AM To: r-help@r-project.org Subject: Re: [R] Simple loop On Tue, May 03, 2011 at 12:04:47PM -0700, William Dunlap wrote:

[R] Panels order in lattice graphs

2011-05-04 Thread Cristina Silva
Hi all, In lattice graphs, panels are drawn from left to right and bottom to top. The flag as.table=TRUE changes to left to right and top to bottom. Is there any way to change to first top to bottom and then left to right? didn´t find anything neither in Help pages nor Lattice book.

[R] xtable without a loop alongside a ggplot

2011-05-04 Thread Justin Haynes
I would like to create a table of my points and identify which 'quadrant' of a plot they are in with the 'origin' at the means. the kicker is i would like to display it right next to or below a ggplot of the data. Maybe xtable isnt the right thing to use, but its the only thing i can think of.

Re: [R] Error Rscript: No such file or directory

2011-05-04 Thread Duncan Murdoch
On 04/05/2011 11:49 AM, xavier abulker wrote: Hello, I'm trying to build a simple cpp file using the R CMD SHLIB command and I always receive the same error message: cygwin warning: MS-DOS style path detected: C:/PROGRA~1/R/R-212~1.1/etc/i386/Makeconf Preferred POSIX equivalent is:

Re: [R] Error Rscript: No such file or directory

2011-05-04 Thread xavier abulker
Thanks a lot Duncan, I have removed the system variable incorrectly setup and it works now. Have a nice day Xavier From: Duncan Murdoch murdoch.dun...@gmail.com Cc: r-help@r-project.org Sent: Wed, 4 May, 2011 17:37:55 Subject: Re: [R] Error Rscript: No such

[R] tryCatch?

2011-05-04 Thread Mikkel Grum
I would like to do inserts into a database table, but do updates in the fairly rare cases in which the inserts fail. I thought tryCatch might be the way to do it, but I honestly do not understand the help file for tryCatch at all. I thought something like this might work: for (i in seq(along =

Re: [R] Simple loop

2011-05-04 Thread Petr Savicky
On Wed, May 04, 2011 at 08:52:07AM -0700, William Dunlap wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Petr Savicky Sent: Wednesday, May 04, 2011 12:51 AM To: r-help@r-project.org Subject: Re: [R] Simple loop

Re: [R] tryCatch?

2011-05-04 Thread Prof Brian Ripley
Start with try(): you may find it easier to understand. if(inherits(try(call1), try-error)) call2 so in your case if(inherits(try(sqlSave(pg, tbl[i, ], tbl, append = TRUE, rownames = FALSE sqlUpdate(pg, tbl[i, ], index = key) or some such. On Wed, 4 May 2011,

[R] problem with package adapt for R in Mac

2011-05-04 Thread Matías Ramírez Salgado
The package is called: adapt (adapt_1.0-4.tgz is the version of package for mac, search it in google) is not in all CRAN servers, because is an old package. I tried to install it from package manager as CRAN (binaries), Other repositorys, and as a local source packages, and allways returns this

Re: [R] problem with package adapt for R in Mac

2011-05-04 Thread Uwe Ligges
On 04.05.2011 19:22, Matías Ramírez Salgado wrote: The package is called: adapt (adapt_1.0-4.tgz is the version of package for mac, search it in google) is not in all CRAN servers, because is an old package. Yes, there is a reason why the package was archived (which I said already!): It does

Re: [R] Superscript number before letter

2011-05-04 Thread David Winsemius
On May 4, 2011, at 8:43 AM, Richard M. Heiberger wrote: David, That is not clear from the original request. The request was for {}^18*O I see your point. The use of the phantom would also allow a pre- superscript with no preceding delta. Thanks for bearing with my obtuseness. --

Re: [R] natural cubic splines

2011-05-04 Thread David Winsemius
On May 4, 2011, at 8:34 AM, Joshua Wiley wrote: Hi Ondrej, What documentation have you looked at? Does this help at all? require(splines) ?ns ## one example summary(lm(y ~ ns(x, df = 3), data = data.frame(y = runif(100), x = rbinom(100, 9, .25)^2))) ## built in examples example(ns) Also,

Re: [R] tryCatch?

2011-05-04 Thread Mikkel Grum
Beautiful Prof. This worked: for (i in seq(along = tbl$key)) { if (inherits(try(sqlSave(pg, tbl[i, ], tbl, append = TRUE, rownames = FALSE), silent = TRUE), try-error, TRUE)) sqlUpdate(pg, tbl[i, ], index = key1) ) } --- On Wed, 5/4/11, Prof Brian Ripley

[R] combine lattice plot and standard R plot

2011-05-04 Thread Lucia Cañas
Dear R users, I would like to combine lattice plot (xyplot) and standard R plot (plot and plotCI) in an unique figure. I use the function par() to combine plot and plotCI and I use the function print() to combine xyplot. I tried to use these functions to combine xyplot and plotCI and plots

Re: [R] combine lattice plot and standard R plot

2011-05-04 Thread Jonathan Daily
If you read the help documentation, lattice is not really compatible with standard graphics. library(lattice) ?lattice 2011/5/4 Lucia Cañas lucia.ca...@co.ieo.es: Dear R users, I would like to combine lattice plot (xyplot) and standard R plot (plot and plotCI) in an unique figure. I use

Re: [R] combine lattice plot and standard R plot

2011-05-04 Thread Scott Chamberlain
What about the example in gridExtra package: require(ggplot2); require(lattice); require(gridExtra) grid.arrange(qplot(1:10), xyplot(1:10~1:10), tableGrob(head(iris)), nrow=2, as.table=TRUE, main=test main, sub=textGrob(test sub, gp=gpar(font=2))) On Wednesday, May 4, 2011 at 1:44 PM, Jonathan

[R] Regexp question

2011-05-04 Thread johannes rara
I have a string like this st - SELECT COUNT(empid), COUNT(mgrid), COUNT(empname), COUNT(salary), FROM Employees How can I remove the last comma before the FROM statement? -J __ R-help@r-project.org mailing list

[R] merging multiple columns from two dataframes

2011-05-04 Thread 1Rnwb
Hello, I have data in a dataframe with 139104 rows which is multiple of 96x1449. i have a phenotype file which contains the phenotype information for the 96 samples. the snp name is repeated 1449X96 samples. I haveto merge the two dataframes based on sid and sen. this is how my two dataframes

[R] join tables in R

2011-05-04 Thread Roccato Alfredo (UniCredit)
I'd to match-merge 2 tables in such a manner that I keep all the rows in table 1, but not the rows that are in both table 1 and 2. Thank you for your help, Alfredo master - data.frame(ID=2001:2011) train - data.frame(ID=2004:2006) valid - ??? in this example table valid should have the

Re: [R] nls problem with R

2011-05-04 Thread Mike Marchywka
Date: Wed, 4 May 2011 07:07:44 -0700 From: sterles...@hotmail.com To: r-help@r-project.org Subject: Re: [R] nls problem with R Thanks Andrew. I am sorry for some typos that I omit some numbers of T2. Based on your suggestion,I think the problem is in the initial values. And I will read

[R] Saving Values in a Vector from a For Loop

2011-05-04 Thread blutack
Hi, I have a created a function, but now I need to call it about a hundred times and store the results as a vector. I think doing a for loop would work, but I cant work out how to save the values generated from the function as a vector. Any ideas? Thanks. -- View this message in context:

[R] Instrumental variable quantile estimation of spatial autoregressive models

2011-05-04 Thread Marie-Line Glaesener
Dear all, I would like to implement a spatial quantile regression using instrumental variable estimation (according to Su and Yang (2007), Instrumental variable quantile estimation of spatial autoregressive models, SMU economics statistis working paper series, 2007, 05-2007, p.35 ). I am

[R] Distance measure in heatmaps

2011-05-04 Thread Muhidini
Hi, First- I am relatively new to R and would therefore greatly appreciate any kind of help from you ! I am currently trying to use R to make a heatmap using the Pheatmap package. My question: I want to include a different distance measure than the ones given by the implemented hclust

[R] Uniform Gaussian Kernel

2011-05-04 Thread blutack
I have a vector with lots of different numbers. I need to make a graph showing the Uniform Distribution of the figures. I have created a graph showing all the different values, but now want individual Gaussian Kernel round each point. This is what I have but each time it comes up with an error as

Re: [R] Superscript number before letter

2011-05-04 Thread Janhal
Yes thats it :-) Thank you very much! Janine -- View this message in context: http://r.789695.n4.nabble.com/Superscript-number-before-letter-tp3495577p3495812.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

Re: [R] Storing data from a test as a vector or matrix

2011-05-04 Thread wwreith
SA gives the output: Response IPS1 : Df Sum Sq Mean Sq F value Pr(F) as.factor(WSD) 3 3.3136 1.10455 23.047 5.19e-12 *** Residuals 129 6.1823 0.04793 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 . . . There are 11 more just

Re: [R] Storing data from a test as a vector or matrix

2011-05-04 Thread wwreith
I figured out that attributes is the command that I was trying to find. It allowed me to find out that I was needing to use stats not Df or Pillai etc. Following command worked. S1-as.vector(S$stats[1,]) However when I try the same thing with summary.aov it is not working. SA-summary.aov(M)

  1   2   >