Thanks David for the paper, I understand the theory.
But my question is about R only: the vector of coefficients that R outputs
in lars(), does it apply against the original variable y or against
(y-y_bar). I have put in intercept=T as well in my lars() model.
I need this information to calculate
On 05/05/2013, at 03:05 AM, Gilson Carvalho wrote:
>
> I am trying to perform CCA and RDA several times to compare their
> performances. To do that I put this two functions of vegan package inside a
> for loop. Unfortunatly when I am trying to perform the anova(by=axis) to test
> the individu
Hi All,
I am trying to perform CCA and RDA several times to compare their
performances. To do that I put this two functions of vegan package inside a
for loop. Unfortunatly when I am trying to perform the anova(by=axis) to
test the individual axes the function don't recognize the species
data.fram
On May 4, 2013, at 3:24 PM, Jose Narillos de Santos wrote:
> Hi all,
>
> I want to know if the is a possibility to increase the size of the default
> numbers that appears on axis on a plot, I don´t mean a labbel I mean the
> number on axis and yaxis.
>
> You see when in excel you select a plot
I am running caret for model exploration. I developed my code a number of
months ago and I've been running it with no issues. Recently, I updated my
version of caret however, and now I am getting a new error. I'm wondering
if this is due to the new release.
The error I am getting is when I am runn
Hi all,
I want to know if the is a possibility to increase the size of the default
numbers that appears on axis on a plot, I don´t mean a labbel I mean the
number on axis and yaxis.
You see when in excel you select a plot and select all and increase the
font letters or mark bold all the graph par
On 05/05/2013 04:23 AM, Jose Narillos de Santos wrote:
Hi I have this data and I want to draw a plot of a line BUR the thing is I
want that this line appear smooth ¿someone knows if there is an special
function?
Thanks in advance.
Vertice 02/05/2013 1w 0,083 2w 0,092 3w 0,1 1m 0,113 2m
Hi everyone,
>From the answers I've received to the question below, it occurs to me there
may be more than inefficient programming on my part involved: The 'apply'
code described below quickly takes up the 18 GB of memory I have available,
which leaves my machine functioning at a crawl for the a
You ask for help with smooth lines, but show no indication what you are doing
with this data. Where is your example code that plots this data with not-smooth
lines? I think the HTML email has messed up your data anyway. Next time post
plain text email (See the Posting Guide mentioned at the bott
(I neglected to use reply-all.)
-- Forwarded message --
From: David Romano
Date: Sat, May 4, 2013 at 11:25 AM
Subject: Re: [R] how to parallelize 'apply' across multiple cores on a Mac
To: Charles Berry
On Sat, May 4, 2013 at 9:32 AM, Charles Berry wrote:
> David,
>
> If you i
Hi I have this data and I want to draw a plot of a line BUR the thing is I
want that this line appear smooth ¿someone knows if there is an special
function?
Thanks in advance.
Vertice 02/05/2013 1w 0,083 2w 0,092 3w 0,1 1m 0,113 2m 0,159 3m
0,201 4m 0,236 5m 0,272 6m 0,302 7m 0,333 8
Hi,
You can get the table results in a matrix:
Unel<-sort(unique(unlist(dat1)))
res2<- t(mapply(function(x,y) ifelse(x %in%
names(y),y[as.character(x)],0),list(Unel),lst1))
dimnames(res2)<- list(rownames(dat1),Unel)
head(results)
# A M P
#1 38 0 0
#2 38 0 0
#3 38 0 0
#4 34 2 2
#5 38 0 0
#6
Hi,
I rectified my error (thanks David for pointing it out)
Now I have been able to run the code:
data=read.table("data.txt", header=T)
> l=data$LOSS
> h=data$HPI
> u=data$UE
> g=data$GDP
>
> matrix=cbind(g,h,u)
> lasso=lars(matrix,l)
>
The final set of coefficients for the regression is the last
The offset(x) term means to include x (or whatever is in the parentheses)
into the model as is, without computing a slope for that term. You could
also include offset( 1 * x ) instead which might make this a bit more
explicit (but would not actually make any difference). Since x by itself
is also
Hi,
lst1<-apply(dat1,1,table)
res<-which(unlist(lapply(lst1,function(x) {indx<-names(x)%in%"A"; if(any(indx))
x[indx]<19 else NA}),use.names=FALSE))
head(res)
#[1] 19 39 40 41 43 45
length(res)
#[1] 1481
lst1[c(19,39,40)]
#$`19`
#
# A P
# 1 37
#
#$`39`
#
# A P
#12 26
#
#$`40`
#
# A P
#16
In addition to David's great answer also read
?contr.SAS
?contr.treatment
?relevel
and also sections 4 and 11.1.1 in An Introduction to R.
On Fri, May 3, 2013 at 9:47 PM, David Winsemius wrote:
>
> On May 3, 2013, at 3:32 PM, Iuri Gavronski wrote:
>
> > Hi,
> > I would like to know the criteri
Parameters are different from functions, and offset is a function. Kindly read
the help for that function and the references given there.
---
Jeff NewmillerThe . . Go Live...
DCN:
David Romano stanford.edu> writes:
>
> Hi everyone,
>
> I'm trying to use apply (with a call to zoo's rollapply within) on the
> columns of a 1.5Kx165K matrix, and I'd like to make use of the other cores
> on my machine to speed it up. (And hopefully also leave more memory free: I
> find that a
On May 3, 2013, at 10:58 PM, peter dalgaard wrote:
>
> On May 3, 2013, at 21:36 , David Winsemius wrote:
>
>>
>> On May 3, 2013, at 10:46 AM, peter dalgaard wrote:
>>>
>>>
>>> Because comparison with an unknown value yields an unknown result.
>>
>> Anything else would violate the Second La
On May 4, 2013, at 6:09 AM, Preetam Pal wrote:
> Hi all,
> I have a data set containing variables LOSS, GDP, HPI and UE.
> (I have attached it in case it is required).
>
> Having renamed the variables as l,g,h and u, I wish to run a Lasso
> Regression with l as the dependent variable and all th
Hi,
If it is to get the first 18 genes (based on rownumbers)
dat1<- read.table("apcall.txt",header=TRUE,sep="",stringsAsFactors=FALSE)
dat1[1:18,]
#or
dat1[as.numeric(rownames(dat1))<19,]
A.K.
>The row numbers represent 7129 different genes
>The 38 columns represent 38 different patients
>
Hi,
I was wondering if there is a package/function for Panel Granger
non-causality tests? I am interested in Toda-Yamamoto like procedure for
panel models.
Thank you,
--
View this message in context:
http://r.789695.n4.nabble.com/Panel-Granger-Non-Causality-Tests-in-R-tp4666316.html
Sent fr
Hi,
One more possibility:
dateA<-seq.Date(as.Date("1jan2012",format="%d%b%Y"),as.Date("14Feb2013",format="%d%b%Y"),by="day")
dateB<-seq.Date(as.Date("1Mar2012",format="%d%b%Y"),as.Date("30Nov2012",format="%d%b%Y"),by="day")
set.seed(15)
A<-data.frame(dateA,value=cumsum(rnorm(411)))
set.seed(2
I am trying to use the github_install from the devtool package.
I get the following error on
install_github("roxygen"):
" Peer certificate cannot be authenticated with given CA certificates"
This is probably some github security issue (and not strictly
a devtools thing). Has anyone seen this bef
Hi all,
I have a data set containing variables LOSS, GDP, HPI and UE.
(I have attached it in case it is required).
Having renamed the variables as l,g,h and u, I wish to run a Lasso
Regression with l as the dependent variable and all the other 3 as the
independent variables.
data=read.table("dat
Hello,
Thanks.
But the parameter offset is new to me.
Please kindly explain why setting offset to x will give a significant test
of whether the slope coefficient is different from one.
(I checked the ?lm but still do not understand it well)
Thanks again
Elaine
On Wed, May 1, 2013 at 11:12 AM,
Hello,
I forgot a third possibility
plot(as.ts(cbind(x, y)), plot.type = "single", col = 1:2)
Rui Barradas
Em 04-05-2013 10:37, Rui Barradas escreveu:
Hello,
Where is the reproducible example? See
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
Any
On 04-05-2013, at 11:24, Preetam Pal wrote:
> Hi all,
>
> (Kindly ignore the previous mail.)
>
> I have 2 time series variables : X_t and Y_t where t=0,1,2,...,100
> I want a plot containing the line charts( time plotted on the x-axis) for
> both these variables.
> How can I do this or w
Hello,
Where is the reproducible example? See
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
Anyway, see if this does what you need.
set.seed(1)
x <- cumsum(rnorm(100))
y <- cumsum(rnorm(100))
ylm <- range(c(x, y))
plot(as.ts(x), ylim = ylm)
lines(y, co
Hi all,
(Kindly ignore the previous mail.)
I have 2 time series variables : X_t and Y_t where t=0,1,2,...,100
I want a plot containing the line charts( time plotted on the x-axis) for
both these variables.
How can I do this or where do I look for the information on this?
Appreciate your he
Hi all,
I have 2 time series variables : X_t and Y_t where t=0,1,2,...,100
I want a plot containing the line charts( time plotted on the x-axis) for
both these variables.
--
Preetam Pal
(+91)-9432212774
M-Stat 2nd Year, Room No. N-114
Statistics
31 matches
Mail list logo
