9902, Adjusted R-squared: 0.9902
F-statistic: 1.009e+05 on 1 and 998 DF, p-value: < 2.2e-16
--
Agustin Lobo
aloboa...@gmail.com
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do re
ot found" from within system()?
Any help much appreciated,
Thanks
--
Agustin Lobo
aloboa...@gmail.com
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
76,0.46,0.85,0.48,0.78,0.87,0.22,0.51,0.03,0.39,0.44)
midata <- data.frame(treatments,people,y)
midata[1:10,]
I would do as follows:
Test among treatments
friedman.test(y ~ treatments|people,data=midata)
Test among subjects
friedman.test(y ~ people|treatments,data=midata)
Is this correct?
Than
ects p1 and p2 within one single object (as
in practice I have many objects generated within a for() loop).
The following used to work:
ps <- c(p1,p2)
grid.arrange(ps,ncol=1,nrow=2)
but does not work any more.
How should I combine p1 and p2 into one single object that would be
accepted by grid.arran
h the same length as delme2 so that I can put it in a
data.frame?
Thanks!
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: [EMAIL PROTECTED]
http://www.ija.csic.es
8 12 16
I can do:
> mitest[cbind(c(1,2,3),c(2,3,2))]
[1] 5 10 7
but using the names does not work:
> mitest[cbind(c("a","b","c"),c("b","c","b"))]
[1] NA NA NA NA NA NA
despite the fact that
> mitest["a","b"]
lazuanna wrote:
Try this:
mapply(function(x, y)mitest[x, y], c("a", "b", "c"), c("b", "c", "b"))
or
diag(`[`(mitest, i = c("a", "b", "c"), j = c("b", "c", "b")))
O
Is there an specific function
to eliminate blank space within
strings?
Thanks
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: agustin.l...@ija.csic.es
http://www.ija.
Could this be a problem with
an inappropriate use of approxfun() ?
Data can be found here:
http://sites.google.com/site/eospansite/dummy/areaproblem.rda
Thanks
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
Lluis Sole Sabaris s/n
08028 Barcelona
B
5 5 142 95 112 176 161.5 133 0.42 B
6 6 89 113 112 125 128.0 195 0.42 B
...
so that I could use
xyplot(RefVal ~ XR | Band, data=datosX)
but do not get to that. Don't even know if the problem is at the melt() or
cast() level. I've also tried reshape(). Any hel
linux-gnu/4.2 -lgfortran -lm -L/usr/lib/R/lib -lR
/usr/bin/ld: cannot find -lRlapack
collect2: ld returned 1 exit status
Despite the fact that I do have Rlapack:
~$ R CMD config LAPACK_LIBS
-L/usr/lib/R/lib -lRlapack
Am I doing something wrong?
Thanks
Agus
--
Dr. Agustin Lobo
Institut de Ciencie
I wish to write "" using paste(), but
> paste("\"","Hola","\"",sep="")
[1] "\"Hola\""
>
while the same approach works with ''
> paste("\'","Hola","\'",
Whenever using lda(...,CV=TRUE), or if using lda.predict(), a
classification is returned.
Which distance is used in order to allocate a given individual to a
class, once the original data have been transformed into discriminant
components?
Thanks
Agustin Lobo
aloboa...@gmail.com
ets methods base
other attached packages:
[1] sp_0.9-99
loaded via a namespace (and not attached):
[1] grid_2.15.1lattice_0.20-6 tools_2.15.1
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
Lluis Sole Sabaris s/n
08
with identify(x,y,labels), but I do not want to get the labels
permanently plotted on the graphic)
Thanks
--
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
Lluis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
e-mail
x86_64-apple-darwin10.8.0 (64-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] lubridate_1.3.0
loaded via a namespace (and not attached):
If I do:
> a
[1] "2012_10_01_14_13_32.445"
> a2
[1] "2012_10_01_14_13_32.500"
> strptime(a,format="%Y_%M_%d_%H_%M_%S")-strptime(a2,format="%Y_%M_%d_%H_%M_%S")
Time difference of 0 secs
Is there any time object in R that would deal with thousands o
Sorry for the typo, I meant thousandths, and thanks for pointing out
the %OS format, which I had overlooked
Agus
On Tue, Oct 9, 2012 at 2:51 PM, Prof Brian Ripley wrote:
> On 09/10/2012 09:54, Agustin Lobo wrote:
>>
>> If I do:
>>>
>>> a
>>
>> [1]
d not attached):
[1] colorspace_1.2-4 dichromat_2.0-0digest_0.6.4
grid_3.0.2 gtable_0.1.2 labeling_0.2
[7] lattice_0.20-24munsell_0.4.2 proto_0.3-10
RColorBrewer_1.0-5 reshape2_1.2.2 scales_0.2.3
[13] stringr_0.6.2 tools_3.0.2
--
Agustin Lobo
aloboa...@gmail.co
I've found a problem when using
categorical variables in lmp() from package lmPerm
According to help(lmp): "This function will behave identically to lm()
if the following parameters are set: perm="", seq=TRUE,
center=FALSE.")
But not in the case of including categorical variables:
require(lmPerm)
I have a problem with duplicated date_time stamps that I do not see as
duplicated.
I read a file with observations taken every 30 minutes:
> aur2009=read.csv(paste(datadir,"AUR_ECPP_2009.csv",sep="/"),sep=";",stringsAsFactors=F)
> aur2009[1:3,1:5]
Date.Time E_filled E_filled_flag LE_filled
rptime’ does not validate such times (it does not assume a
> specific timezone), but conversion by ‘as.POSIXct’) will do so.
>
>
>
> On Wed, May 18, 2011 at 3:53 PM, Agustin Lobo
> wrote:
>>
>> I have a problem with duplicated date_time stamps that I do not
Hi!
I've tried to install playwith on Windows XP but,
> install.packages("cairoDevice")
trying URL
'http://cran.at.r-project.org/bin/windows/contrib/2.13/cairoDevice_2.15.zip'
Content type 'application/zip' length 52371 bytes (51 Kb)
opened URL
downloaded 51 Kb
package 'cairoDevice' successfully
which
the inverse function of as.POSIXct()?
I've tried as.numeric and unclass() using a origin= argument, but this does not
work.
Thanks
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
F
it be possible that future versions provide the same labeling?
rect.hclust() is useful to display the classes over the dendrogram but
the output from cutree() is often more handy.
Thanks
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
Lluis Sole Sabaris S/N
08028
; a <- chron("20100506",format="ymd")
> a
[1] 100506
> years(a)
[1] 2010
Levels: 2010
but instead have to cast to Date to apply julian:
> julian(a)
Error in names(d) : 'd' is missing
> julian(as.Date(a))
[1] 14735
attr(,"origin")
[1] "1970-01-
rmutation of the characters "d", "m",
or "y" delimited by a separator (possibly null), e.g., "m/d/y",
"d-m-y", "ymd", are all valid; "
Thanks!
Agus
On Fri, Jun 8, 2012 at 12:54 PM, Gabor Grothendieck
wrote:
> On Fri, Jun 8, 20
Hi!
I'm displaying a contingency table with heatmap():
> svm.predPix.tabla
svm.predPix CC DD LL NN NN2
CC 22 0 3 8 3
DD 0 27 0 1 0
LL 1 1 90 3 7
NN 2 0 1 11 4
NN2 0 0 5 1 20
> heatmap(svm.predPix.tabla[5:1,], Rowv=NA, Colv=NA,c
Error in plot.window(xlim = xlim, ylim = ylim, log = log, yaxs =
pars$yaxs) :
need finite 'ylim' values
Calls: plot ... boxplot -> boxplot.default -> do.call -> bxp -> plot.window
why do I get the warning in window()? My command follows what is
mentioned in the manual:
dow(...) : invalid 'ylim' value
Any help appreciated,
Thanks
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
Lluis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
e-mail agustin.l...@ictja.csic.es
https://sites.google.com
Thanks for your interest.
I've put the dataframe alyL32007 in
http://dl.dropbox.com/u/3180464/alyL32007.rda
ready to be used with load()
Agus
On Mon, May 28, 2012 at 6:11 PM, Gabor Grothendieck
wrote:
> On Mon, May 28, 2012 at 5:35 AM, Agustin Lobo
> wrote:
>> I'm do
) :
Incompatible methods ("Ops.POSIXt", "Ops.Date") for "<="
> str(delme)
‘zoo’ series (without observations)
> delme <- window(alyL32007z, start = c(2007, 3), end = c(2007, 12))
> str(delme)
‘zoo’ series (without observations)
Is this a problem with t
Month" "Day" "Hour" "doy2" ...
Index: POSIXct[1:12001], format: "2007-03-04 00:00:00" "2007-03-04
00:30:00" "2007-03-04 01:00:00" "2007-03-04 01:30:00" "2007-03-04
02:00:00" "2007-03-04 02:30:00"
="ph"
variable=="temperature"| variable== "Ba",mean,na.rm=T)
Is there a more compact form? something like select=c(ph, temperature, Ba)
Thanks
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel.
()
Any ideas on how to write an script (could be done
in R itself?) that would automatically modify
the current script into a "htmlize-friendly" script?
One problem is that while R accepts png()
providing a default name for the png file,
htmlize() fails unless a name is provided.
Than
: Factor w/ 4 levels "0","1","2","3": 4 1 4 4 3 NA 3
3 3 4 ...
> table(b2$p_desc_no3)
0 1 2
42 44 66
and
> levels(b2$p_desc_no3)[levels(b2$position_geo)=="3"] <- NA
which does not result into error but leaves b2$p_desc_no
aux <- origen
origen <- imagen
imagen <- aux
}
m <- match(v, origen, 0)
#print("match finished")
if(resto==0) v <- v*0
v[m > 0] <- imagen[m]
v
}
Agus
--
Dr. Agustin Lobo
Institut de Ciencies
- Phil Spector
Statistical Computing Facility
Department of Statistics
UC Berkeley
[EMAIL PROTECTED]
On Fri, 27 Jun 2008, Agustin Lobo wrote:
Hi!
Given a vector (or a factor within a d
umn of mydata would be the class.
But I'm looking for a nicer code and a nicer graphic
using the approach of the lattice package.
Thanks!
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fa
In the code of my previous message
the barplot should be:
barplot(apply(sel[,-1],2,mean))
instead of
barplot(sel)
Sorry for the confusion.
Agus
Mensaje original
Asunto: barplot as Trellis graphic
Fecha: Wed, 26 Mar 2008 22:24:04 +0100
De: Agustin Lobo <[EMAIL PROTEC
00),rep("y",100)),f=factor(c(datos[,3],datos[,3])))
and then
require(lattice)
barchart(V1~VAR|f,data=datos2)
I get horizontal lines in the bars that I do not understand, though.
Agus
Deepayan Sarkar escribió:
> On 3/26/08, Agustin Lobo <[EMAIL PROTECTED]> wrote:
>> De
guide to Trellis graphics?
Agus
Charilaos Skiadas escribió:
>
> On Mar 27, 2008, at 1:47 PM, Agustin Lobo wrote:
>
>> Thanks, it was a matter of reshaping the data matrix as I usually have
>> it, ie:
>> datos <-
>> data.frame(x=abs(round(rnorm(100,10,5))),y=abs
e care of recalculating
the factor levels automatically?
Thanks!
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: [EMAIL PROTECTED]
http://www.ija.csic.es/gt/obster
Thanks, problem solved thanks to your hints, but, within subset(),
drop=T still keeps
the unused levels (I'm using 2.6.1 on win).
Agus
Stefan Grosse escribió:
> On Tuesday 01 April 2008 09:43:00 am Agustin Lobo wrote:
> AL> I'm doing:
> AL> etni <- subset(etni, NAMECO
hanks!
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: [EMAIL PROTECTED]
http://www.ija.csic.es/gt/obster
__
R-help@r-project.o
e emf as vector option, and emf files are poorly imported
by OO.
Agus
Stefan Grosse escribió:
> On Wed, 16 Apr 2008 13:20:35 +0200 Agustin Lobo wrote:
> AL>I have problems with emf files because OpenOffice does
> AL>a poor job at importing figures in this format. On the other
> AL&g
;> I've seen so many problems with SVG renderers (especially related to
>> fonts) that I would prefer not to rely on it.
>>
>> If this is Windows, the preferred route seems to be to use EPS with a
>> preview (which GSView can add for you).
>>
>> On We
Stefan,
Stefan Grosse escribió:
> On Wed, 16 Apr 2008 18:51:30 +0200 Agustin Lobo wrote:
> AL> 2. I need the figure displayed on the screen, I'm using
> AL> Impress (the equivalent of ppt). Thus EPS is not an option,
>
> 1. So what please is wrong with png? I am using i
Thanks, actually I had thought on Lynx as an option, so
your opinion encourages me to try.
But wanted to solve the problem within the OO environment
as well.
Agus
Marc Schwartz escribió:
> See comments inline:
>
> Agustin Lobo wrote:
>> Thanks for your answers, let me summariz
Dear list:
Is there any way of getting the equivalent to what you get with bty="l"
for the right and bottom axes? and the equivalent to
bty="7" for the upper and left axes?
Thanks!
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC
I've set
layout(matrix(c(1,3,2,4),ncol=2))
and
par(mar=c(0, 4, 4, 0) + 0.0) and
par(mar=c(0, 0, 4, 3) + 0.0)
for the upper-left and upper-right plots, but
still get an space between both plots,
is there any way to completely eliminate that space?
Thanks
Agus
--
Dr. Agustin Lobo
Instit
0, 4, 3) + 0.0)
>>
>> for the upper-left and upper-right plots, but
>> still get an space between both plots,
>>
>> is there any way to completely eliminate that space?
>>
>> Thanks
>>
>> Agus
>> --
>> Dr. Agustin Lobo
>&
problems at importing emf files by OpenOffice
have been reported since so long ago that there is little hope that this
problem will
ever be solved within Open Office).
Thank you all.
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
080
another problem) would be making
the axes coincident with the rectangular plot region, that is, if ploting
plot(0:10,0:10)
get (0,0) just on the bottom down corner of the rectangle. I've tried
with several par() and could not find the way.
Thanks
Agus
--
Dr. Agustin Lobo
Institut de Cienci
ate your help.
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: [EMAIL PROTECTED]
http://www.ija.csic.es/gt/obster
__
R-help@r-
hanks
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: [EMAIL PROTECTED]
http://www.ija.csic.es/gt/obster
__
R-help@r-project.o
atlab,
and you can state diary on and
diary off to control what is being saved to
the file.
Thanks
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: [EMAIL
Windows)
user interface you will find for R.
HTH Vincent
Le lun. 19 mai à 03:24, Agustin Lobo a écrit :
Hi!
Is it possible to set a file to which both
commands and output would get automatically
saved? I've tried with sink(), but only get
the output. I mean something
like a combined
ore
same result.
On 05/19/08 09:06, Agustin Lobo wrote:
Hi!
Is it possible to start R by clicking
the .RData file in linux as in Windows?
I've tried with ubuntu hardy using
the right button and selecting R, but does
not work. Is there any way to set it up?
Marianne Promberger Wrote:
Yo
]
(801) 408-8111
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Agustin Lobo
Sent: Monday, May 19, 2008 1:24 AM
To: r-help@r-project.org
Subject: [R] Log or diary file
Hi!
Is it possible to set a file to which both commands and
output would get
m and ylim from 0 to 30, i.e.,
plot(runif(10,0,10),runif(10,10,30),xlim=c(0,30),ylim=c(0,30))
but this is a waste of space.
Thanks
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 93411001
liver the information you want into a form you can incorporate elsewhere,
Richard
Agustin Lobo wrote:
Hi!
If I'm plotting 2 variables with the same units
(i.e. wing length of 2 species), how can I force that X
and Y axes have the same size for length unit? For example, if
X ranges from 0 to 1
videoplay?docid=4237353244338529080&sourceid=searchfeed
)
Thanks
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: [EMAIL PROTECTED]
htt
" package:
http://cran.r-project.org/web/packages/animation/index.html
http://animation.yihui.name/
Regards,
Yihui
On Fri, Sep 12, 2008 at 8:35 PM, Agustin Lobo <[EMAIL PROTECTED]> wrote:
Hi!
Is there any way of having a greyed ("ghosted") text
(i.e, 2006) in the background of a plot?
I
Hi!
it seems that the directions in
http://wiki.r-project.org/rwiki/doku.php?id=getting-started:installation:eeepc
refer to installing R on a eeepc 701 and actually
fail for a 901 at the first step.
Does anyone have any experience on installing R on a eeepc 901?
Thanks
--
Dr. Agustin Lobo
I'll be erasing the Xandros OS in the internal SSD.
So I'd appreciate help to install
and test R on Xandros to compare.
Thanks
Agus
Agustin,
which EEE PC 901 is it? The one for Windoze or Linux (Xandros)?
el
on 9/16/08 6:50 PM Agustin Lobo said the following:
Hi!
it seems that t
rwiki/doku.php?id=getting-started:installation:eeepc
refer to installing R on a eeepc 701 and actually
Yes, it does :)
fail for a 901 at the first step.
What exactly failed? You are actually do not need the KDE desktop.
Greetings Sigbert
--
Dr. Agustin Lobo
Institut de Ciencies de la T
(centered and possibly scaled) data matrix..."
The parenthesis implies some ambiguity, but I do interpret the sentence
as indicating that the calculation should always be done using a
centered data matrix.
Finally, all the examples in the help page use centering (or scaling,
which implie
ly appreciated, I'm really
astonished. I think it might involve an encoding problem
with the non-ascii characters, but don't get to see it.
I'm using R 2.8.1 on ubuntu 8.04 (in english; And R is in English too)
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terr
Thu, 19 Mar 2009 20:04:48 +0100
From: Agustin Lobo
Reply-To: agustin.l...@ija.csic.es
To: r-help@r-project.org
Hi!
I've done a merging between 2 dataframes using merge():
delme <-
merge(miDUNS50peqB,Bnomscodmunicipis,by.x="POBLACION",by.y="NOMMUNI",all.x=T,sort=F)
elements 6,10 and 3
Is this right? Am I understanding correctly?
Thanks
Agus
Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: [EMAIL PROTECTED]
http://www.ija.csic.es/gt/obster
___
What's the most correct way of doing the equivalent to
apply(x[,2:5],1,sum)
if x is dataframe in which the only numeric fields are
in columns 2:5 ?
(using apply returns a character vector)
Thanks
Agus
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (C
79 215 84 58 127 134 15 122 101 99 ...
$ 5-9 : num 76 180 32 56 81 106 10 112 128 96 ...
$ 10-14: num 68 145 39 46 78 81 8 92 142 107 ...
$ 15-19: num 73 170 49 52 103 77 10 116 129 129 ...
> str(data.matrix(x2))
num [1:5284, 1:4] 79 215 84 58 127 134 15 122 101 99 ...
- attr(*, "
2:34 +0100
From: Agustin Lobo <[EMAIL PROTECTED]>
Reply-To: [EMAIL PROTECTED]
To: sig-geo <[EMAIL PROTECTED]>
Hi!
I have problems at reading a shp file with readOGR(), names with
non-ascii characters (accents ',`) in the dbf table are imported as i.e.
Barber\x85 del Vall\x8as
ch 3" (h= 44.8 midpoint = 6.8, x.member = 14 )
`--leaf "Branch 4" (h= 54.7 midpoint = 7.65, x.member = 20 )
correspond to the classes 1, 2, 3... in dendclass, where
dendclas <- cutree(hc,h=70)
?
If they do not, how can I get a vector with the original individuals
for eac
left and '2' for the tree on
> the right.
>
> hope the above helps.
>
> yaomin
>
> --
> Yaomin Xu
> Statstical Computing and Bioinformatics
> The Cleveland Clinic Foundation
>
>
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Alm
;Arkansas" "Tennessee" "Georgia"
> "Colorado" "Texas"
>> labels(dend2$lower[[4]])
> [1] "Idaho" "Nebraska" "Kentucky" "Montana"
> "Ohio" "
77 matches
Mail list logo