Re: [R] R-patched-2.10.1-50.1.x86_64 installlation error

2010-02-15 Thread Detlef Steuer
Hi! libreadline6 exists only inside Factory, not in 11.1. Are you sure you clicked on the correct link? Otherwise I have no quick idea how a readline6 could slip into the build process. Report back, please, if the problem persists. Detlef On Tue, 16 Feb 2010 04:11:48 +0100 wrote: > here is t

Re: [R] Total and heading of portfoilo table

2010-02-15 Thread Madhavi Bhave
Hi!   I am not expert in R, but perhaps you can try the following -   X = as.numeric(read.csv('quantity.csv')) Y = read.csv('equity_price.csv') Y = Y[, -1]   Z = X*Y   port_val = NULL   for(i in 1 : nrow(Z)) {   port_val[i] = sum(Z[i,])   }   write.csv(data.frame(Z, port_val = port_val), 'PORTFOLIO

[R] Total and heading of portfoilo table

2010-02-15 Thread Sarah Sanchez
Dear R helpers, I have two input files as 'quantity.csv' and 'equity_price.csv' as (for example) given below. 'quantity.csv' GOOG YHOO 1000 100 'equity_price.csv' sr_no   GOOG_price   YHOO_price 1    15.22 536.40 2    15.07 53

[R] survival - ratio likelihood for ridge coxph()

2010-02-15 Thread Damjan Krstajic
It seems to me that R returns the unpenalized log-likelihood for the ratio likelihood test when ridge regression Cox proportional model is implemented. Is this as expected? In the example below, if I am not mistaken, fit$loglik[2] is unpenalized log-likelihood for the final estimates of coeffi

Re: [R] error message when downloading packages using the OS X shell

2010-02-15 Thread Ista Zahn
Hi Jessica, As far as I know that's the expected behavior (at least it's not OS X-specific: I get the same on OpenSuse 11.2). What happens if you do download.packages("ape", destdir = "~/Desktop/") ? Also, out of curiosity, why are you using download.packages() instead of install.packages()? Be

[R] HELP on Non-Linera Mixed-Effect model

2010-02-15 Thread vaibhav dua
Hi, I'm trying to fit nonlinear mixed effects model using nlme function but getting an error message. Here is what I have: fitted_model = nlme(scores~spline(b1,b2,b3,kt,time), fixed = list(b1~1, b2~1, b3~1, kt~1), random = b1+b2+b3~1, groups= ~id, da

Re: [R] error message when downloading packages using the OS X shell

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 10:51 PM, Jessica Joganic wrote: I apologize for not including my entire script. What I typed into the shell was: download.packages(ape) to which R responded with a Tcl/Tk interface allowing me to set my CRAN. After I did so it proceeded to spit out the following error

[R] Outliers detection using wavelets

2010-02-15 Thread assaedi76 assaedi76
Hi R users   Thanks in advance   I have time series data sets. Is there any package in R that helps detect the outliers  using wavelet methods.   Thanks   email: ( assaed...@yahoo.com ) [[alternative HTML version deleted]] __ R-help@r-p

Re: [R] error message when downloading packages using the OS X shell

2010-02-15 Thread Jessica Joganic
I apologize for not including my entire script. What I typed into the shell was: *download.packages(ape)* to which R responded with a Tcl/Tk interface allowing me to set my CRAN. After I did so it proceeded to spit out the following error: *Loading Tcl/Tk interface ... done* Error in **file.info

Re: [R] difftimes; histogram; memory problems

2010-02-15 Thread Gabor Grothendieck
Just one further point. If you do run out of memory using #2 then try this which is the same as #2 but adds a dbname argument to force the computation to be done from disk rather than memory. sqldf("select d1.x - d2.x, count(*) from d1, d2 group by d1.x - d2.x", dbname = tempfile()) On Mon, Feb 1

Re: [R] difftimes; histogram; memory problems

2010-02-15 Thread Moshe Olshansky
Hi Jonathan, If minDate = min(Condition1) - max(Condition2) and maxDate = max(Condition1) - min(Condition2) then all your differences would be between minDay and maxDay, and hopefully this is not a very big range (unless you are going many thousands years into the past or the future). So basica

Re: [R] difftimes; histogram; memory problems

2010-02-15 Thread Gabor Grothendieck
Here are two approaches to try: > # test data > d1 <- data.frame(x = Sys.Date() + 1:3) > d2 <- data.frame(x = Sys.Date() - 1:3) > # 1. you might not have enough memory for this but its short > table(outer(1:3, -(1:3), "-")) 2 3 4 5 6 1 2 3 2 1 > # 2. this one performs all the operations outsid

