[R] Re ading Image Files

2009-07-23 Thread cotixan
Hello, I'm rather new to R and I want to do some image analysis. Is there a way to read jpeg files into a matrix like matlab's imread? -- View this message in context: http://www.nabble.com/Reading-Image-Files-tp24619129p24619129.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] IDE's/Editors for R-- Was: Multi-line comments?

2009-07-23 Thread Romain Francois
On 07/22/2009 06:49 PM, Liviu Andronic wrote: Hello, On Wed, Jul 22, 2009 at 6:25 PM, Bert Guntergunter.ber...@gene.com wrote: You can find a list of IDE's/R code editors for R here: http://www.sciviews.org/_rgui/projects/Editors.html However, this is somewhat dated, and you may find others

[R] Antw: How to extract the upper xlim and ylim of my plot?

2009-07-23 Thread Frank Bloos
You may want to use the function corner.label from the plotrix-package. Frank Mark Na mtb...@gmail.com 21.07.2009 23:03 Dear R-helpers, I wish to place some text in a plot, at approx 10% of my upper xlim and approx 90% of my upper ylim, i.e. plot(log(all$SR,10)~log(all$AREA,10))

Re: [R] Re placing null values (#NULL!)

2009-07-23 Thread Michael Knudsen
On Wed, Jul 22, 2009 at 10:18 PM, Josh Rollj_r...@hotmail.com wrote: In `[-.factor`(`*tmp*`, Props_$pct_vacant == #NULL!, value = 0) :   invalid factor level, NAs generated Thats what made me wonder if the #NULL! value was being treated differently than a typical NULL value.  Thoughts? I

[R] [R-pkgs] plotrix v2.6-4

2009-07-23 Thread Jim Lemon
Hi all, After a lot of rewriting, I am happy to say that plotrix_2.6-4 is on CRAN. The reason I am announcing this is that the newish hierobarp function has had a complete rewrite. This was the result of my learning that the old version was horrendously slow on large data sets. Let me know if

[R] simple question about exporting data...

2009-07-23 Thread Rbeginner
I just have a simple question. I'm using the predict function, and I get about 1000 values. Now, how would I make it into a data frame or a list so that I can easily use Excel to graph the data? Right now, it's in horizontal lines, and it would be a pain to copy and paste and make it into a

[R] How to install packages without Tcl/Tk interface

2009-07-23 Thread Matej Kovacic
Hi, how to install packages without Tcl/Tk interface? For instance, if I run install.packages('R.oo'), I need to select Cran repository via graphic interface. How to avoid this? bye, Matej __ R-help@r-project.org mailing list

[R] Network from package functions

2009-07-23 Thread Andrej-Nikolai Spiess
Dear R-helpers, does anyone know of some package/function that can build a network from the functions that are implemented in a package, i.e. visualize the cross-references from one function to another in the same or some dependent package? An example would be a function like 'nls' on top of the

Re: [R] Memory errors when using QCA package

2009-07-23 Thread Adrian Dusa
Allan Engelhardt allane at cybaea.com writes: It is a little stupid, but the length of a vector is limited to 2^31-1 entries on any platform. A matrix is stored as a vector, so the product of all dimensions is also limited to 2^31-1. Allan. Matthew Gwynne wrote: Hi, I have been

[R] SystemFit

2009-07-23 Thread Ferdogan
Hi, I have two products which are substitudes. I try to fix a system as below to mydata. Demand1 = A1 -B1*Price1 + C1*Price2 Demand2 = A2 +B2*Price1 - C2*Price2 I would expect C1 B2 to be symmetric, If they are truly substitude. How can I enforce this symmetry when creating a system of

Re: [R] How to install packages without Tcl/Tk interface

2009-07-23 Thread Romain Francois
On 07/23/2009 08:58 AM, Matej Kovacic wrote: Hi, how to install packages without Tcl/Tk interface? For instance, if I run install.packages('R.oo'), I need to select Cran repository via graphic interface. How to avoid this? bye, Matej Use the repos argument of install.packages:

Re: [R] A question on operation on list

2009-07-23 Thread Gavin Simpson
On Wed, 2009-07-22 at 14:07 -0700, megh wrote: Thanks for your suggestions. I need one more thing : x = y = vector(list) for (i in 1:5) x[[i]] = rnorm(2); y[[i]] = rnorm(2) Here I want to get t(x[[i]]) %*% y[[i]] for each i. Can anyone please help me? Two ways: set.seed(123) x = y =

[R] Odp: simple question about exporting data...

2009-07-23 Thread Petr PIKAL
Hi Strange. Before starting with R about 10 years ago I used Excel for for graphing too. As soon as I learned few tricks and read intro documents I abandoned Excel for graphing almost completely. Now even for simple visualisations of data I copy them ***from*** spreadsheet to R. But as you

Re: [R] Network from package functions

2009-07-23 Thread Romain Francois
Hi, This post might get you started: http://romainfrancois.blog.free.fr/index.php?post/2009/03/07/What-functions-are-called-by-my-function tail( callees( nls ), 20 ) functions names new.env nls.control nlsModel 812

Re: [R] param:qsec in ggplot2

2009-07-23 Thread ONKELINX, Thierry
Dear Mo You only need to use the name of a variable in your dataframe. The example below lets the size of the points depend on the variable carat. library(ggplot2) ggplot(diamonds, aes(x = table, y = price, size = carat)) + geom_point() HTH, Thierry

Re: [R] param:qsec in ggplot2

2009-07-23 Thread Matthieu Dubois
Hi, qseq is a variable from the data.frame mtcars, that has been used to create the initial plot p. library(ggplot2) p - ggplot(mtcars, aes(wt, mpg)) p + geom_point(aes(size = qsec) Just replace qseq by your variable of interest. Regards, Matthieu

[R] xyplot axis scaling with dates/times

2009-07-23 Thread David Carslaw
Dear R users, I'm trying to get a good x-scale and labels on a plot like the one below. library(lattice) ## make almost a year of hourly data: mydat - data.frame(dates = Sys.time() + 3600 * (1:7000), y = runif(7000)) ## plot it xyplot(y ~ dates, data = mydat, type = l) Only one x-label is

Re: [R] contour plot

2009-07-23 Thread axionator
That problem is that for every y I get a different x and thus cant create an array as in the help page examples. I was just wondering, whether there is a possibility to create such a contour plot easily without having to spend too much time with setting up a suitable matrix manually. Armin

Re: [R] Antw: How to extract the upper xlim and ylim of my plot?

2009-07-23 Thread Uwe Ligges
Frank Bloos wrote: You may want to use the function corner.label from the plotrix-package. Frank Mark Na mtb...@gmail.com 21.07.2009 23:03 Dear R-helpers, I wish to place some text in a plot, at approx 10% of my upper xlim and approx 90% of my upper ylim, i.e.

Re: [R] contour plot

2009-07-23 Thread Petr PIKAL
Hi maybe package akima with function interp can help you Regards Petr r-help-boun...@r-project.org napsal dne 23.07.2009 11:04:25: That problem is that for every y I get a different x and thus cant create an array as in the help page examples. I was just wondering, whether there is a

Re: [R] contour plot

2009-07-23 Thread ONKELINX, Thierry
Have a look at the reshape and ggplot2 package. ds - matrix(rnorm(100), nrow = 10) library(reshape) molten - melt(data = ds) library(ggplot2) ggplot(molten, aes(x = X1, y = X2, z = value)) + geom_contour() HTH, Thierry

Re: [R] What works on windows

2009-07-23 Thread Uwe Ligges
David Winsemius wrote: On Jul 22, 2009, at 3:02 PM, malcolm Crouch wrote: Hi , How can you tell which packages work on windows ? This is the place to find listing for automated compiler check results: http://cran.stat.ucla.edu/bin/windows/contrib/checkSummaryWin.html Yes, and

[R] simple lme question

2009-07-23 Thread Buckmaster, Sarah
Hi everyone, I am trying to analyse my data from a small plant experiment (for a meeting tomorrow afternoon) and am a beginner to R so I apologise if this is a very basic question. I carried out a plant experiment examining plant interactions between two species (A and B) under different

[R] source compile and problem with Hmisc

2009-07-23 Thread Stefan
Hey, I want to compile R version 2.9.1 on debian stable but the make command stops with this message: ... Building/Updating help pages for package 'tcltk' Formats: text html latex example make[2]: Leaving directory `/usr/src/R/R-2.9.1/src/library' make[1]: Leaving directory

Re: [R] source compile and problem with Hmisc

2009-07-23 Thread Uwe Ligges
Stefan wrote: Hey, I want to compile R version 2.9.1 on debian stable but the make command stops with this message: ... Building/Updating help pages for package 'tcltk' Formats: text html latex example make[2]: Leaving directory `/usr/src/R/R-2.9.1/src/library' make[1]: Leaving directory

Re: [R] source compile and problem with Hmisc

2009-07-23 Thread smu
On Thu, Jul 23, 2009 at 12:16:36PM +0200, Stefan wrote: I want to compile R version 2.9.1 on debian stable but the make command stops with this message: ... begin installing recommended package VR Error in library(Hmisc, verbose = FALSE) : there is no package called 'Hmisc' Execution

[R] Counter

2009-07-23 Thread amor Gandhi
Hi everyone,   Is there any counter function in R for the following purpose: x - matrix(c(1,1,0,2,1,0,0,2,0,1,2,1,2,1,0,1),nrow=4) As I would like to know how many zeros, ones, and twos in each row of x?   Many thank in advance, Amor     [[alternative HTML version deleted]]

Re: [R] Counter

2009-07-23 Thread Uwe Ligges
apply(x, 1, table) Uwe Ligges amor Gandhi wrote: Hi everyone, Is there any counter function in R for the following purpose: x - matrix(c(1,1,0,2,1,0,0,2,0,1,2,1,2,1,0,1),nrow=4) As I would like to know how many zeros, ones, and twos in each row of x? Many thank in advance, Amor

Re: [R] Counter

2009-07-23 Thread Philipp Pagel
On Thu, Jul 23, 2009 at 10:44:16AM +, amor Gandhi wrote: Hi everyone, Is there any counter function in R for the following purpose: x - matrix(c(1,1,0,2,1,0,0,2,0,1,2,1,2,1,0,1),nrow=4) As I would like to know how many zeros, ones, and twos in each row of x? Yes there is: As you are

[R] Non-negative solutions to complicated equations

2009-07-23 Thread stenort
Hi all, I have a system of 3 equations with many defined parameters and 3 variables I need to find solutions to. I actually know the solutions I'm aiming for (.0061,0.287,0.070) but R tends to give me (0,0,0). I tried the BB package but don't really follow how to refine my solutions from that;

[R] R interface to PSwarm

2009-07-23 Thread Hans W Borchers
Dear list: Being a bit unsatisfied with global optimization approaches in R such as SANN (in 'optim') or DEoptim, I looked for alternatives on the Web, such as PSwarm http://www.norg.uminho.pt/aivaz/pswarm/ or PIKAIA. There is an R interface to PSwarm (version 1.4) on its home page which I was

Re: [R] Re ading Image Files

2009-07-23 Thread Paul Hiemstra
cotixan wrote: Hello, I'm rather new to R and I want to do some image analysis. Is there a way to read jpeg files into a matrix like matlab's imread? http://cran.r-project.org/web/packages/rimage/index.html -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences

[R] Another question on gplot

2009-07-23 Thread RON70
I am getting following strange error : library(sn) library(ggplot2) dat1 - as.matrix(rsn(1000, 0, 1, 0))[,1] ggplot() + geom_histogram(aes(x = dat1, y = ..density.., fill = ..count..)) + xlab(Distribution) + scale_y_continuous() Error: No data for layer Can anyone please tell

[R] SweaveOpts(eval=false) not working

2009-07-23 Thread Rainer M Krug
Hi I hope this is the right mailing list - if not, could you please refer me to a mora appropriate? My question: I am using sweave (in LyX with beamer) for a lecture and I would like to evaluate the R code chunks for the handout, but not for the presentation. I thought that I can use

Re: [R] Another question on gplot

2009-07-23 Thread ONKELINX, Thierry
Well, you did not specify the dataset. You want: ds - data.frame(dat1 = rsn(1000, 0, 1, 0)) ggplot(ds, aes(x = dat1)) + geom_histogram(aes(y = ..density.., fill = ..count..)) + xlab(Distribution) + scale_y_continuous() HTH, Thierry

Re: [R] SweaveOpts(eval=false) not working

2009-07-23 Thread Duncan Murdoch
On 23/07/2009 8:58 AM, Rainer M Krug wrote: Hi I hope this is the right mailing list - if not, could you please refer me to a mora appropriate? My question: I am using sweave (in LyX with beamer) for a lecture and I would like to evaluate the R code chunks for the handout, but not for the

Re: [R] SweaveOpts(eval=false) not working - CORRECTION

2009-07-23 Thread Rainer M Krug
On Thu, Jul 23, 2009 at 3:04 PM, Duncan Murdochmurd...@stats.uwo.ca wrote: On 23/07/2009 8:58 AM, Rainer M Krug wrote: Hi I hope this is the right mailing list - if not, could you please refer me to a mora appropriate? My question: I am using sweave (in LyX with beamer) for a lecture and

[R] RMySQL Error

2009-07-23 Thread Brock Tibert
Hi Everyone, I am relatively new to R, but I want to try to learn how to use it and potentially replace the other commercial packages I have learned to use over the years. That said, I need to get at data in databases, more specifically, MySQL. I installed the RMySQL package without error, but

Re: [R] How to list R object properties save workspace?

2009-07-23 Thread gug
Hi Matej, You can also try: sapply(ls(), function(x) object.size(get(x))) or eapply(.GlobalEnv, object.size) which list all objects - as with ls() - but giving their sizes. I'm also quite new to R so am not sure which other properties you could hope to get out of it, but by substituting

Re: [R] RMySQL Error

2009-07-23 Thread Cedrick Johnson
Make a copy of the libmySQL.dll and put it in your R_HOME\bin directory and try again.. Just to be on the safe side, I created an env variable called R_HOME and added R_HOME\bin to my path. HTH -c Brock Tibert wrote: Hi Everyone, I am relatively new to R, but I want to try

Re: [R] simple lme question

2009-07-23 Thread Ben Bolker
Sarah B wrote: Hi everyone, I am trying to analyse my data from a small plant experiment (for a meeting tomorrow afternoon) and am a beginner to R so I apologise if this is a very basic question. I carried out a plant experiment examining plant interactions between two species (A

Re: [R] Counter

2009-07-23 Thread Jorge Ivan Velez
Try also: t( apply(x, 1, function(x) table( factor(x, levels = 0:2) ) ) ) # 0 1 2 # [1,] 1 2 1 # [2,] 1 3 0 # [3,] 3 0 1 # [4,] 0 2 2 See ?apply, ?table and ?factor for more details, examples and information. HTH, Jorge On Thu, Jul 23, 2009 at 6:44 AM, amor Gandhi amorigan...@yahoo.de

[R] ROCR - confidence interval for Sens and Spec

2009-07-23 Thread Rasanga Ruwanthi
Dear List,   I am new to ROC analysis and the package ROCR. I want to compute the confidence intervals of sensitivity and specificity for a given cutoff value. I have used the following to calculate sensitivity and specificity:   data(ROCR.simple) pred - prediction(ROCR.simple$predictions,

[R] How to get w in a SVR_package e1071

2009-07-23 Thread marlene marchena
Hi R user, I need some help about how to calculate w in a SVR in package e1071. I have a regression y_i=f(x_i)+e where f(*x*)=(w,phi(x))+b then go on with the SVR calculation I know that w*=Sum_i=1^n [(á_i - á*_i)K(x,x_i) ] where á_i and á*_i are the lagrangian multipliers of the dual form.

Re: [R] simple question about exporting data...

2009-07-23 Thread Greg Snow
Doing the computations in R then the graphs in Excel reminds me of the maxim: Measure with a micrometer Mark with chalk Cut with an ax Before continuing you should really read: http://www.burns-stat.com/pages/Tutor/spreadsheet_addiction.html and look at the final graphic in:

[R] exactly overlaid semi-transparent lines

2009-07-23 Thread Daniel Farewell
I'd like to plot n (say n = 10) semi-transparent lines in such a way that if all n happen to exactly overlay each other, the resulting line is completely opaque. In principle, I believe that this is what setting alpha = 1/n should accomplish. In practice, different values of n produce

[R] remove multiple columns by name from dataframe

2009-07-23 Thread Anne Skoeries
Hi there, I'm trying to remove multiple columns by name from a data.frame. As a result I need to get back the modified data.frame without the removed columns. My columns I want to delete are listed in a vector called delete. data - read.csv2(data.csv) delete - c(col1, col2, col3) newData

[R] rows missing after dataset loaded to R

2009-07-23 Thread Rnewbie
Dear all, When I loaded a dataset (a txt file), which is structured in a tabular format, to R by using read.delim, I found some rows were missing. The column number was correct. These missing rows are no apparently different from the other rows, and for some unknown reasons these missing rows

[R] constrained OLS again

2009-07-23 Thread Bastian.Offermann
Hello, i have a short question concerning constrained OLS regressions using ?lm . i would like the value of my only independent variable to be constrained to 1. ## data capg - rnorm(100) invest - rnorm(100) ## regression reg- lm(capg ~ invest) how do i use ?offset in that regard?

Re: [R] remove multiple columns by name from dataframe

2009-07-23 Thread Jorge Ivan Velez
Dear Anne, Check out http://www.nabble.com/correct-way-to-subset-a-vector-to24412577.html#a24412577 for alternatives. HTH, Jorge On Thu, Jul 23, 2009 at 9:49 AM, Anne Skoeries h...@anne-skoeries.dewrote: Hi there, I'm trying to remove multiple columns by name from a data.frame. As a

Re: [R] remove multiple columns by name from dataframe

2009-07-23 Thread Sarah Goslee
There are several ways to do it. Here's one: fakedata - data.frame(col1 = runif(10), col2 = runif(10), col3 = runif(10)) delete - c(col1, col3) fakedata2 - fakedata[, !(colnames(fakedata) %in% delete), drop=FALSE] Sarah On Thu, Jul 23, 2009 at 9:49 AM, Anne Skoeriesh...@anne-skoeries.de wrote:

Re: [R] remove multiple columns by name from dataframe

2009-07-23 Thread jim holtman
try this: data - data[, !(names(data) %in% delete)] On Thu, Jul 23, 2009 at 9:49 AM, Anne Skoeriesh...@anne-skoeries.de wrote: Hi there, I'm trying to remove multiple columns by name from a data.frame. As a result I need to get back the modified data.frame without the removed columns. My

Re: [R] rows missing after dataset loaded to R

2009-07-23 Thread jim holtman
try using the options: quote='', comment.char='' You might have a comment character (#) or unbalanced quote marks. Look closely at the rows that your missing, especially several rows prior to them. On Thu, Jul 23, 2009 at 9:35 AM, Rnewbiexua...@yahoo.com wrote: Dear all, When I loaded a

[R] goodfit() in vcd package: computation of chi-squared

2009-07-23 Thread Boris.Vasiliev
I have troubles understanding how goodfit() function in the vcd package computes the Pearson coefficient. Can anybody provide more information on the computation? In particular, for HorseKicks data in vcd package, goodfit() yields oo - goodfit(HorseKicks,type=poisson,method=MinChisq)

Re: [R] Automatic differentiation in R

2009-07-23 Thread Prof. John C Nash
Gabor G. wrote R does not currently have AD (except for the Ryacas package which can do true AD for certain simple one line functions, i..e. input the function and output a function representing its derivative); however, for specific problems one can get close using deriv and

Re: [R] Find multiple elements in a vector

2009-07-23 Thread Philip Twumasi-Ankrah
This should work which((x==2)|(x==3)) --Quotable Quotes- A Smile costs Nothing But Rewards Everything   - Anonymous Happiness is not perfected until it is shared  

Re: [R] remove multiple columns by name from dataframe

2009-07-23 Thread Anne Skoeries
That works perfekt! Thanks so much! -- Anne Skoeries __ 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-guide.html and provide commented, minimal,

Re: [R] simple question about exporting data...

2009-07-23 Thread Marc Schwartz
On Jul 23, 2009, at 8:59 AM, Greg Snow wrote: Doing the computations in R then the graphs in Excel reminds me of the maxim: Measure with a micrometer Mark with chalk Cut with an ax Definitely a fortunes candidate... Marc Schwartz __

[R] setting up LMER for repeated measures and how do I get a p value for my fixed effect, group?

2009-07-23 Thread John Sorkin
R 2.8.1 Windows XP I am trying to analyze repeated measures data (the data are listed at the end of this Email message) and I need help to make sure that I have properly specified my model, and would like to know why lmer does not return a p value for Group, my fixed effect. My subjects are

[R] help with randomisation test...

2009-07-23 Thread Anne Kempel
Dear R-people, I hope asking this is not too cheeky, but I do have a R Problem. I hope that some of you like to play around with R and can help me. Its like this. I have several plant species (A,B,C) and 10 replicates per species. 5 plants per species are damaged, 5 not. I let a caterpillar

Re: [R] Find multiple elements in a vector

2009-07-23 Thread Keith Jewell
Try which(x %in% y) HTH KJ Philip Twumasi-Ankrah nana_kwadwo_der...@yahoo.com wrote in message news:519822.56196...@web39502.mail.mud.yahoo.com... This should work which((x==2)|(x==3)) --Quotable Quotes- A Smile costs Nothing But Rewards Everything - Anonymous

Re: [R] Automatic differentiation in R

2009-07-23 Thread Ben Bolker
nashjc wrote: Gabor G. wrote R does not currently have AD (except for the Ryacas package which can do true AD for certain simple one line functions, i..e. input the function and output a function representing its derivative); however, for specific problems one can get close

Re: [R] Non-negative solutions to complicated equations

2009-07-23 Thread Ravi Varadhan
Hi, I found at least 2 typos in your function. Can you check your equations carefully to make sure that they are correct and then send me a reproducible version of the code that I can cut and paste into my R session? Ravi.

Re: [R] Split plot analysis problems

2009-07-23 Thread Jean-Paul Maalouf
Thanks Mark and Richard for your propositions on how to extract residuals. However, I've tried both solutions on my model, and I got different residuals : If we consider the within residuals : Mark's solution gave me a vector of 24 residuals :

Re: [R] Non-negative solutions to complicated equations

2009-07-23 Thread stenort
Hi Ravi, apologies for the typos. I believe these have been corrected, please see the new version of the question on the forum page Ravi Varadhan wrote: Hi, I found at least 2 typos in your function. Can you check your equations carefully to make sure that they are correct and then send

Re: [R] Automatic differentiation in R

2009-07-23 Thread Gabor Grothendieck
On Thu, Jul 23, 2009 at 11:07 AM, Ben Bolkerbol...@ufl.edu wrote: nashjc wrote:   Gabor G. wrote    R does not currently have AD (except for the Ryacas package    which can do true AD for certain simple one line functions, i..e.    input the function and output a function representing its

[R] error message: .Random.seed is not an integer vector but of type 'list'

2009-07-23 Thread Jim Bouldin
I'm trying to run this simple random sample procedure and keep getting the error message shown. I don't understand this; I've designated x as a numeric vector, so what is going on here? Thanks. x = as.vector(c(1:12));x [1] 1 2 3 4 5 6 7 8 9 10 11 12 mode(x) [1] numeric sample(x, 3)

[R] tweak a bit the legend wiht ggplot2

2009-07-23 Thread Benoit Boulinguiez
Hi all, I've this ggplot2 graph (http://www.4shared.com/file/120101043/f3e9350/isotherm.html). I'd like to tweak the legend for the adsorbent part. I'd like to have full square of the colors instead of the full circles. I guess that I've to work on the scale_color_hue or fill somewhere but I

[R] Allan Variance in R

2009-07-23 Thread Javi Hidalgo
Dear R users, I really like how R works :) I would like to know if there is a way to work with the Allan Variance in R. http://en.wikipedia.org/wiki/Allan_variance In general, if there is a package for working in Signal Analysis. Thank you, Javier.

[R] Allan Variance in R

2009-07-23 Thread Javi Hidalgo
Dear R users, I really like how R works :) I would like to know if there is a way to work with the Allan Variance in R. http://en.wikipedia.org/wiki/Allan_variance In general, if there is a package for working in Signal Analysis. Thank you, Javier.

Re: [R] error message: .Random.seed is not an integer vector but of type 'list'

2009-07-23 Thread Uwe Ligges
Jim Bouldin wrote: I'm trying to run this simple random sample procedure and keep getting the error message shown. I don't understand this; I've designated x as a numeric vector, so what is going on here? Thanks. x = as.vector(c(1:12));x [1] 1 2 3 4 5 6 7 8 9 10 11 12 mode(x)

[R] Changing ts times to dates

2009-07-23 Thread Chuse chuse
Dear all, Ive just started with R and I have question: how can you change a time from a ts object .i.e 2009.004 to 2009-01-01? Is there any function for this? I tried around with as.Date... but it hasnt worked. Thank you beforehand. Chuse. [[alternative HTML version deleted]]

Re: [R] error message: .Random.seed is not an integer vector but of type 'list'

2009-07-23 Thread Jim Bouldin
Thank you. However, when I tried that, I got this message: Warning message: In rm(.Random.seed) : variable .Random.seed was not found Jim Bouldin wrote: I'm trying to run this simple random sample procedure and keep getting the error message shown. I don't understand this; I've

Re: [R] error message: .Random.seed is not an integer vector but of type 'list'

2009-07-23 Thread Uwe Ligges
Jim Bouldin wrote: Thank you. However, when I tried that, I got this message: Warning message: In rm(.Random.seed) : variable .Random.seed was not found In that case, have you attached some package that has its own .Random.seed? Try to find where the current .random.seed comes from R

[R] Activation Functions in Package Neural

2009-07-23 Thread Brigid Mooney
Hi, I am trying to build a VERY basic neural network as a practice before hopefully increasing my scope. To do so, I have been using package neural and the MLP related functions (mlp and mlptrain) within that package. So far, I have created a basic network, but I have been unable to change the

Re: [R] Changing ts times to dates

2009-07-23 Thread Gabor Grothendieck
You may need to modify this slightly depending on the specifics of the fraction calculation but try this to get started: x - 2009.004 dt - as.Date(paste(floor(x) + 0:1, 1, 1, sep = -)) dt[1] + as.numeric(diff(dt)) * (x %% 1) - 1 # 2009-01-01 On Thu, Jul 23, 2009 at 11:48 AM, Chuse

[R] Early-Bird Deadline, Schedule: Data Mining Conference in California - Salford Systems

2009-07-23 Thread Lisa Solomon
Salford Data Mining Conference August 23rd - 25th, 2009 San Diego, California Early-Bird Registration Deadline: July 24th, 2009 http://SalfordDataMining.com/agenda.php Post-Conference Training August 26th - 28th, 2009 Our conference will offer prominent keynote speakers and informative real

Re: [R] ROCR - confidence interval for Sens and Spec

2009-07-23 Thread Frank E Harrell Jr
Rasanga Ruwanthi wrote: Dear List, I am new to ROC analysis and the package ROCR. I want to compute the confidence intervals of sensitivity and specificity for a given cutoff value. I have used the following to calculate sensitivity and specificity: data(ROCR.simple) pred -

Re: [R] error message: .Random.seed is not an integer vector but of type 'list'

2009-07-23 Thread Jim Bouldin
Jim Bouldin wrote: Thank you. However, when I tried that, I got this message: Warning message: In rm(.Random.seed) : variable .Random.seed was not found In that case, have you attached some package that has its own .Random.seed? Try to find where the current .random.seed

Re: [R] error message: .Random.seed is not an integer vector but

2009-07-23 Thread Ted Harding
On 23-Jul-09 15:30:05, Jim Bouldin wrote: I'm trying to run this simple random sample procedure and keep getting the error message shown. I don't understand this; I've Thanks. x = as.vector(c(1:12));x [1] 1 2 3 4 5 6 7 8 9 10 11 12 mode(x) [1] numeric sample(x, 3) Error in

Re: [R] error message: .Random.seed is not an integer vector but

2009-07-23 Thread Ted Harding
On 23-Jul-09 16:08:17, Jim Bouldin wrote: Jim Bouldin wrote: Thank you. However, when I tried that, I got this message: Warning message: In rm(.Random.seed) : variable .Random.seed was not found In that case, have you attached some package that has its own .Random.seed? Try to find

[R] read.csv greater than

2009-07-23 Thread stephen sefick
I have a csv file that I am trying to read in and know that values 0 are erroneous - is there a way to read only value grater than 0. thanks, -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us

Re: [R] read.csv greater than

2009-07-23 Thread stephen sefick
here is the solution x - subset(read.csv(foo.csv), VALUE0) On Thu, Jul 23, 2009 at 11:28 AM, stephen sefickssef...@gmail.com wrote: I have a csv file that I am trying to read in and know that values 0 are erroneous - is there a way to read only value grater than 0. thanks, -- Stephen

[R] using k-means clustering in conjunction with heatmap.2 function

2009-07-23 Thread Suzanne Matthews
Hello, I am trying to create a heatmap that clusters based on a k-means scheme rather than a hierarchical clustering scheme. Suppose I have the following input data, located in sample.table: x1 x2 x3 x4 x1 17.198 16.306 16.806 16.374 x2 14.554 10.866 15.780 14.596 x3 14.374 14.118 14.569 17.352

Re: [R] help with randomisation test...

2009-07-23 Thread Greg Snow
I am not sure that I fully understand what all you want to do (and I don't understand why you need the correlation and if a correlation based on 3 pairs is even meaningful), but here is a first stab at what you are trying to do: tmp - Species Control_CR Damage_DR A 10 2 A 9 3 A 7 4 A 9 2 A 8 3

Re: [R] error message: .Random.seed is not an integer vector but

2009-07-23 Thread Jim Bouldin
Thanks much Ted. I actually had just tried what you suggest here before you posted, and resolved the problem. Thanks also for the other tips. I wrote x = as.vector(c(1:12)) because I thought that the mode of x might be the problem, the error message pointing to .Random.seed notwithstanding.

[R] How to get w in SVR with e1071 package

2009-07-23 Thread marlene marchena
Hi all, I need some help about how to calculate w in a SVR in package e1071. I have a regression y_i=f(x_i)+e where f(*x*)=(w,phi(x))+b then go on with the SVR calculation I know that w*=Sum_i=1^n [(á_i - á*_i)K(x,x_i) ] where á_i and á*_i are the lagrangian multipliers of the dual

Re: [R] how to predict dynamic model in R

2009-07-23 Thread Hongwei Dong
Hi, Gabor and Other R users, I'm re-posting my script and the results I got. here is the dynamic model I used to estimate in-sample model (1996-2006) and it works: fit-dyn$lm(Y~lag(Y,-1)+z+x+lag(x,-1)+lag(x,-2)+lag(x,-3)+lag(x,-4)) Then I used this model to do out sample forecast with

Re: [R] error message: .Random.seed is not an integer vector but

2009-07-23 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Jim Bouldin Sent: Thursday, July 23, 2009 9:49 AM To: ted.hard...@manchester.ac.uk; ted.hard...@manchester.ac.uk; r-h...@r- project.org Subject: Re: [R] error message:

Re: [R] simple question about exporting data...

2009-07-23 Thread John Kane
I vote for it. --- On Thu, 7/23/09, Marc Schwartz marc_schwa...@me.com wrote: From: Marc Schwartz marc_schwa...@me.com On Jul 23, 2009, at 8:59 AM, Greg Snow wrote: Doing the computations in R then the graphs in Excel reminds me of the maxim: Measure with a micrometer Mark with

[R] How to pass a character argument which contains expressions to arg.names in barplot?

2009-07-23 Thread jcano
Hi all Can anybody help me with this? I am trying to include in an automatic way the argument in arg.names in a barplot. I generate the labels I want to appear below the bars with a for loop, and they contain subscripts, so I need to use expression anch-0.05 esp-4 for (i in

Re: [R] how to predict dynamic model in R

2009-07-23 Thread Gabor Grothendieck
Please provide your code in a reproducible form (as requested previously). Here we fit with first 9 points and then add a point for prediction. (Of course your model can only predict the current value of Y so you may have to rethink your model even aside from the implementation if you really

[R] Random # generator accuracy

2009-07-23 Thread Jim Bouldin
Dan Nordlund wrote: It would be necessary to see the code for your 'brief test' before anyone could meaningfully comment on your results. But your results for a single test could have been a valid random result. I've re-created what I did below. The problem appears to be with the weighting

[R] howto create a list row-by-row as input to function call?

2009-07-23 Thread Mark Knecht
Hi, I'm having trouble within my function CalcPos to get it to call CalcHorz with values from each row. I *think* it's calling CalcHorz with the final values of the inputs and not the values from each row. How can I do this properly in R? The values aa,bb,cc,dd are inputs. CalcPos first

[R] How to perform a calculation in each element of my list?

2009-07-23 Thread Mark Na
Hi R-helpers, I have a list containing 10 elements, each of which is a dataframe. I wish to add a new column to each list element (dataframe) containing the product of the last two columns of each dataframe. I'd appreciate any pointers, thanks! Mark Na [[alternative HTML version

[R] dimension trouble for a matrix

2009-07-23 Thread Erin Hodgess
Dear R People: I'm having trouble with something that should be very simple. I'm setting up a matrix outside of a loop and writing items into it during the loop. Here is the output: glob3b(sites.info) dim 27 3 [1] /raid1/osg-app Error in xy[i, ] : incorrect number of dimensions Here is the

Re: [R] How to perform a calculation in each element of my list?

2009-07-23 Thread Jorge Ivan Velez
Dear Mark, Try this: lapply(yourlistofdataframes, function(d){ k - ncol(d) d$product - d[,k-1] * d[,k] d } ) HTH, Jorge On Thu, Jul 23, 2009 at 2:05 PM, Mark Na mtb...@gmail.com wrote: Hi R-helpers, I

Re: [R] How to perform a calculation in each element of my list?

2009-07-23 Thread Daniel Malter
Hi Mark, study the following example. The two simulated dataframes are put in a list called listdata. The loop iterates throught the elements in the list and multiplies the last column listdata[[i]][,length(listdata[[i]])] by the column before the last listdata[[i]][,length(listdata[[i]])-1]. This

Re: [R] dimension trouble for a matrix

2009-07-23 Thread Steve Lianoglou
Hi Erin, On Jul 23, 2009, at 2:12 PM, Erin Hodgess wrote: Dear R People: I'm having trouble with something that should be very simple. I'm setting up a matrix outside of a loop and writing items into it during the loop. Here is the output: glob3b(sites.info) dim 27 3 [1] /raid1/osg-app

Re: [R] Random # generator accuracy

2009-07-23 Thread Greg Snow
Try adding replace=TRUE to your call to sample, then you will get numbers closer to what you are expecting. -- 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

  1   2   >