Re: [R] xyplot and lwd

2013-09-05 Thread Daniel Hornung
On Wed, Sep 4, 2013 at 1:45 PM, Daniel Hornung daniel.horn...@ds.mpg.dewrote: Hello, can it be that xyplot does not support the lwd argument? At least here, the following still shows thin lines, as opposed to the regular plot command: xyplot(Sepal.Length ~ Sepal.Width, data

Re: [R] Question about R2 in pls package

2013-09-05 Thread Bjørn-Helge Mevik
Euna Jeong eaje...@gmail.com writes: I have questions about R2 used in pls (or multivariate analysis). Is R2 same with the square of the PCC (Pearson Correlation Coefficient)? If you read the manual for R2 in the pls package, it will tell you how R2 is calculated there, and that for

Re: [R] Problem with installing the TRR package

2013-09-05 Thread BLANDENIER Lucien
Do you know if R-3.0.1 is available for Linux Mint? Do you know how I can check it? De : John Kane [jrkrid...@inbox.com] Envoyé : mercredi 4 septembre 2013 17:31 À : BLANDENIER Lucien; r-help@R-project.org Objet : RE: [R] Problem with installing the TRR

Re: [R] Problem with installing the TRR package

2013-09-05 Thread szehnder
I do not know Linux Mint, but what is always possible is to build R from sources on a UNIX system. Gesendet über den BlackBerry® Service von E-Plus. -Original Message- From: BLANDENIER Lucien lucien.blanden...@unine.ch Sender: r-help-bounces@r-project.orgDate: Thu, 5 Sep 2013 07:01:48

Re: [R] Problem with installing the TRR package

2013-09-05 Thread Berend Hasselman
On 05-09-2013, at 09:01, BLANDENIER Lucien lucien.blanden...@unine.ch wrote: Do you know if R-3.0.1 is available for Linux Mint? Do you know how I can check it? Look on CRAN: http://cran.r-project.org/bin/linux/ubuntu/ Since Mint is derived from Ubuntu this should work. Berend

[R] sparse PCA using nsprcomp package

2013-09-05 Thread array chip
Hi all, I am using nsprcomp() from nsprcomp package to run sparse PCA. The output is very much like regular PCA by prcomp() in that it provides sdev for standard deviation of principle components (PC). For regular PCA by prcomp(), we can easily calculate the percent of total variance

[R] R Help

2013-09-05 Thread Mª Teresa Martinez Soriano
Hi everone and thanks for this service, I have a dataset which look like: X. IE.2003 IE.2004 IE.2005 IE.2006 IE.2007 IE.2008 IE.2009 IE.2010 14560 118958 187 475 571 76410471203 715 807 12737 105571 935 942 917 948 991 861 NA

Re: [R] R Help

2013-09-05 Thread Pascal Oettli
Hello, Please stop to use R help for the subject of your mail. You already used it several times. You have been asked to stop to send e-mail in HTML. You also have been asked to use dput() when you want to submit data to this list. For the current problem, there is no reproducible code, as it

Re: [R] New Version of R 3.0.1 problems with installing Rcmdr

2013-09-05 Thread Uwe Ligges
On 05.09.2013 04:19, alanidris wrote: I have been using R version 2.15.1 happly along side R Commander. I then tried to go through a fresh install using the latest version of R, R 3.0.1. The trouble started when I wanted to install Rcmdr, I kept getting an error message about previous installs

[R] optim evils

2013-09-05 Thread Michael Meyer
Thanks for all replies. The problem occurred in the following context: A Gaussian one dimensional mixture (number of constituents, locations, variances all unknown) is to be fitted to data (as starting value to or in lieu of mixtools). A likelihood maximization is performed. I'll try to

[R] Capturing warnings with capture.output

2013-09-05 Thread Sebastien Bihorel
Dear R-users, I would like to follow-up on a old thread by Hadley Wickham about capturing warnings with capture.output. My goal is to evaluate a function call and capture the results of the function call plus warning calls and messages if a warning is returned. For instance, in the following

Re: [R] Read a Google Spreadsheet?

2013-09-05 Thread Ista Zahn
I've also never had a problem with both 32 and 64 bit java installed. Best, Ista On Wed, Sep 4, 2013 at 9:50 PM, Joshua Wiley jwiley.ps...@gmail.com wrote: Hi Spencer, It really is not very hard, and I have never had issue with it:

Re: [R] xyplot and lwd

2013-09-05 Thread David Winsemius
On Sep 5, 2013, at 2:54 AM, Daniel Hornung wrote: On Wed, Sep 4, 2013 at 1:45 PM, Daniel Hornung daniel.horn...@ds.mpg.dewrote: Hello, can it be that xyplot does not support the lwd argument? The lattice plotting system uses the grid plotting engine and does accepts some base par-type

Re: [R] xyplot and lwd

2013-09-05 Thread Daniel Hornung
On Thursday, September 05, 2013 13:40:00 David Winsemius wrote: can it be that xyplot does not support the lwd argument? The lattice plotting system uses the grid plotting engine and does accepts some base par-type arguments but not all. You may need to read more about lattice and grid:

Re: [R] Attribute Length Error when Trying plm Regression

2013-09-05 Thread Millo Giovanni
Dear Laura, as Arun said it is difficult to help w/o a reproducible example. However this is most likely to be an indexing problem, as he suggests; the output of traceback() is far from useless here, because it shows that the problem occurs in the data transformation step. The latter, which is by

[R] Problem with converting F to FALSE

2013-09-05 Thread Venkata Kirankumar
Hi, I have a peculier problem in R-Project that is when my CSV file have one column with all values as 'F' the R-Project converting this 'F' to FALSE. Can some one please suggest how to stop this convertion. Because I want to use 'F' in my calculations and show it in screen. for example my data is

Re: [R] Problem with converting F to FALSE

2013-09-05 Thread philippe massicotte
Look at colClasses in ?read.csv Date: Thu, 5 Sep 2013 18:14:49 +0530 From: kiran4u2...@gmail.com To: r-help@r-project.org Subject: [R] Problem with converting F to FALSE Hi, I have a peculier problem in R-Project that is when my CSV file have one column with all values as 'F' the

Re: [R] Problem with converting F to FALSE

2013-09-05 Thread arun
Hi, Try: dat1-read.table(text= sex  group F   1 F   2 F   3 ,sep=,header=TRUE,colClasses=c(character,numeric))  dat1 #  sex group #1   F 1 #2   F 2 #3   F 3 #if you are using read.csv()   dat2-read.csv(new1.csv,sep=,header=TRUE,colClasses=c(character,numeric))  dat2 # 

Re: [R] Problem with converting F to FALSE

2013-09-05 Thread Joshua Wiley
Hi, You can either manually specify colClasses or the asis argument. See ?read.csv for more details. If you just had those two columns, something like: read.table(header = TRUE, text = sex group F 1 T 2 , colClasses = c(character, integer)) Cheers, Josh

[R] Y-axis labels as decimal numbers

2013-09-05 Thread mohan . radhakrishnan
Hi, I am able to create a graph with this code but the decimal numbers are not plotted accurately because the ylim values are not set properly. x-axis is proper. How do I accurately set the ylim for duration.1 column ? Thanks, Mohan set1$duration- as.POSIXct(paste('2013-08-24',

[R] ggplot2: connecting medians of boxes using facet_wrap. Changing median marker.

2013-09-05 Thread Anna Zakrisson Braeunlich
# Dear all, # Thank you for taking your time. # What I would like to do: # (Run the code below to see what I am referring to) # I want lines connecting the medians of of the boxes. I do not want a function, just a simple, # straight connection between points. I also would like the lines to be

[R] Using qmap module

2013-09-05 Thread Jaya Pudashine
Dear all, Is there anybody who is using qmap module for bias correction of RCM data. I have few queries about that. Please let me know if anybody is using qmap module in R. Thanking in Advance -- Jaya Pudashine M.Eng. Water Engineering and Management Asian Institute of Technology, Thailand

Re: [R] Poly Correlations

2013-09-05 Thread Ista Zahn
Hi Michael, See comments in line. On Wed, Sep 4, 2013 at 10:18 PM, Michael Hacker mhac...@nycap.rr.com wrote: Dear Colleagues, I'm working on a Delphi study comparing perceptions of high school technology teachers and university engineering educators about the importance of concepts about

Re: [R] Problem with converting F to FALSE

2013-09-05 Thread Keith Jewell
Depending what you're doing with the data, you might want colClasses=c(factor,numeric) On 05/09/2013 13:58, Joshua Wiley wrote: Hi, You can either manually specify colClasses or the asis argument. See ?read.csv for more details. If you just had those two columns, something like:

[R] plot densities outside axis

2013-09-05 Thread Dustin Fife
I've been working on a way to visualize a spearman correlation. That seemed pretty simple: generate skewed data x = rnorm(100)^2 y = .6*x + rnorm(100, 0, sqrt(1-.6^2)) plot(x,y) regular plot plot(rank(x),rank(y), xaxt=n, yaxt=n) ### spearman-like plot make axis labels axis(1,

Re: [R] Histogram

2013-09-05 Thread David Carlson
I wasn't suggesting that much detail, but I think the addition of one sentence in the last paragraph of the Details section would make it the meaning of the number is a suggestion only clearer. These functions provide a suggested number of bins that may be modified to produce 'round' breakpoints

Re: [R] xyplot and lwd

2013-09-05 Thread Bert Gunter
Daniel: I wondered if that might be what you meant ... To amplify a bit on David's response, the answer is that you do **not** have separate control over the line width of characters -- lwd controls the width of lines in a graph (exactly as it does in base graphics! ), so you misunderstood the

Re: [R] optim evils

2013-09-05 Thread Martin Maechler
Thanks for all replies. The problem occurred in the following context: A Gaussian one dimensional mixture (number of constituents, locations, variances all unknown) is to be fitted to data (as starting value to or in lieu of mixtools). A likelihood maximization is performed. Cool. That

Re: [R] Histogram

2013-09-05 Thread Duncan Murdoch
On 05/09/2013 10:17 AM, David Carlson wrote: I wasn't suggesting that much detail, but I think the addition of one sentence in the last paragraph of the Details section would make it the meaning of the number is a suggestion only clearer. These functions provide a suggested number of bins that

Re: [R] ggplot2: connecting medians of boxes using facet_wrap. Changing median marker.

2013-09-05 Thread Ista Zahn
On Thu, Sep 5, 2013 at 11:08 AM, Anna Zakrisson Braeunlich anna.zakris...@su.se wrote: Hi and thank you for the help and for the fast reply! A. One thing still is a problem. I have prior my first mail tried to not fill the boxes. The result is a different plot. If chosing to do:

Re: [R] ggplot2: connecting medians of boxes using facet_wrap. Changing median marker.

2013-09-05 Thread Ista Zahn
Hi Anna, On Thu, Sep 5, 2013 at 6:13 AM, Anna Zakrisson Braeunlich anna.zakris...@su.se wrote: # Dear all, # Thank you for taking your time. # What I would like to do: # (Run the code below to see what I am referring to) # I want lines connecting the medians of of the boxes. I do not

Re: [R] ggplot2: connecting medians of boxes using facet_wrap. Changing median marker.

2013-09-05 Thread Anna Zakrisson Braeunlich
Hi and thank you for the help and for the fast reply! A. One thing still is a problem. I have prior my first mail tried to not fill the boxes. The result is a different plot. If chosing to do: ggplot(mydata, aes(x = week, y = var1)) instead of ggplot(mydata, aes(x = week, y = var1,

[R] setStatusBar function gives error message in R 3.01 under Windows 7

2013-09-05 Thread Gwen D. LaSelva
I am running R 3.01 under 64-bit Windows 7. When I try to set the status bar, I get an error message. For example: text-hello setStatusBar(text) Error in .Call(setStatusBar, text) : first argument must be a string (of length 1) or native symbol reference The related function,

Re: [R] Problem with converting F to FALSE

2013-09-05 Thread Venkata Kirankumar
Thank you everyone I got actual values in my view On Thu, Sep 5, 2013 at 7:31 PM, Keith Jewell keith.jew...@campdenbri.co.ukwrote: Depending what you're doing with the data, you might want colClasses=c(factor,**numeric) On 05/09/2013 13:58, Joshua Wiley wrote: Hi, You can either

Re: [R] Y-axis labels as decimal numbers

2013-09-05 Thread jim holtman
So what is wrong with the y-axis? When I run your script, things seem right. Can you explain what it is that you want. Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it. On Thu, Sep 5, 2013 at 9:22 AM,

Re: [R] optim evils

2013-09-05 Thread Bert Gunter
Slight correction: On Thu, Sep 5, 2013 at 7:48 AM, Bert Gunter bgun...@gene.com wrote: Michael: Your parameter specification is probably over-determined, so that you have an infinite set of parameter **values** that give essentially the same solution within numerical error. I would

Re: [R] binary symmetric matrix combination

2013-09-05 Thread arun
Also, some of the steps could be reduced by: names1-unique(c(colnames(m1),colnames(m2),colnames(m3),colnames(m4))) Out3-matrix(0,length(names1),length(names1),dimnames=list(names1,names1)) lst1-sapply(paste0(m,1:4),function(x) {x1- get(x); x2-paste0(colnames(x1)[col(x1)],rownames(x1)[row(x1)]);

[R] xyplot() with discontinuous x-axis variable

2013-09-05 Thread Rich Shepard
My xyplot() with superposed multiple condiions looks better with lines than with points (it's easier to see changes over time with the lines). But, there are gaps in the years (the x axis) for which there are data to be plotted. For example, there are data for years 2004-2006 and 2010-2012, but

Re: [R] Poly Correlations

2013-09-05 Thread John Fox
Dear Michael, Please see comments below, interspersed with your questions: On Wed, 4 Sep 2013 22:18:57 -0400 Michael Hacker mhac...@nycap.rr.com wrote: Dear Colleagues, I'm working on a Delphi study comparing perceptions of high school technology teachers and university engineering

Re: [R] ANN: Course Data Mining with R in Spain

2013-09-05 Thread ernesto villarino
We do it tomorrow if you have some time.good night! On Sep 3, 2013 5:04 PM, Soledad De Esteban Trivigno [via R] ml-node+s789695n4675267...@n4.nabble.com wrote: Dear colleague: Registration is open for the course CLASSIFICATION AND REGRESSION TREES AND NEURAL NETWORKS WITH R - Second

[R] Matrix randomization

2013-09-05 Thread levan mumladze
Dear friends I am very new to R and no experience in programming at all. I need to generate random swaps of binary matrix in a way that row and column totals remain constant. and for each derived matrix calculate site weighted richness. Could you suggest most appropriate function (package) to do

Re: [R] sparse PCA using nsprcomp package

2013-09-05 Thread Christian Sigg
Hi John I am currently traveling and have sporadic net access, I therefore can only answer briefly. It's also quite late, I hope what follows still makes sense... For regular PCA by prcomp(), we can easily calculate the percent of total variance explained by the first k PCs by using

Re: [R] Assessing temporal correlation in GAM with irregular time steps

2013-09-05 Thread Gavin Simpson
On 3 September 2013 16:10, Worthington, Thomas A thomas.worthing...@okstate.edu wrote: Dear Gavin Thank you for the very detailed response. I had started to go down the route of fitting a correlation structure via gamm. I tried applying your code to my data but returned the error Error in

[R] get syntax for inherited environments

2013-09-05 Thread ivo welch
quick question. how do I search up the calling environments until I find a variable? x=function() { m=22; y() } y=function() { z() } z=function() { mget(m, inherits=TRUE, ifnotfound=m not found) } x() $m [1] m not found from the perspective of z(), function x is not an enclosing

Re: [R] sparse PCA using nsprcomp package

2013-09-05 Thread array chip
HI Christian, Thanks so much for the detailed explanation! I look forward to the new release of nsprcomp package! At the meantime, I will use the function below for calculation of adjusted standard deviation. I have 2 more questions, hope you can shed some lights on: 1). Assume now I can

Re: [R] binary symmetric matrix combination

2013-09-05 Thread arun
HI, No problem. I think you didn't run the `vecOut` after adding the new matrix.  `lst1` is based on `vecOut` For example: m5- as.matrix(read.table(text=y1 e6 l16  y1 0 1 1 e6 1 0 1 l16 1 1 0,sep=,header=TRUE)) names1-unique(c(colnames(m1),colnames(m2),colnames(m3),colnames(m4), colnames(m5)))

Re: [R] get syntax for inherited environments

2013-09-05 Thread Bert Gunter
?sys.parent (and friends) -- Bert On Thu, Sep 5, 2013 at 4:37 PM, ivo welch ivo.we...@anderson.ucla.eduwrote: quick question. how do I search up the calling environments until I find a variable? x=function() { m=22; y() } y=function() { z() } z=function() { mget(m, inherits=TRUE,

Re: [R] Y-axis labels as decimal numbers

2013-09-05 Thread Jim Lemon
On 09/05/2013 11:22 PM, mohan.radhakrish...@polarisft.com wrote: Hi, I am able to create a graph with this code but the decimal numbers are not plotted accurately because the ylim values are not set properly. x-axis is proper. How do I accurately set the ylim for duration.1 column

Re: [R] Looping an lapply linear regression function

2013-09-05 Thread Flavio Barros
Hello Arun. Can you provide some data? To help you better i will need a complete reproducible example ok? On Thu, Sep 5, 2013 at 1:49 PM, arun smartpink...@yahoo.com wrote: HI, May be this helps: set.seed(28) dat1-

Re: [R] setStatusBar function gives error message in R 3.01 under Windows 7

2013-09-05 Thread Duncan Murdoch
On 05/09/2013 10:20 AM, Gwen D. LaSelva wrote: I am running R 3.01 under 64-bit Windows 7. When I try to set the status bar, I get an error message. For example: text-hello setStatusBar(text) Error in .Call(setStatusBar, text) : first argument must be a string (of length 1) or native

Re: [R] setStatusBar function gives error message in R 3.01 under Windows 7

2013-09-05 Thread Duncan Murdoch
On 05/09/2013 10:20 AM, Gwen D. LaSelva wrote: I am running R 3.01 under 64-bit Windows 7. When I try to set the status bar, I get an error message. For example: text-hello setStatusBar(text) Error in .Call(setStatusBar, text) : first argument must be a string (of length 1) or native

[R] applying a univariate function for each response in a multivariate linear model (mlm)

2013-09-05 Thread Michael Friendly
After fitting a multivariate linear model (mlm), I'd like to be able to run or apply a standard univariate stats:::*.lm function to each of the response variables, within a function -- i.e., by operating on the mlm object, rather than re-running the univariate models separately manually. An

Re: [R] Looping an lapply linear regression function

2013-09-05 Thread arun
HI, May be this helps:  set.seed(28)  dat1- setNames(as.data.frame(matrix(sample(1:40,10*5,replace=TRUE),ncol=5)),letters[1:5]) indx-as.data.frame(combn(names(dat1),2),stringsAsFactors=FALSE) res-t(sapply(indx,function(x) {x1-cbind(dat1[x[1]],dat1[x[2]]);summary(lm(x1[,1]~x1[,2]))$coef[,4]}))  

Re: [R] optim evils

2013-09-05 Thread Bert Gunter
Michael: Your parameter specification is probably over-determined, so that you have an infinite set of parameters that give essentially the same solution within numerical error. I would venture to guess that this will not be fixable with alternative optimizers. It is up to you to provide a

Re: [R] Question about R2 in pls package

2013-09-05 Thread Euna Jeong
Thank you David, it is a good resouce to find related subjects. I need time to understand the topics. Thank you again Bjørn-Helge for your response. Euna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] read.dta()

2013-09-05 Thread Debasish Roy
I've been using R 3.0.1 version. I tried to read a file named  abc.dta() I used the command  X - read.dta(abc.dta) and it gave me             Error: could not find function read.dta Can anyone help me what could be the problem and how to fix it ? Thanks, Deb. [[alternative HTML

Re: [R] read.dta()

2013-09-05 Thread Christopher W. Ryan
I don't know about 3.0.1, but the 2.15.x that I'm still using requires the foreign package--that's where the read.dta command resides. library(foreign) --Chris Ryan SUNY Upstate Medical University Binghamton, NY USA Debasish Roy wrote: I've been using R 3.0.1 version. I tried to read a file

Re: [R] Y-axis labels as decimal numbers

2013-09-05 Thread mohan . radhakrishnan
Hi, set1$duration.2- as.POSIXct(paste('2013-08-24', set1$duration)) plot(set1$duration.2,set1$duration.1,type=b,col = blue, ylab=, xaxt = 'n', xlab=,las=2,lwd=2.5, lty=1,cex.axis=2.5) # now plot you times axis(1, at = set1$duration.2, labels = set1$duration, las = 2,cex.axis=2.5)

Re: [R] plot densities outside axis

2013-09-05 Thread Pascal Oettli
Hello, Using a web search engine, I found, for example: http://www.unt.edu/benchmarks/archives/2003/february03/rss.htm http://sas-and-r.blogspot.jp/2012_09_01_archive.html Hope this helps, Pascal 2013/9/5 Dustin Fife df...@ou.edu I've been working on a way to visualize a spearman