Euna Jeong 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 _training_ data it is in
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 p
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
Sender: r-help-bounces@r-project.orgDate: Thu, 5 Sep 2013 07:01:48
To: John Kane; r-help@R-pro
On 05-09-2013, at 09:01, BLANDENIER Lucien 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
> ___
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 explai
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 541
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
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 o
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 destill
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
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 wrote:
> Hi Spencer,
>
> It really is not very hard, and I have never had issue with it:
>
> http://www.oracle.com/technetwork/java/javase/downloads/jdk7-downloads-188026
On Sep 5, 2013, at 2:54 AM, Daniel Hornung wrote:
On Wed, Sep 4, 2013 at 1:45 PM, Daniel Hornung
wrote:
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 arguments but not all. Yo
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 g
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
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
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' th
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","numer
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
read.
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', s
# 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
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
Co
Hi Michael,
See comments in line.
On Wed, Sep 4, 2013 at 10:18 PM, Michael Hacker 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 engineering f
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:
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(
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' breakpoin
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 lw
> 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.
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 tha
On Thu, Sep 5, 2013 at 11:08 AM, Anna Zakrisson Braeunlich
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:
>
> ggplot(mydata, aes(x = w
Hi Anna,
On Thu, Sep 5, 2013 at 6:13 AM, Anna Zakrisson Braeunlich
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 want a
> funct
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, fill=organ
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, setWi
Thank you everyone
I got actual values in my view
On Thu, Sep 5, 2013 at 7:31 PM, Keith Jewell
wrote:
> 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 spec
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, wrote:
Slight correction:
On Thu, Sep 5, 2013 at 7:48 AM, Bert Gunter 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 venture to guess
HI,
No problem.
I am using linux, may be there is a difference.
You could try this:
for(i in seq_along(lst1))
write.table(lst1[[i]],file="Crish1.txt",append=TRUE,quote=FALSE,sep=",")
#Will get some warning message.
You can save .txt to .csv and open it in excel. The column headings would be
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(x
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
Dear Michael,
Please see comments below, interspersed with your questions:
On Wed, 4 Sep 2013 22:18:57 -0400
"Michael Hacker" wrote:
> Dear Colleagues,
>
>
>
> I'm working on a Delphi study comparing perceptions of high school
> technology teachers and university engineering educators about
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
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 t
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 cumsum(
On 3 September 2013 16:10, Worthington, Thomas A
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 corCAR1(~ID | SiteCode1
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 en
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 cal
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
?sys.parent (and friends)
-- Bert
On Thu, Sep 5, 2013 at 4:37 PM, ivo welch wrote:
> 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 f
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
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 wrote:
> 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
Hi,
May be this helps:
m1<- as.matrix(read.table(text="
y1 g24
y1 0 1
g24 1 0
",sep="",header=TRUE))
m2<-as.matrix(read.table(text="y1 c1 c2 l17
y1 0 1 1 1
c1 1 0 1 1
c2 1 1 0 1
l17 1 1 1 0",sep="",header=TRUE))
m3<- as.matrix(read.table(text="y1 h4 s2 s30
y1 0 1 1 1
h4 1 0 1 1
s2 1
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 nati
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 nati
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 ex
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]}))
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 sensibl
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
https://stat.et
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
Hi Rich
This is a case of where data needs to be supplied to go further.
Regards
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behal
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
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)
text(
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
> I've been working on a way to visualize a spearman correlation. That se
61 matches
Mail list logo