[R] R-patched-2.10.1-50.1.x86_64 installlation error

2010-02-15 Thread mauede
here is the message I get upon trying to upgrade R from the provided RPM package for SuSE 11.1 YaST2 conflicts list - generated 2010-02-15 14:00:50 nothing provides libreadline.so.6()(64bit) needed by R-patched-2.10.1-50.1.x86_64 [ ] do not install R-patched-2.10.1-50.1.x86_64

Re: [R] error message when downloading packages using the OS X shell

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 8:58 PM, Jessica Joganic wrote: Hi Fellow R Users, I recently upgraded to Mac OS X 10.5 (Snow Leopard) and had some issues downloading and running R 2.10.1. I fixed the tcl/tk problem I was originally having but it was replaced with another. I run R out of the shell (t

Re: [R] projections

2010-02-15 Thread Ray Brownrigg
If you use map() with the projection= option, then see ?mapproject, which allows you to project arbitrary (long/lat) points with the previously defined projection. Then you use points() to overlay them onto the map. HTH, Ray Brownrigg On Tue, 16 Feb 2010, Shaun Langley wrote: > I have lat/lo

Re: [R] Hypercube in R

2010-02-15 Thread Hartmut Oldenbürger
Hi Andrej, Uwe, Bill isn't the central keyword 'tesseract', on which and it's extensions we have extensive literature? e.g., starting from http://en.wikipedia.org/wiki/Tesseract I remember an article by Martin Gardner on the construction principles, too.

[R] projections

2010-02-15 Thread Shaun Langley
I have lat/lon coordinates in Michigan GeoRef that I need to re-project so they can be shown over a Michigan county map generated from the Maps library. I'm very confused as to how to do this. Can anyone help? Thanks, Shaun [[alternative HTML version deleted]]

Re: [R] difftimes; histogram; memory problems

2010-02-15 Thread Jonathan
Let me fix a couple of typos in that email: Hi All: Let's say I have two dataframes (Condition1 and Condition2); each being on the order of 12,000 and 16,000 rows; 1 column. The entries contain dates. I'd like to calculate, for each possible pair of dates (that is: Condition1[1:12,000] and Cond

[R] error message when downloading packages using the OS X shell

2010-02-15 Thread Jessica Joganic
Hi Fellow R Users, I recently upgraded to Mac OS X 10.5 (Snow Leopard) and had some issues downloading and running R 2.10.1. I fixed the tcl/tk problem I was originally having but it was replaced with another. I run R out of the shell (terminal) and when I ask it to download.packages() it gives me

[R] tm package

2010-02-15 Thread David Neu
Hi, I'm using version 0.5.1 of tm package with R 2.10.1. It looks to me as if after the following reuters21578 <- Corpus(DirSource(corpusDir), readerControl = list(reader = readReut21578XMLasPlain)) reuters21578 <- tm_map(reuters21578, stripWhitespace) reuters21578 <- tm_map(reuters

[R] Help with lsei() from package limSolve()

2010-02-15 Thread Liam Blanckenberg
Hi all, The problem I am trying to solve is the following: I have a 'trace' of half-hourly data across an entire year - i.e. I have 17,500 observations (one for every half-hour of every day of a year). The data follows a daily, weekly and seasonal pattern that is important to preserve. What I want

[R] difftimes; histogram; memory problems

2010-02-15 Thread Jonathan
Hi All: Let's say I have two dataframes (Condition1 and Condition2); each being on the order of 12,000 and 16,000 rows; 1 column. The entries contain dates. I'd like to calculate, for each possible pair of dates (that is: Condition1[1:10,000] and Condition2[1:10,000], the number of days differen

Re: [R] Predicted points in splines

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 6:24 PM, Oddie wrote: I have a data set (date, concentration) were i applied the spline interpolation. I would like to use, though, the Spline () command R is a functional language. We don't "do" commands. You give us something, we return something. in order to get

[R] Predicted points in splines

2010-02-15 Thread Oddie
I have a data set (date, concentration) were i applied the spline interpolation. I would like to use, though, the Spline () command in order to get the predicted values that correspond not only to points used but also to the points inbetween. Is there a way to set the time step for which i will ge

Re: [R] What symbols are loaded when require()ing?

2010-02-15 Thread Greg Snow
The ls() or objects() functions will be a place to start. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of Gui

[R] a fast way to get AIC, BIC and logLik?

2010-02-15 Thread Esmail
Hello all, Disclaimer, I know very little about statistics :-| I am running a series of regressions for friends, the goal being to determine the AIC, BIC and logLik values from the result/summary of each regression. We are doing this as part of an evolutionary algorithm, so it is important that

Re: [R] Error analysis for circular data

2010-02-15 Thread William Revelle
See the packages: CircSpatial CircStats circular At 10:58 PM +0100 2/15/10, Francisco Javier Santos Alamillos wrote: Dear R users, I would like to know if it is possible to calculate the Mean Error (ME), the Root Mean-squared error (RMSE) and absolute error (MAE) for two temporal series of

Re: [R] How to check the two different nulls?

2010-02-15 Thread Erik Iverson
blue sky wrote: x=list(a=1,b=NULL) is.null(x$b) is.null(x$c) Both the above two commands give me TRUE, but in the first one, b is NULL, in the second one, c doesn't exist. Are there functions that can help me distinguish the two different nulls? __ R-

[R] How to check the two different nulls?

2010-02-15 Thread blue sky
x=list(a=1,b=NULL) is.null(x$b) is.null(x$c) Both the above two commands give me TRUE, but in the first one, b is NULL, in the second one, c doesn't exist. Are there functions that can help me distinguish the two different nulls? __ R-help@r-project.org

Re: [R] plot more series or more colors

2010-02-15 Thread Jim Lemon
On 02/16/2010 09:50 AM, joonR wrote: Hey, I need either to plot more than one series of data on the same plot, or to color the points of a single series with different colors according to some parameter, how can I do that? (data are 2 coordinates point that I want to plot in 2D) Hi joonR, Loo

Re: [R] plot more series or more colors

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 5:50 PM, joonR wrote: Hey, I need either to plot more than one series of data on the same plot, or to color the points of a single series with different colors according to some parameter, how can I do that? (data are 2 coordinates point that I want to plot in 2D) ?pl

[R] plot more series or more colors

2010-02-15 Thread joonR
Hey, I need either to plot more than one series of data on the same plot, or to color the points of a single series with different colors according to some parameter, how can I do that? (data are 2 coordinates point that I want to plot in 2D) THANKS! -- View this message in context: http://n4.

Re: [R] Problems with boxplot in ggplot2:qplot

2010-02-15 Thread Dimitri Shvorob
library(sqldf) library(ggplot2) t = data.frame(t = seq.Date(as.Date("2009-01-01"), to = as.Date("2009-12-01"), by = "month")) x = data.frame(x = rnorm(5)) df = sqldf("select * from t, x") qplot(factor(df$t), df$x, geom = "boxplot") + theme_bw() qplot(factor(df$t), df$x, geom = "boxplot") + the

Re: [R] help with EXPASY HTML form submission in RCurl package

2010-02-15 Thread Sunando Roy
Hi Duncan, I finally got this to work, i.e. to extract my file two pages down. Thanks a lot for your help. I think there is some problem with the HTML source code of "http://www.expasy.ch/tools/protscale.html";. When I use the getHtmlFormDescription function the output that I get is not what I wou

Re: [R] Confidence intervals nls

2010-02-15 Thread Walmes Zeviani
The problem was because attr(predict(model, list(DOY=days))) didn't provide any gradient, it returned NULL. To correct this we need to specify the gradient by using deriv3() function. Then, when you use predict(model, list(DOY=days)) will be provided the gradient at this new experimental points.

Re: [R] Using text put into a dialog box

2010-02-15 Thread Greg Snow
You may want to use tkwait.window(tt) If the below is in a function, then tkwait.window(tt) will wait until the window goes away (the user clicking on the button) before continuing with the code, then you can return the result after that. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Stati

[R] Variance functions in nlme: varFixed with nlme() function.

2010-02-15 Thread Ines Naya
Dear R-help suscribers, I am doing a meta-analysis of sea urchin growth data in R. I am fitting a non-linear growth function using nlme(). Most of my observations are means, and I want to give them weights according to the number of individuals that were used to obtain those means (as I do not h

Re: [R] Non-metric multi-dimensional scaling with "NA"s / "0"s

2010-02-15 Thread Sarah Goslee
What do you mean by "ignored"? Do you want to drop those sites from the ordination? If so, you need to subset your data beforehand. Or if you mean that you have presence data, and not really absences (common in ecology), and want to ignore joint absences, then you need to select a dissimilarity me

[R] Error analysis for circular data

2010-02-15 Thread Francisco Javier Santos Alamillos
Dear R users, I would like to know if it is possible to calculate the Mean Error (ME), the Root Mean-squared error (RMSE) and absolute error (MAE) for two temporal series of directional data. Where Can I get documentation about it? Thanks in advance [[alternative HTML version deleted]

Re: [R] Shapiro-Wilk for levels of factor

2010-02-15 Thread Greg Snow
Achim showed you how, but you might want to consider why. If you are trying to learn more about your data, then plots or other strategies may work better than the test. If you are testing for normality in order to meet the assumptions of a test, then the test may not be accomplishing what you t

Re: [R] Extract values from a predict() result... how?

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 4:35 PM, Jay wrote: Hello, silly question I suppose, but somehow I can't manage to extract the probabilities from a glm.predict() result: str(res) Named num [1:9] 0.00814 0.01877 0.025 0.02941 0.03563 ... - attr(*, "names")= chr [1:9] "1" "2" "3" "4" ... I got from: #

[R] Non-metric multi-dimensional scaling with "NA"s / "0"s

2010-02-15 Thread kellys17
Hi, I am tryng to use Non-metric multi-dimensional (MDS) to illustrate the difference in mixed-species flock composition between two species of birds. It is part of a thesis I am currently working on. My data consists of 18 columns and 62 rows, with numbers ranging from 0 to 7. I have car

[R] Extract values from a predict() result... how?

2010-02-15 Thread Jay
Hello, silly question I suppose, but somehow I can't manage to extract the probabilities from a glm.predict() result: > str(res) Named num [1:9] 0.00814 0.01877 0.025 0.02941 0.03563 ... - attr(*, "names")= chr [1:9] "1" "2" "3" "4" ... I got from: # A Gamma example, from McCullagh & Nelder (

Re: [R] Problems with boxplot in ggplot2:qplot

2010-02-15 Thread hadley wickham
Without a reproducible example, it's impossible to give you any more suggestions. Hadley On Mon, Feb 15, 2010 at 2:16 PM, Dimitri Shvorob wrote: > > Trying > > + scale_x_date(format = "%b") > > produces a peculiar result: Apr and Dec facets are labeled "Jan", remaining > labels are blank. > -- >

Re: [R] (no subject)

2010-02-15 Thread Liviu Andronic
Hello You are more likely to get help if you provide a subject more informative than "no subject". On 2/15/10, hussain abu-saaq wrote: > is there a good notes for loops(if else and while) in r. > Check this [1][3], but also the R manuals [2]. [1] http://www.r-project.org/doc/Rnews/Rnews_2008-

Re: [R] Problems with boxplot in ggplot2:qplot

2010-02-15 Thread Dimitri Shvorob
Trying + scale_x_date(format = "%b") produces a peculiar result: Apr and Dec facets are labeled "Jan", remaining labels are blank. -- View this message in context: http://n4.nabble.com/Problems-with-boxplot-in-ggplot2-qplot-tp1555338p1556573.html Sent from the R help mailing list archive at N

Re: [R] Problems with boxplot in ggplot2:qplot

2010-02-15 Thread Dimitri Shvorob
Thank you, Hadley. I try jpeg(file, width = 800, height = 600, quality = 100) qplot(factor(closed$close.month), closed$closing.balance, geom = "boxplot", main = "Monthly distributions of closing balances", xlab = "Month", ylab = "Balance, USD") + theme_bw() + scale_x_date(major = "months",

Re: [R] Confidence intervals nls

2010-02-15 Thread FishR
Dear Peter Thank you that is an error in the code but unfortunately the problem is still apparent. I think it is related to calculating the gradient Tom -- View this message in context: http://n4.nabble.com/Confidence-intervals-nls-tp1556487p1556585.html Sent from the R help mailing list arch

[R] (no subject)

2010-02-15 Thread hussain abu-saaq
is there a good notes for loops(if else and while) in r. I have a code written in Matlab and I tried to rewrite it in R but it does not give me the write result . for example in matlab the result is (1 2 3 4) In r it is give my only 2 One more question if I have variable,lets say x and I use

Re: [R] Confidence intervals nls

2010-02-15 Thread Peter Ehlers
On 2010-02-15 12:01, FishR wrote: Dear All I am quite new to R and would appreciate some help fitting 95% confidence intervals to a nls function. I have the data DOY CET 90 5.9 91 8 92 8.4 93 7.7 95 6.6 96 6.8 97 7.1 98 9.7 99 12.3 100 12.

Re: [R] tree-drawing in R ?

2010-02-15 Thread William Revelle
try Rgraphviz diagram At 6:07 PM + 2/15/10, Oliver Kullmann wrote: Hello, I wonder whether binary (rooted) trees with simple labels (say, integers) can be drawn by some R-package? Couldn't find one. Just to make sure (since "trees" can be many things): I mean those computer-science creatu

Re: [R] creating functions question

2010-02-15 Thread AC Del Re
Perfect! Thank you, Dimitris. AC On Mon, Feb 15, 2010 at 1:45 PM, Dimitris Rizopoulos wrote: > try this: > > MRfit <- function (...) { >    models <- list(...) >>    do.call(anova, models) >> } >> >> >> I hope it helps. >> >> Best, >> Dimitris >> >> >> AC Del Re wrote: >>> >>> Hi All, >>>

[R] Adjusted means and generalized chain block designs

2010-02-15 Thread Giles
Dear Colleagues, John Mandel ( Chain block designs with two-way elimination of heterogeneity. Biometrics 10, 251-272 ,1954). extended the class of chain block designs (Youden & Conner (1953) to elimination of both row and column (blocks) effects. These experimental designs can be useful in engin

Re: [R] creating functions question

2010-02-15 Thread Dimitris Rizopoulos
try this: MRfit <- function (...) { models <- list(...) do.call(anova, models) } I hope it helps. Best, Dimitris AC Del Re wrote: Hi All, I am interested in creating a function that will take x number of lm objects and automate the comparison of each model (using anova). Here is a

Re: [R] using step() with package geepack

2010-02-15 Thread Steve Taylor
Is anyone else using the 'geepack' package? >>> From: "Steve Taylor" To: Date: 11/Feb/2010 11:28 a.m. Subject: [R] using step() with package geepack I'm using the package geepack to fit GEE models. Does anyone know of methods for add1 and drop1 for a 'geeglm' model object, or perhaps a method

[R] creating functions question

2010-02-15 Thread AC Del Re
Hi All, I am interested in creating a function that will take x number of lm objects and automate the comparison of each model (using anova). Here is a simple example (the actual function will involve more than what Im presenting but is irrelevant for the example): # sample data: id<-rep(1:20)

Re: [R] Selective load of .First() function just for Rgui.exe

2010-02-15 Thread Steve Taylor
Here's another way to have different behaviour only for interactive sessions: if (interactive()) { winMenuAdd(menuname, NULL, NULL) # etc. } >>> From: Duncan Murdoch To: CC: Date: 13/Feb/2010 6:15 a.m. Subject: Re: [R] Selective load of .First() function just for Rgui.exe On 12/02/2010 11:

Re: [R] Problems with boxplot in ggplot2:qplot

2010-02-15 Thread hadley wickham
Hi Dimitri, Have you looked at the examples for scale_x_date - http://had.co.nz/ggplot2/scale_date.html? They show you how to both set the limits and control the labels. Hadley On Sun, Feb 14, 2010 at 1:34 PM, Dimitri Shvorob wrote: > > ... Unfortunately, a problem remains: I cannot label x ti

[R] Confidence intervals nls

2010-02-15 Thread FishR
Dear All I am quite new to R and would appreciate some help fitting 95% confidence intervals to a nls function. I have the data DOY CET 90 5.9 91 8 92 8.4 93 7.7 95 6.6 96 6.8 97 7.1 98 9.7 99 12.3 100 12.8 102 11 103 9.3 104 9.8

Re: [R] Sampling from Bivariate Uniform Distribution

2010-02-15 Thread Haneef_An
When I wrap those values in to a matrix will it be still independent ? ( non zero correlation). Can I do this for any multivariate distribution which has the univariate form? Thank you for the response. Haneef -- View this message in context: http://n4.nabble.com/Sampling-from-Bivariate-Unif

Re: [R] Non-monotonic spline using splinefun(method = "monoH.FC")

2010-02-15 Thread Tim Heaton
Hi, Thanks for the reply but I think you are confusing monotonic and strictly increasing/decreasing. I also just used the y-value of the last knot as a simple example as that is not the bit where it goes wrong. It will still produce a non-monotonic spline if you use for example x <- 1:8 y <- c(-

Re: [R] Adressing multiple cores (CPUs)

2010-02-15 Thread Cedrick W. Johnson
Since you're on Windows, try the 'snow' package... the other suggestions are highly suggested and recommended as well hth c thomas.schwan...@mvv.de wrote: Dear all, I'm sitting here just in front of my new p...@work and wonder about the following question: * How can I adress multipl

Re: [R] Which method is called when print(a_list)?

2010-02-15 Thread Erik Iverson
blue sky wrote: I don't find print.list. Could somebody let me know which method is called when I run command print(a_list), where a_list is a list? Is 'print.default' used for printing a list? Yes. You can always debug functions to investigate what's going on, too. See ?debug. __

[R] Which method is called when print(a_list)?

2010-02-15 Thread blue sky
I don't find print.list. Could somebody let me know which method is called when I run command print(a_list), where a_list is a list? Is 'print.default' used for printing a list? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] rlaplace using rmutil - HELP

2010-02-15 Thread Peter Ehlers
On 2010-02-15 9:38, vaibhav dua wrote: Hi Peter, Thanks for quick reply. Like I said that I'm aware that it needs a vector but I want to use diagonal matrix. As you have noted, It could be possible that rlaplace reading dispDiag as a vector and not matrix like c(3,0,0,0,.20,0,0,0,.1) so I tho

Re: [R] Non-monotonic spline using splinefun(method = "monoH.FC")

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 1:06 PM, Tim Heaton wrote: Hi, Thanks for the reply but I think you are confusing monotonic and strictly increasing/decreasing. I also just used the y-value of the last knot as a simple example as that is not the bit where it goes wrong. It will still produce a non-mon

[R] tree-drawing in R ?

2010-02-15 Thread Oliver Kullmann
Hello, I wonder whether binary (rooted) trees with simple labels (say, integers) can be drawn by some R-package? Couldn't find one. Just to make sure (since "trees" can be many things): I mean those computer-science creatures, with roots and subtrees. Thanks for your consideration in any case Ol

Re: [R] Which method is called in command like "class(x)='something'"?

2010-02-15 Thread Barry Rowlingson
On Mon, Feb 15, 2010 at 5:07 PM, blue sky wrote: >> x=3 >> `class<-`(x,'something')#this command prints > [1] 3 > attr(,"class") > [1] "something" >> x=3 >> class(x)='something'#this command doesn't print anything > > The first of the above two commands print the content of 'x' but the > second do

Re: [R] Non-monotonic spline using splinefun(method = "monoH.FC")

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 10:59 AM, Tim Heaton wrote: Hi, In my version of R, the stats package splinefun code for fitting a Fritsch and Carlson monotonic spline does not appear to guarantee a monotonic result. If two adjoining sections both have over/undershoot the way the resulting adjustment of a

[R] R still accepts the undocumented option --with-iconv=no

2010-02-15 Thread Dr. David Kirkby
This is not so much a user help question, as a bug report (I think!) Previous versions of R allowed the option --with-iconv=no to used to disable the use of iconv. That option is no longer documented. However, using that option on the latest does not generate an error message, but is silently i

[R] Non-monotonic spline using splinefun(method = "monoH.FC")

2010-02-15 Thread Tim Heaton
Hi, In my version of R, the stats package splinefun code for fitting a Fritsch and Carlson monotonic spline does not appear to guarantee a monotonic result. If two adjoining sections both have over/undershoot the way the resulting adjustment of alpha and beta is performed can give modified values

[R] Which method is called in command like "class(x)='something'"?

2010-02-15 Thread blue sky
> x=3 > `class<-`(x,'something')#this command prints [1] 3 attr(,"class") [1] "something" > x=3 > class(x)='something'#this command doesn't print anything The first of the above two commands print the content of 'x' but the second doesn't, although both of them set the argument 'x'. I'm wondering

Re: [R] Mailing List

2010-02-15 Thread stephen sefick
Go through the procedure outlined where you signed up. I can't take you off of the list. On Mon, Feb 15, 2010 at 10:39 AM, Audra Taylor wrote: > I would like to be removed from the R mailing list > > > >        [[alternative HTML version deleted]] > >

Re: [R] density estimates for fixed points

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 10:02 AM, geir wrote: Sorry, i did indeed mean function - not package. kde2d lets you determine the grid - you choose "n" and "lims", which gives you and n x n matrix for which the density is estimated (Seems like the output is meant for graphical purposes). But i have

[R] Mailing List

2010-02-15 Thread Audra Taylor
I would like to be removed from the R mailing list [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-g

Re: [R] lattice/ylim: how to fix ylim[1], but have ylim[2] dynamically calculated?

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 2:29 AM, Johannes Graumann wrote: Rolf Turner wrote: On 15/02/2010, at 9:40 AM, Johannes Graumann wrote (In response to some advice from David Winsemius): I am quite certain that this is the most elaborately worded version of "RTFM" I have ever come across. I no

Re: [R] rlaplace using rmutil - HELP

2010-02-15 Thread vaibhav dua
Hi Peter, Thanks for quick reply. Like I said that I'm aware that it needs a vector but I want to use diagonal matrix. As you have noted, It could be possible that rlaplace reading dispDiag as a vector and not matrix like c(3,0,0,0,.20,0,0,0,.1) so I thought to change the code little bit: d

Re: [R] argh .. if/else .. why?

2010-02-15 Thread Gabor Grothendieck
On Mon, Feb 15, 2010 at 11:24 AM, David Winsemius wrote: > > On Feb 15, 2010, at 11:01 AM, hadley wickham wrote: > >>> I, personally, utilize the ifelse(test,statement,statement) function when >>> possible over the methodology outlined. >> >> if + else and ifelse perform quite different tasks, and

Re: [R] get current window (device) size

2010-02-15 Thread Edmond.Ng
Hi David, Thank you for your advice. I have just figured why neither the suggestion of 'par("din")' or "dev.size()" works for me. Having specified 'rescale="fixed"' in a preceding windows() command seems to prevent R from reporting the current window size. See here: windows(width=7.5,height=6,r

Re: [R] argh .. if/else .. why?

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 11:01 AM, hadley wickham wrote: I, personally, utilize the ifelse(test,statement,statement) function when possible over the methodology outlined. if + else and ifelse perform quite different tasks, and in general can not (and should not) be exchanged. In particular, note

Re: [R] argh .. if/else .. why?

2010-02-15 Thread hadley wickham
> I, personally, utilize the ifelse(test,statement,statement) function when > possible over the methodology outlined. if + else and ifelse perform quite different tasks, and in general can not (and should not) be exchanged. In particular, note that for ifelse, "the class attribute of the result

Re: [R] Adressing multiple cores (CPUs)

2010-02-15 Thread Sharpie
Benilton Carvalho-2 wrote: > > http://cran.r-project.org/web/views/HighPerformanceComputing.html > Also, of the links at the bottom of that page, I found the journal article: "State of the Art in Parallel Computing with R" To be a particularly helpful and easy to read overview of the avail

Re: [R] rlaplace using rmutil - HELP

2010-02-15 Thread Peter Ehlers
On 2010-02-15 8:26, vaibhav dua wrote: Hi, I'm using rlaplace distribution to draw non normal random deviates. I'm using rmutil package for this. I know rlaplace needs location and dispersion parameters but I want to pass diag matrix as dispersion parameters and I'm not getting the correct va

Re: [R] Which is the proper mail list?

2010-02-15 Thread Sharpie
Roy Mendelssohn wrote: > > Hi All: > > Even after reading the web pages, it is unclear to me exactly how it is > decided what goes into the base distribution, and more to the point, which > list I should post to if I want to suggest that a functionality now done > by a contributed package be a

[R] Which is the proper mail list?

2010-02-15 Thread Roy Mendelssohn
Hi All: Even after reading the web pages, it is unclear to me exactly how it is decided what goes into the base distribution, and more to the point, which list I should post to if I want to suggest that a functionality now done by a contributed package be added to the base packages. Thanks fo

[R] rlaplace using rmutil - HELP

2010-02-15 Thread vaibhav dua
Hi, I'm using rlaplace distribution to draw non normal random deviates. I'm using rmutil package for this. I know rlaplace needs location and dispersion parameters but I want to pass diag matrix as dispersion parameters and I'm not getting the correct values. Here is my code snippet: disp

Re: [R] Is this a bug? (layout and par)

2010-02-15 Thread Peter Ehlers
On 2010-02-15 7:06, Martin Ivanov wrote: Dear R users, I think I have spotted a bug in R, but as I am not sure, I first post it to you. Here is a minimal working example: default<- par(no.readonly=TRUE); par(font=1, adj=0.5, cex=1, cex.lab=1, cex.axis=0.5, font.axis=2, lend=2, family="Times"

Re: [R] executable R script under xp (to avoid migration toward Matlab or C++)

2010-02-15 Thread PtitBleu
Hello everyone, Barry : Thanks for the link to Rscript. I will check this possibility (the os was not given in the text, but in subject :-)) Romain : I think my skills at computing are too low for your solution (I only know a bit of R and of Quick Basic ... nothing with C, C++, ...) and I also d

Re: [R] get current window (device) size

2010-02-15 Thread David Winsemius
On Feb 15, 2010, at 9:49 AM, wrote: Dear Prof Ripley, Thank you very much for your advice. In fact, I have already tried 'par("din")'. It does not seems to pick up the changes made to the window. For example, I issued the par command after the plot: par("din") [1] 7.5 7.0 [changed window si

Re: [R] density estimates for fixed points

2010-02-15 Thread geir
Sorry, i did indeed mean function - not package. kde2d lets you determine the grid - you choose "n" and "lims", which gives you and n x n matrix for which the density is estimated (Seems like the output is meant for graphical purposes). But i have a matrix consisting of points (x,y) where i want

Re: [R] get current window (device) size

2010-02-15 Thread Edmond.Ng
Dear Prof Ripley, Thank you very much for your advice. In fact, I have already tried 'par("din")'. It does not seems to pick up the changes made to the window. For example, I issued the par command after the plot: > par("din") [1] 7.5 7.0 [changed window size] > par("din") [1] 7.5 7.0 The res

Re: [R] argh .. if/else .. why?

2010-02-15 Thread Peter Dalgaard
Esmail wrote: > Ah .. so the fact that this is inside a block (in this case a function) > makes a difference .. I see .. well this is certainly different from other > languages I have worked with, but the explanation is clear. Yes. I think most people involved have realized that it was a design b

Re: [R] Adressing multiple cores (CPUs)

2010-02-15 Thread Benilton Carvalho
http://cran.r-project.org/web/views/HighPerformanceComputing.html 2010/2/15 : > Dear all, > > I'm sitting here just in front of my new p...@work and wonder about the > following question: >        * How can I adress multiple CPUs (cores) out of R to speed up the > simulations I run? >        *

[R] (no subject)

2010-02-15 Thread mauede
here is the message I get upon trying to upgrade R from the provided RPM package for SuSE 11.1 (http://cran.r-project.org): YaST2 conflicts list - generated 2010-02-15 14:00:50 nothing provides libreadline.so.6()(64bit) needed by R-patched-2.10.1-50.1.x86_64 [ ] do not install R-

[R] Adressing multiple cores (CPUs)

2010-02-15 Thread thomas.schwander
Dear all, I'm sitting here just in front of my new p...@work and wonder about the following question: * How can I adress multiple CPUs (cores) out of R to speed up the simulations I run? * What are the prerequisites to do so? Maybe anyone could give me a hint where to start read

Re: [R] argh .. if/else .. why?

2010-02-15 Thread Esmail
On 15-Feb-10 09:27, lee.ad...@luminant.com wrote: I, personally, utilize the ifelse(test,statement,statement) function when possible over the methodology outlined. Haven't used that construct before, I'll check it out. Thanks. And mystery solved. Esmail

Re: [R] argh .. if/else .. why?

2010-02-15 Thread Esmail
On 15-Feb-10 09:13, jim holtman wrote: This is a common question. The problem is basically when input commands at the command line, when you hit the carriage return it is assuming there is a complete command. Within a function, or a block like the 'for' loop, the statement is not complete until

Re: [R] get current window (device) size

2010-02-15 Thread Prof Brian Ripley
On Mon, 15 Feb 2010, edmond...@lshtm.ac.uk wrote: Dear R-listers, I am writing a function that will generate a simple plot for the user . The size of the window (that contains the plot) may not be what the user wants. Supposing the user adjusts the window to a size of her/his choice, I wou

Re: [R] Problem with R on USB-drive

2010-02-15 Thread Peter Dalgaard
Lauri Nikkinen wrote: > The answer is: > >> .libPaths() > [1] "E:/R/R-2.10.1/library" > > stats package is found under > > .../library/ Sorry about the delay Lauri, I don't have the time to skim r-help very well these days, so your note sort of whizzed by. However, I'm not sure I get the p

[R] installation of package 'e1071' had non-zero exit status

2010-02-15 Thread mauede
I tried twice to install package "e1071"as it provides the Hamming distance. The installation failed twice. I am running R version 2.9.0 (2009-04-17). Any suggestion is welcome. Thank you, Maura tutti i telefonini TIM! [[alternative HTML version deleted]]

  1   2   >