Hi Anna,
>> How can I change the barplot so that the left hand axis scales from 0 to
>> 15 and the right hand
>> axis from 0 to 5?
Try this:
par(mfrow=c(1,1), mai=c(1.0,1.0,1.0,1.0))
Plot1<-barplot(rbind(Y1,Y2), beside=T, axes=T, names.arg=c("a","b"),
ylim=c(0,15), xlim=c(1,9), space=c(0,1), c
On 17/11/2010, at 2:37 PM, Brant Inman wrote:
> R-helpers,
>
> I have had difficulty installing the "pcvsuite" package on R version 2.12.0
> (2010-10-15). The pcvsuite package is not available on CRAN, but is located
> for download at the following website at the University of Washington:
>
try to quote the special LaTeX symbols, eg
## add more here if needed
sani <- function(x) gsub("_", "_", x)
so your example:
\documentclass[a4paper]{article}
<>=
#mytitlevar <- "Stuff" # case 1, everything is find
mytitlevar <- "Stuff_first" # case 2, f is turned into sub-text
sani <- functi
Hi
I got an table that contains this colum and i wonder if there is any nice
way to extract all rows that contains all dose that start with A02, so all
the rows that have A020 and A021 etc etc...
| A010 |
| A010 |
| A010 |
| A020 |
| A020 |
| A020 |
| A020 |
| A020 |
| A021 |
| A021 |
Thx for
dear Nick,
I do not know your data, however it seems to me that the pattern of relative
abundance
of salamanders should not exhibit a sudden change, but rather a gradual change.
If this is the case, have a look to the segmented package and references
therin. In
particular have a look to the rele
On 17-Nov-10 00:02:39, José Fernando Zea Castro wrote:
> Hello.
> First, I'm thankful about your wonderful project.
>
> However, I have serious worries about the reliability of R.
> I found the next bug which I consider important because in
> my job everytime We work with datanames like next. Plea
Hi everyone,
I wonder, if anyone can help me with this annoying issue, although it is
related to Tinn-R and Windows Vista...
I am running R 2.11.1 and Tinn-R 2.3.5.2 together. I saved 2 R-hotkeys:
ALT+A for sending the whole content of one file and ALT+S for a selection
from one file from Tinn-R to
Hi,
Someone recently told me that R is not working well with other programs in
a clustered environment. I believe we are using gensoft and working in a
redhat environment. What was also mentioned was that there are problems when
running multiple R instance at the same time.
Personally I haven't
Hi,
I try to make a xyplot withou colors. I try:
pdf(file="test.pdf")
trellis.device(color=F)
Depth <- equal.count(quakes$depth, number=8, overlap=.1)
xyplot(lat ~ long | Depth, data = quakes)
dev.off()
the graphic is showed without colors in a X device but it dont work
saving directly on pdf d
Hi,
A reproducible example would have been nice...
Here is one:
vect <- c(rep("| A010 |",3), rep("| A020 |",5), rep("| A021 |",2))
df <- data.frame(a=vect, b=rnorm(10), d=rnorm(10))
df[grep(pattern="A02", df$a), ]
HTH,
Ivan
Le 11/17/2010 10:16, Joel a écrit :
Hi
I got an table that contains
Thx alot works perfect!
//Joel
--
View this message in context:
http://r.789695.n4.nabble.com/Extract-elements-that-starts-with-tp3046332p3046420.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https
On 11/17/2010 03:06 PM, anna.richa...@csiro.au wrote:
Hi,
I hope this is a simple question. I am having trouble changing the scale of a
secondary y-axis on a barplot. When I run the code below the limits set for the
first axis are always applied to the second axis as well. I am using the lates
Hello everyone.
Consider the following number below:
0.9333 0.9333 0.06667
Is it possible to print them with less digits like 0.93 0.93 0.06 (but without
losing the R's accuracy).
I would like to thank you in advance
Best Regards
Alex
[[altern
Try this:
print(c(0.9333, 0.9333, 0.06667), 1)
On Wed, Nov 17, 2010 at 8:55 AM, Alaios wrote:
> Hello everyone.
> Consider the following number below:
>
> 0.9333 0.9333 0.06667
>
> Is it possible to print them with less digits li
Hi,
I am using the clusters function in the evd package and was wondering if
anyone knew how to set the row names as date/time objects. My data has 1
column of date and time and another of wave height (H).
My two options are:
1 - When in import the data using read.table make sure that the row.n
Thanks a lot.
Worked :)
--- On Wed, 11/17/10, Henrique Dallazuanna wrote:
From: Henrique Dallazuanna
Subject: Re: [R] reduce print accuracy
To: "Alaios"
Cc: "Rhelp"
Date: Wednesday, November 17, 2010, 10:58 AM
Try this:
print(c(0.9333, 0.9333, 0.06667), 1)
Hi!
I was wondering if there are any other functions for numerical integration,
besides 'integrate' from the stats package, but which wouldn't require the
integrand to be vectorized. Oh, and must be capable of integrating over
(-inf,+inf).
Thanks in advance,
Eduardo Horta
[[alternative
Hello everoyne,
If you have ever used matlab you should know the variable editor. You click
over the value of a variable in the workspace and it opens like a excel sheet.
Do you know if there is something like that in R . This will make easier for me
to understand what values are stored in a mat
Take a look on ?edit
edit(iris)
On Wed, Nov 17, 2010 at 10:11 AM, Alaios wrote:
> Hello everoyne,
> If you have ever used matlab you should know the variable editor. You click
> over the value of a variable in the workspace and it opens like a excel
> sheet. Do you know if there is something l
Dear all at the R-project help list.
I have run into a problem when it comes to getting values for "total inertia"
and "R-squared" for my DCA using "decorana" in the package "vegan". I have
tried the "goodness" function, but the reply indicates that it does not work
with "decorana" class object
On Wed, Nov 17, 2010 at 1:11 PM, Alaios wrote:
> Hello everoyne,
> If you have ever used matlab you should know the variable editor. You click
> over the value of a variable in the workspace and it opens like a excel
> sheet. Do you know if there is something like that in R . This will make
> e
Hi all,
I am trying to do this type of graph (code taken from the example of xYplot
help).
dfr <- expand.grid(month=1:12, continent=c('Europe','USA'),
sex=c('female','male'))
set.seed(1)
dfr <- upData(dfr,
y=month/10 + 1*(sex=='female') + 2*(continent=='Europe') +
Dear all,
I am having a hard time to figure out a suitable test for the match
between two nominal classifications of the same set of data.
I have used hierarchical clustering with multiple methods (ward,
k-means,...) to classify my dat into a set number of classesa, and I
would like to compare
Hello is there in R any operator that give you all the data of a matrix
for example in matlab
x(2,3) returns the 2ndth row and 3rdth column
x(2,:) returns all the columns of the 2nd row.
In R now I would like to print all the
CRagent[[i]][2]
CRagent[[:]][2] doesnot work of course. Other opti
On Wed, 17 Nov 2010 04:11:17 -0800, Alaios wrote:
> Hello everoyne,
> If you have ever used matlab you should know the variable editor.
> You click over the value of a variable in the workspace and it opens
> like a excel sheet. Do you know if there is something like that in R .
> This will make
Hi Alex,
Is that what you're looking for:
> df <- data.frame(a=LETTERS[1:5], b=rnorm(5))
> df
a b
1 A -0.2401323
2 B -0.9414998
3 C 0.4289836
4 D 1.9802749
5 E -0.6993612
> df[3,2]
[1] 0.4289836
> df[3,]
a b
3 C 0.4289836
> df[,2]
[1] -0.2401323 -0.9414998 0.4289836 1.980
There are really no set ways to determine a changepoint, since a
changepoint depends completely on what you decide. Recursive partitioning
will fit a best changepoint, but it will pretty much always fit one. This
function can be found in the package rpart:
> fit <- rpart(count ~ year, control =
Thanks a lot :)
I am using Rkward but I didnot about that feature
Regards
Alex
--- On Wed, 11/17/10, Lee Hachadoorian wrote:
From: Lee Hachadoorian
Subject: Re: [R] Variable Editor
To: r-h...@stat.math.ethz.ch
Date: Wednesday, November 17, 2010, 1:38 PM
On Wed, 17 Nov 2010 04:11:17 -0800, Al
I've noticed it, but I haven't looked into it much since I rarely work
on Vista. I have found that opening R before I open Tinn-R tends to
work better than using Tinn-R to open the preferred GUI.
Benjamin
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-pro
Thanks a lot for your helpful answer.
In my case now I have implemented some "struct" with the following structure
Mystruct.Map
Mystruct.xy
If I do Mystruct[2].$xy I get correctly the xy values of the second item.
I want now to print all the $xy fields of the Mystruct[[]]
I tried
Mystruct[data.f
I have a file, each column of which is a separate year, and each row of each
column is mean precipitation for that month. Looks like this (except it goes
back to 1964).
monthX2000 X2001 X2002 X2003 X2004 X2005 X2006 X2007 X2008
X2009
11.600 1.010 4.320 2.110 0.
A list contains several matrices. Over all matrices (list elements) I'd like to
access one matrix cell:
m <- matrix(1:9, nrow=3, dimnames=list(LETTERS[1:3], letters[1:3]))
l <- list(m1=m, m2=m*2, m3=m*3)
l[[3]] # works
l[[3]][1:2, ] # works
l[[1:3]][1, 1] # does not work
How can I slice all C-c
Well, without knowing what is the structure of Mystruct.Map and
Mystruct.xy, it's really difficult to know how to help you.
And what is Mystruct by the way!?
Could you send us the output from dput(your_objects) and/or
str(your_objects)?
I think you're really confused about the structure of you
> To: carus...@gmail.com
> From: jda...@usgs.gov
> Date: Wed, 17 Nov 2010 08:45:01 -0500
> CC: r-help@r-project.org; r-help-boun...@r-project.org
> Subject: Re: [R] Population abundance, change point
>
> There are really no set ways to determine a changepoint, since a
> changepoint depends comp
try this:
m <- matrix(1:9, nrow=3, dimnames=list(LETTERS[1:3], letters[1:3]))
l <- list(m1=m, m2=m*2, m3=m*3)
lapply(l, "[", 1, 1)
# or
sapply(l, "[", 1, 1)
I hope it helps.
Best,
Dimitris
On 11/17/2010 3:01 PM, soeren.vo...@eawag.ch wrote:
A list contains several matrices. Over all matr
Indeed I have looked into various non-standard changepoint analysis
methods. I figured the OP was more interested in traditional methods since
you have to spend less time justifying your methodology. Wavelets are one
potential nontraditional method, as is Significant Zero Crossings (R
package S
one approach is the following:
# say 'Data' is your data frame
DataNew <- reshape(Data, direction = "long", varying = list(2:length(Data)))
DataNew$year <- rep(2000:2009, each = 12)
DataNew
I hope it helps.
Best,
Dimitris
On 11/17/2010 3:03 PM, Graves, Gregory wrote:
I have a file, each col
It seems that I am confusing something:
List of 50
$ :List of 2
..$ CRmap: logi [1:100, 1:100] NA NA NA NA NA NA ...
..$ xy : num [1:2] 21 11
$ :List of 2
..$ CRmap: logi [1:100, 1:100] NA NA NA NA NA NA ...
..$ xy : num [1:2] 80 68
--- On Wed, 11/17/10, Ivan Calandra wrote:
Fro
Please see below.
On Wed, 2010-11-17 at 04:41 -0500, Ted Harding wrote:
> On 17-Nov-10 00:02:39, José Fernando Zea Castro wrote:
> > Hello.
> > First, I'm thankful about your wonderful project.
> >
> > However, I have serious worries about the reliability of R.
> > I found the next bug which I co
Hello everyone.
In matlab (again) there is a fucntion find that returns you the indexes where
the condition in find was met. I want the same functionality in R i.e
find(Mydata>2) to return all the indexes where the condition is met. Do you
know something like that?
Also when I try to search in
Maybe
lapply(l,function(x){x[1,1]})
or
unlist(lapply(l,function(x){x[1,1]}))
does what you want?
Benno
Am 17.Nov.2010 um 15:01 schrieb soeren.vo...@eawag.ch:
> m <- matrix(1:9, nrow=3, dimnames=list(LETTERS[1:3], letters[1:3]))
> l <- list(m1=m, m2=m*2, m3=m*3)
> l[[
There are several statistics used to compare nominal classifications, or
_partitions_ of a data set. A partition isn't quite the same in this
context because partitioned data are not restricted to a fixed number of
classes. However, the statistics used to compare partitions should also
work for the
Sorry I am a bit late to this discussion, but I can't see if you ever
got an answer. Anyhow, on your first question:
On 31/10/10 19:14, Lorenzo Isella wrote:
Dear All,
I have some questions about probit regressions.
I saw a nice introduction at
http://bit.ly/bU9xL5
and I mainly have two ques
I am trying to download and open an on-line netcdf file.
I'm using Windows XP and R 2.11.1
Here's my script
library(ncdf)
link <-
"http://ibis.grdl.noaa.gov/SAT/SeaLevelRise/slr/slr_sla_gbl_free_all_66.nc";
dest <- "C:/temp/slr_sla_gbl_free_all_66.nc"
download.file(url=link,destfile=d
Dear Gregory,
Is there an easier, cleaner way to do this? Thanks.
There are of course several ways...
(assuming yearmonth to be a data.frame)
--- 1 ---
year <- colnames (yearmonth) [-1]
year <- gsub ("^[^[:digit:]]*([[:digit:]]*[^[:digit:]]*$)", "\\1", year)
year <- as.numeric (year)
month
Hi Gregory,
is this what you want? Ok, not the most elegant way...
# using 'melt' from the 'reshape' package
library(reshape)
Data <- data.frame(month = 1:12,
x2002 = runif(12),
x2003 = runif(12),
x2004 = runif(12),
x2005 = ru
I understand now!
But I don't know of any easy way to do it.
You would expect that something like these would work:
Mystruct[[1:2]][["xy"]] ## but error: subscript out of bounds, which
makes sense because "[[" accept only a single subscript
Mystruct[1:2][["xy"]] ## but NULL, which makes sen
Thanks Mat,
I have in the meantime identified the Rand index, but not the others. I
will also have a look at profdpm, that did not pop-up in my searches.
Indeed, the interpretation is going to be critical... Could you please
elaborate on what you mean by the bootstrap process?
Thanks a lot for
On Wed, Nov 17, 2010 at 2:54 PM, D Kelly O'Day wrote:
>
> I am trying to download and open an on-line netcdf file.
>
> I'm using Windows XP and R 2.11.1
>
> Here's my script
>
> library(ncdf)
> link <-
> "http://ibis.grdl.noaa.gov/SAT/SeaLevelRise/slr/slr_sla_gbl_free_all_66.nc";
> dest <- "C:
Dickison, Daniel carnegielearning.com> writes:
>
> The documentation for agrep says it uses the Levenshtein edit distance,
> but it seems to get this wrong in certain cases when there is a
> combination of deletions and substitutions. For example:
>
> > agrep("abcd", "abcxyz", max.distance=1)
Hi again!
The "Introduction to R" or other documentation might help you to get
started.
And you should REALLY read the posting guide and provide a reproducible
example.
In this case, some sample data corresponding to the type of objects you
have, plus the result you expect would be enough (I
Elegant or not ... This approach worked as intended, the objective being to
create individual 'years' each of which was 24 months long (not 12) such that
each 'year' reflected the antecedent rainfall condition. THANKS.
I needed to do this to help us identify years that were similar to one anoth
Hi Alex,
This ought to do what you are after. For searching, I found the most
helpful thing to be thinking about what I was really after---this
helps with generating several key words to look for, which is more
likely to turn up results. I also included some functions/packages to
help with R spe
Dear people at the R-help-list.
I have run into a problem in the package "vegan", while using the function
"vegan",
and would greatly appreciate any advice how to solve the problem.
I have put in my data R, using the "decorana" function in the package "vegan"
and received nice coordinates for th
Em 16/11/2010 04:55, Tal Galili escreveu:
Hello Cesar,
Thank you for the reply.
Another question I have is if it is possible to detect the library path of
an old R install, from the terminal of the new R install.
Cheers,
Tal,
An at large reply would be no. The old install can be in a complet
Following up again. I found on the forums for the Google Apps API this
thread that seems to be about a similar issue:
http://www.google.com/support/forum/p/apps-apis/thread?tid=1c22cb44eb5cbba6&hl=en&search_impression_id=ab161b010ecf8803%3A12c5a65ce83&search_source=related_question
It's using Java
On Wed, 17 Nov 2010, Barry Rowlingson wrote:
On Wed, Nov 17, 2010 at 2:54 PM, D Kelly O'Day wrote:
I am trying to download and open an on-line netcdf file.
I'm using Windows XP and R 2.11.1
Here's my script
library(ncdf)
link <-
"http://ibis.grdl.noaa.gov/SAT/SeaLevelRise/slr/slr_sla_gbl
hi all,
after fitting a multiple linear regression
model <- lm(y ~ a + b+ c+d)
i wanted to plot diagnostics
plot(model)
but get the error message
Error in object$coefficients : $ operator is invalid for atomic vectors.
which does not make a lot of sense, since there is no "$" - i am working with
a
On 17.11.2010 10:11, Bernd Jagla wrote:
Hi,
Someone recently told me that R is not working well with other programs in
a clustered environment.
Why not, was is the problem supposed to be?
> I believe we are using gensoft and working in a
redhat environment. What was also mentioned was tha
On Nov 17, 2010, at 7:33 AM, Martin Tomko wrote:
> Dear all,
> I am having a hard time to figure out a suitable test for the match between
> two nominal classifications of the same set of data.
> I have used hierarchical clustering with multiple methods (ward, k-means,...)
> to classify my dat i
On 17/11/2010 10:28 AM, Manderscheid Katharina wrote:
hi all,
after fitting a multiple linear regression
model<- lm(y ~ a + b+ c+d)
i wanted to plot diagnostics
plot(model)
but get the error message
Error in object$coefficients : $ operator is invalid for atomic vectors.
which does not make a lot
> rowsum(value, paste(factor1, factor2, factor3))
That is dangerous in general, and always inefficient. Imagine factor1
is c("a", "a b") and factor2 is ("b c", "c"). Use interaction with
drop = T.
Hadley
--
Assistant Professor / Dobelman Family Junior Chair
Department of Statistics / Rice Uni
On Nov 17, 2010, at 4:49 AM, Ronaldo Reis Junior wrote:
Hi,
I try to make a xyplot withou colors. I try:
pdf(file="test.pdf")
trellis.device(color=F)
Depth <- equal.count(quakes$depth, number=8, overlap=.1)
xyplot(lat ~ long | Depth, data = quakes)
dev.off()
the graphic is showed without col
Follows is the exact solution to this:
v <- NULL
#note that decreasing is FALSE so preceding year preceeds
for(i in 2:46) {
kk <- melt(yearmonth[, c(1, i, i+1)], id.vars="month", variable_name="year")
v[[i-1]] <- kk[order(kk$year, decreasing=FALSE),]
}
x <- do.call(cbind, v)
write.table(x,"
Barry & Prof Ripley:
Thanks!!
mode="wb" solve the problem.
BTW, it's he.
D Kelly O'Day
--
View this message in context:
http://r.789695.n4.nabble.com/Problem-downloading-and-opening-netcdf-file-tp3046897p3047002.html
Sent from the R help mailing list archive at Nabble.com.
___
Dear R helpers,
I have two time series, and I want to perform the cross-spectral analysis
for these time series.
I would like to know whether there is any R function to generate
cross-spectrum such as co-periodogram and quadrature periodogram.
Can anyone give me an advice for cross spectral analy
Hello
How can I get multipanel conditioning graphics with rgl as I do with lattice
"|"
For example I have three variables x, y, z, w. Where x,y,z are continuous
and w is categorical or discrete.
I want to use plot3d(x,y,z) for each value of w in a panel
something like plot3d(z~x*y|w)
cheers
Hello
How can I get multipanel conditioning graphics with rgl as I do with lattice
"|"
For example I have three variables x, y, z, w. Where x,y,z are continuous
and w is categorical or discrete.
I want to use plot3d(x,y,z) for each value of w in a panel
something like plot3d(z~x*y|w)
cheers
Another useful measure to compare partitions is the adjusted Rand
index which is implemented in the library(e1071) within the
classAgreement function.
If you have your data partitions to be compared in a matricial form
(where each column is a different partition), the syntax is
ARI<-classAgreement(
MyStruct <- list(list(CRmap = rnorm(50), xy = 1:50), list(CRmap =
rnorm(50), xy = 51:100))
I wonder if it would be possible to extend indexing to accept lists.
For example:
MyStruct[[c(1, 2)]] extracts level1[1] and level2[2], so if the vector
was a list, perhaps
MyStruct[[list(c(1, 2), 2)]] ## l
On Nov 17, 2010, at 6:44 AM, Eduardo de Oliveira Horta wrote:
Hi!
I was wondering if there are any other functions for numerical
integration,
besides 'integrate' from the stats package, but which wouldn't
require the
integrand to be vectorized. Oh, and must be capable of integrating
over
Thank you Matta for the great suggestion,
I will try the additional tests. I have just been experimenting with the
e1071 package and the adjustedRand. It works perfectly, The only
outstadning question is interpretation - is there any rule of thumbs for
the level of agreement that needs to be re
Hi, I have a plot and I would like to overlay a PNG image over it. I'm
using the rasterImage function to do this, but the problem I'm facing
is working out the coordinates of the upper right corner of the final
image in user coordinates.
That is I can place the image so the lower left is located a
On Nov 17, 2010, at 9:26 AM, Alaios wrote:
It seems that I am confusing something:
List of 50
$ :List of 2
..$ CRmap: logi [1:100, 1:100] NA NA NA NA NA NA ...
..$ xy : num [1:2] 21 11
$ :List of 2
..$ CRmap: logi [1:100, 1:100] NA NA NA NA NA NA ...
..$ xy : num [1:2] 80 68
la
When a polytomous variable (nominal scale) is used as an explanatory
variable in GLM (generalized linear models), the result object or
summary of glm function does not show the effect of this variable. In
this case, anova function can be used. My question is on the appropriate
procedure(s) for tes
Hello when my code executes I receive the message that were some warnings. I
want to catch warning messages at run time so to print some local variables and
try to understand why this warning happens.
I searched on internet and I tried withCallingHandlers(
which seems to work but as I used Rkwar
On Nov 17, 2010, at 10:37 AM, Graves, Gregory wrote:
Follows is the exact solution to this:
v <- NULL
#note that decreasing is FALSE so preceding year preceeds
for(i in 2:46) {
kk <- melt(yearmonth[, c(1, i, i+1)], id.vars="month",
variable_name="year")
v[[i-1]] <- kk[order(kk$year, decre
Bill -
An excellent point, especially if you're concerned
about efficiency:
Y = list(sqrt,sin,function(u)u/2)
Ybar0 = function(u)mean(sapply(Y,function(fun)fun(u)))
Ybar1 = function(u) rowMeans(sapply(Y, function(fun) fun(u)))
system.time(one <- Vectorize(Ybar0)(seq(0,1,length=1)))
skan gmail.com> writes:
> How can I get multipanel conditioning graphics with rgl as I do with lattice
> "|"
>
> For example I have three variables x, y, z, w. Where x,y,z are continuous
> and w is categorical or discrete.
> I want to use plot3d(x,y,z) for each value of w in a panel
> somethin
Actually, after thinking about this, what I want to do is to fit a
curve (Weibull and LogLogistic) to the first 9 data points in the
series. The numbers represent percent of development, and I don't
"trust" the 1.000 point, so I want to fit these curves to get 2
estimates of the tail, both a thick
Dear All,
I made a function which gives 3 plots in one window(I used
par(mfrow=c(1,3)) in the function).
Using that function 3 times, I want to produce 9 plots in one window.
I tried par(mfrow=c(3,1)) or par(mfrow=c(3,3)) but it didn't work.
For example,
pf <- function(p) {
par(mfrow=c(1,3))
Hi:
Try this:
pf <- function(p) {
plot(c(p:(p+10)),c(1:11))
plot(c(p:(p+10)),c(2:12))
plot(c(p:(p+10)),c(3:13))
}
par(mfrow = c(3, 3))
for(i in 1:3) pf(i)
par(mfrow = c(1, 1))
HTH,
Dennis
On Wed, Nov 17, 2010 at 8:56 AM, Soyeon Kim wrote:
> Dear All,
>
> I made a function which gives 3
Dear R-users,
I am trying to make a plot in R where x and y are plotted in a regular way, but
the x axis corresponds to another set of values.
For example I have x,y and T (all 29 values)
x<- c( -1.31846232, -1.04744756, -0.87034853, -0.72883370, -0.60618971,
-0.49501845, -0.39128988, -0.
I think this is what you want for the 'axis' command:
axis(1, at = x, labels = T, las=2)
On Wed, Nov 17, 2010 at 12:19 PM, Pelt van, Saskia (KNMI)
wrote:
> Dear R-users,
>
> I am trying to make a plot in R where x and y are plotted in a regular way,
> but the x axis corresponds to another set
Hi,
I need to produce an ordinary scatter plot and it is vital that the aspect
ratio equals 1.
I set the axis as:
plot(x, y, type="n", asp=1, ,ylim=c(-80,70),xlim=c(0,100)).
The problem is that I get some 'additional' blank plot area (basically, the
lower bound of xlim becomes quite negative
Don't fixate on avoiding loops. Bury them in a function
so you don't have to see them and then you just want something
that does the right thing quickly enough. (I'm assuming
this is not a homework/puzzle type problem where you are not
allowed to use loops). E.g., the following does the job (wi
Thanks.
I wanted to avoid loops (even inside functions) because I thought they were
much slower than vector operations, and I really need an efficient code (not
an exercise!).
However, after a few tentatives I'm coming to the conclusion that loops are
actually MORE efficient (at least in the type
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of William Dunlap
> Sent: Wednesday, November 17, 2010 9:42 AM
> To: Eduardo de Oliveira Horta
> Cc: r-help@r-project.org
Hello,
is there anyone, who works with modelling of time series of count
data. I try to model traffic accidents with inar(1)-models, but I have
a lot of problems with the R-Code. I would be very grateful, if
someone helped me.
Nazli
__
R-help@r
Hi All,
I am trying to figure out how to get the position of the knots in a pspline
used in a cox model.
my.model = coxph(Surv(agein, ageout, status) ~ pspline(x), mydata) # x being
continuous
How do I find out where the knot of the spline are? I would like to know to
figure out how many case
On Windows R 2.12.0, the following seemed to work for me:
par(pty="s") ## before plotting
plot(1:10,xlim=c(0,20),ylim=c(0,30))
I do not find that setting the asp parameter does anything useful;
indeed, the Help documentation seems to be backwards (asp=x/y not
y/x). I probably misunderstood and/
Hi all. It will be great if some one will help me to solve my home task. So,
the deal : i have .pcap file, i convert it to csv using tcpdump (tcpdump -tt
-n -r x.pcap > x.csv)
CSV file looks like that :
12890084,761659 IP 10.10.20.20.47808 > 10.10.20.255.47808: UDP, length 12
12890084,761659 IP 1
Dear R-user
I used lme to fit a linear mixed model inlcuding weights=varPower().
Additionally I wanted to use glht to calculate Tukey-Kramer multiple
comparision.
error:
> glht(modelF, linfct=mcp(Species="Tukey"))
Error in glht.matrix(model = list(modelStruct = list(reStruct =
list(SubPl
On Wed, Nov 17, 2010 at 9:42 AM, William Dunlap wrote:
> Don't fixate on avoiding loops. Bury them in a function
> so you don't have to see them and then you just want something
> that does the right thing quickly enough. (I'm assuming
> this is not a homework/puzzle type problem where you are
Eduardo:
On Wed, Nov 17, 2010 at 9:53 AM, Eduardo de Oliveira Horta
wrote:
> Thanks.
>
> I wanted to avoid loops (even inside functions) because I thought they were
> much slower than vector operations, and I really need an efficient code (not
> an exercise!).
This is a common miconception. What
I am new comer in R.There r few IDE like Tinn R,VIM etc.I mean How to
use them? Do I need to install R and then install them to use or they
can work alone? Also does one install packages on R or IDEs? Can I
call/use the package from IDEs?
regards
Parth
--
Socrates, proclaimed: "I came to know one
Thanks to everyone.
Paolo
On 16 November 2010 20:14, William Dunlap wrote:
> Have you tried filter()?
> filter(a, rep(1,7)/7)
>
> Bill Dunlap
> Spotfire, TIBCO Software
> wdunlap tibco.com
>
> > -Original Message-
> > From: r-help-boun...@r-project.org
> > [mailto:r-help-boun...@r-p
Dear R users,
Here is the coxme output I obtain on my survival dataset having 3 strains
and 2 infection status (i: infected, ni: non infected)
coxme(Surv(lay) ~ infection*strain, data=datalay, random= ~1 |block)
Cox mixed-effects model fit by maximum likelihood
Data: datalay
n= 1194
Iteration
Each IDE is going to have its own setup that will also usually depend on
what type of computer you are running. EMACS, for instance, requires a
pre-existing R install with the 'ess' package (packages are usually
installed inside R). Geany, on the other hand, requires an R install and,
in some c
Hi List,
I'm hoping to get opinions for enhancing the efficiency of the following
code designed to take a vector of probabilities (outcomes) and calculate a
union of the probability space. As part of the union calculation, combn()
must be used, which returns a matrix, and the parallelized version
1 - 100 of 166 matches
Mail list logo