Dear R users,
is there a package in R, where one can use modified ward's method for
distance measurement??
thnx in advance
regards
eliza
[[alternative HTML version deleted]]
__
R-help@r-project.org
Hello,
I'm new to R-help, although I am a frequent user of its resources. I am
using the package REACTRAN to run some relatively simple transport modeling
(which we hope will evolve greater functionality over time). I would like to
optimize model parameters to match observed field data, and am
Dear Ozgur,
Try
sort(apply(expand.grid(a, b), 1, paste0, collapse = ))
HTH,
Jorge.-
On Fri, Sep 14, 2012 at 8:54 AM, Özgür Asar wrote:
Dear all,
I want to concatenate the elements of two vectors such as
a-c(a1,a2)
b-c(b1,b2)
and obtain
a1b1, a1b2,a2b1,a2b2
I tried the paste and
I can confirm that I have this problem occasionally as well (windows 7
laptop), I see it more often when clicking the button in the GUI that
does the same as ctrl-R.
Here is my sessionInfo:
R version 2.15.1 (2012-06-22)
Platform: x86_64-pc-mingw32/x64 (64-bit)
locale:
[1]
On Fri, Sep 14, 2012 at 1:54 PM, Özgür Asar oa...@metu.edu.tr wrote:
Dear all,
I want to concatenate the elements of two vectors such as
a-c(a1,a2)
b-c(b1,b2)
and obtain
a1b1, a1b2,a2b1,a2b2
I tried the paste and paste0 functions, but they yielded elementwise
concatenation such as
Post on R-sig-mixed-models rather than here for better responses.
-- Bert
On Fri, Sep 14, 2012 at 2:19 AM, barbara costa rbarbar...@gmail.com wrote:
Dear R users,
Does anyone knows how to run a glmm with one fixed factor and 2 random
numeric variables (indices)? Is there any way to force in
Hi,
Try this:
paste0(expand.grid(a,b)$Var1,expand.grid(a,b)$Var2)
[1] a1b1 a2b1 a1b2 a2b2
A.K.
- Original Message -
From: Özgür Asar oa...@metu.edu.tr
To: r-help@r-project.org
Cc:
Sent: Friday, September 14, 2012 8:54 AM
Subject: [R] concatenating two vectors
Dear all,
I want to
On 14-09-2012, at 14:54, Özgür Asar wrote:
Dear all,
I want to concatenate the elements of two vectors such as
a-c(a1,a2)
b-c(b1,b2)
and obtain
a1b1, a1b2,a2b1,a2b2
I tried the paste and paste0 functions, but they yielded elementwise
concatenation such as
This is also an
On Fri, Sep 14, 2012 at 4:33 PM, Jesse Robinson jrobi...@syr.edu wrote:
Hello,
I'm new to R-help, although I am a frequent user of its resources. I am
using the package REACTRAN to run some relatively simple transport modeling
(which we hope will evolve greater functionality over time). I
this function is supposed to canonicalize the language:
--8---cut here---start-8---
canonicalize.language - function (s) {
s - tolower(s)
long - nchar(s) == 5
s[long] - sub(^([a-z]{2})[-_][a-z]{2}$,\\1,s[long])
s[nchar(s) != 2 s != c] - unknown
s
}
Hi everyone, my apologies if I'm overlooking something obvious in the
documentation. I'm relatively inexperienced with the (awesome) ff package.
My goal is to use the read.table.ffdf() function to call the read.fwf()
function and pass through the colClasses and comment.char arguments. The
code
Jeff Newmiller jdnew...@dcn.davis.ca.us
on Wed, 12 Sep 2012 14:57:38 -0700 writes:
It is my normal practice to install R libraries without
root. Just use your own library directory instead of the
system library.
Thanks for the example. Makes it easy to see what you mean.
Yes, if I understand you correctly, you are right:
boxplot() (base) transforms the axes, so ?boxplot.stats, which is the
function that essentially computes the boxplot, does so on the
original data.
bwplot(lattice) transforms the data
Have you looked at the sprintf function? E.g.,
n - 3.14 * 10^(0:4)
x - state.name[35:39]
cat(sprintf(N:%8.2fSTATE:%0.8s eol\n, n, x), sep=)
N:3.14STATE:Ohio eol
N: 31.40STATE:Oklahoma eol
N: 314.00STATE:Oregon eol
N: 3140.00STATE:Pennsylv eol
N:31400.00
First thing to do is to run Rprof and see where the time is going;
here it is from my computer:
self.time self.pct total.time total.pct
tolower4.4239.46 4.42 39.46
sub3.5631.79 3.56 31.79
nchar
Sorry but what is modified Ward's method? I am guessing this has to do
with hierarchical clustering.
Ranjan
On Fri, 14 Sep 2012 16:09:00 + eliza botto
eliza_bo...@hotmail.com wrote:
Dear R users,
is there a package in R, where one can use modified ward's method for
distance
Hello,
Use ?colorRampPalette
nr - nrow(yes.matrix_11)
clrs - colorRampPalette(c(white, black))(nr)
barplot(yes.matrix_11, beside =TRUE, ylim=c(0,250), col = clrs)
[...etc...]
Hope this helps,
Rui Barradas
Em 14-09-2012 12:33, Robert Pazur escreveu:
Hi,
i created a set of graphs and want to
Dear Maitra,
yah!!! i m willing to do hclustering. i initially used simple ward but it
didnt give any satisfying result. I came know that modified ward can solve
the issue. i need to know how to do that?
thnx for your concern
regards
eliza
Date: Fri, 14 Sep 2012 12:11:57 -0500
From:
On Sep 14, 2012, at 1:46 AM, Ledile Mankga wrote:
Dear Sir /Madam
I have a problem in installing quantreg package. I recieve this message
„package Œquantreg‚ is not available (for R version 2.13.1).
Binary packages are only available from some repositories for the current
version of R.
Hi there,
I have a largish optimisation problem (10 years of daily observations).
I want to optimise between 4 and 6 parameters.
I'd like to utilise parallel computing if I can as I will have to run it
with different starting values etc.
I have a quad core PC with 16GB ram running windows 7.
On 14/09/2012 10:44 AM, Jie wrote:
Dear All,
This might be a tiny question but I do not know the reason. On my desktop,
sometimes when I use the mouse or ctrl+A select a piece of R script, and
use ctrl+R to run it,
there is no action. If I redo it, it may work.
On my laptop, it is fine. Both
On Fri, Sep 14, 2012 at 6:00 PM, Bazman76 h_a_patie...@hotmail.com wrote:
Hi there,
I have a largish optimisation problem (10 years of daily observations).
I want to optimise between 4 and 6 parameters.
I'd like to utilise parallel computing if I can as I will have to run it
with different
Dear everyone, as a newbie I fear I have made an easy mistake and am just too
blind to see it. Can anyone help me seeing the (possibly) obvious?
I have a command that creates a map of South America. In a separate
txt-file I stored data about latitude and longitude of languages, and in a
third
Sorry, I forgot to say that I use Rstudio (newest version), but it is the
same problem with Mac and Windows.
--
View this message in context:
http://r.789695.n4.nabble.com/plotting-points-on-a-map-assigning-vectors-to-values-from-dataset-tp4643182p4643183.html
Sent from the R help mailing
Hello R-helpers.
I've tried to recreate a parallel version of tapply() and table() using a
combination of the parallel functions mclapply() and pvec() and papply(), but
haven't been successful. In the end, I'm trying to get a cross tab of two
vectors. I currently (can) use tapply(..., sum)
On 14.09.2012 15:30, ramoss wrote:
It looks like they are all corrupted. I tried several other CRAN sites across
the world. How can we notify the package owner?
Take a look at the packages CRAN webpage:
http://cran.r-project.org/web/packages/xlsx/index.html
--
View this message in
Thanks for that I hadn't realised parallel could run on windows Pc's.
The code is differential evolution but its not part of a package.
Still I would like to be able to use cloud computing if possible, any
thoughts on the easiest way to achieve that using a windows based PC?
Found this blog
HI all-
I have used this .fun in S+ without a problem however, in R when I run this
code to generate multiple graphs:
trendplot-function(datafr,dataf2, abbrev=, titlestr=,
devname=s,filen=,styr=1990,endyr=2012) {
if (!is.null(dev.list())) {dev.off()}
dataf-datafr[datafr$abbrev==abbrev,]
Hello,
It never happened to me, ctrl-r always works.
But I seldom use ctrl-a and even less times click on the run icon. So it
seems to be a conflict with those.
sessionInfo()
R version 2.15.1 (2012-06-22)
Platform: i386-pc-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=Portuguese_Portugal.1252
* jim holtman wubyg...@tznvy.pbz [2012-09-14 13:10:37 -0400]:
more than half the time is in 'tolower' and 'nchar', so it is not all
'sub's problem.
aha, thanks!
This version runs a little faster since it does not need the 'tolower':
canonicalize.language - function (s) {
# s -
I have a large data.frame Z (2,424,185,944 bytes, 10,256,441 rows, 17 columns).
I want to get the result of
table(aggregate(Z$V1, FUN = length, by = list(id=Z$V2))$x)
alas, aggregate has been running for ~30 minute, RSS is 14G, VIRT is
24.3G, and no end in sight.
both V1 and V2 are characters (not
Hi,
On Fri, Sep 14, 2012 at 3:26 PM, Sam Steingold s...@gnu.org wrote:
I have a large data.frame Z (2,424,185,944 bytes, 10,256,441 rows, 17
columns).
I want to get the result of
table(aggregate(Z$V1, FUN = length, by = list(id=Z$V2))$x)
alas, aggregate has been running for ~30 minute, RSS
On Fri, Sep 14, 2012 at 7:22 PM, Bazman76 h_a_patie...@hotmail.com wrote:
Thanks for that I hadn't realised parallel could run on windows Pc's.
The code is differential evolution but its not part of a package.
Still I would like to be able to use cloud computing if possible, any
thoughts on
Hi,
On Fri, Sep 14, 2012 at 2:33 PM, agrins agri...@d.umn.edu wrote:
HI all-
I have used this .fun in S+ without a problem however, in R when I run this
code to generate multiple graphs:
trendplot-function(datafr,dataf2, abbrev=, titlestr=,
devname=s,filen=,styr=1990,endyr=2012) {
if
Using data.table will probably speed lots of things up, but also note that
aggregate(x, FUN=length, by)$x
is a slow way to compute
table(by)
.
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
-Original Message-
From: r-help-boun...@r-project.org
Well it's not pretty, but I guess you could try something like this...
best,
simon
library(mgcv)
dat - gamSim(1,n=200,dist=normal,scale=2)
b - gam(y ~ s(x1)+s(x2),data=dat)
m - 40
x1 - seq(0,1,length=m)
x2 - seq(0,1,length=m)
pd - expand.grid(x1=x1,x2=x2)
fv - predict(b,newdata=pd,se=TRUE)
Hi:
This should give you some idea of what Steve is talking about:
library(data.table)
dt - data.table(x = sample(10, 1000, replace = TRUE),
y = rnorm(1000), key = x)
dt[, .N, by = x]
system.time(dt[, .N, by = x])
...on my system, dual core 8Gb RAM running Win7
Hi all,
I have the following sql query that I am executing on a machine with single
core. I want to know how can I execute the same sqery on a maching that is
running with 4 cores. Please provide me the code.
NEW_TABLE - rhive.query(SELECT A, B, COUNT(C) FROM TABLE_A WHERE
A='01-01-2012')
Also
Hello all,
I am running a logistic model with a continuous variable as first
predictor and factor variable as the second. I would like to include
an interaction term in the following fashion:
outcome ~ continous + factor + normalized(continuous):factor
My problem is that when I do so, the
Hi,
On Fri, Sep 14, 2012 at 4:26 PM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
This should give you some idea of what Steve is talking about:
library(data.table)
dt - data.table(x = sample(10, 1000, replace = TRUE),
y = rnorm(1000), key = x)
dt[, .N, by = x]
Hello All,
I am trying to write a routine that loops over some links and parses those
links using htmlParse. The problem is that one of the links may display a pop
up window on top of that link's web page. If there is a pop up, the routine
bombs and I get an error message that the url
Hello,
Check the packages UScensus2000blkgrp and UScensus2000tract
http://stat.ethz.ch/CRAN/web/packages/UScensus2000blkgrp/index.html
http://stat.ethz.ch/CRAN/web/packages/UScensus2000tract/index.html
Hope this helps,
Rui Barradas
Em 14-09-2012 15:22, Michael Leitson escreveu:
Does anyone
One more i.e cdp.
Best Regards,
Bhupendrasinh Thakre
Sent from my iPhone
On Sep 14, 2012, at 5:30 PM, Rui Barradas ruipbarra...@sapo.pt wrote:
Hello,
Check the packages UScensus2000blkgrp and UScensus2000tract
http://stat.ethz.ch/CRAN/web/packages/UScensus2000blkgrp/index.html
You can alway convert to lower case afterwards with probably a shorter
vector. You did not indicate that you needed that conversion; it only
looked like you did it for the regular expression.
On Fri, Sep 14, 2012 at 3:13 PM, Sam Steingold s...@gnu.org wrote:
* jim holtman wubyg...@tznvy.pbz
Does anyone know if there are any special considerations with Random Forest and
correlated fields or rather derived fields?
For example if we are trying to predict who might leave our company to go work
for another company some of the variables we may look at are below (in addition
to others).
On Thu, 09-Aug-2012 at 03:40PM -0700, Kirk Fleming wrote:
| My data is 50,000 instances of about 200 predictor values, and for
| all 50,000 examples I have the actual class labels (binary). The
| data is quite unbalanced with about 10% or less of the examples
| having a positive outcome and the
I addition to Michael's suggestions, you can also check out this
tutorial which shows how to use lapply into EC2.
http://www.rinfinance.com/agenda/2012/workshop/WhitArmstrong.pdf
Unfortunately, rzmq is not available on windows, so this may not be
the best solution for your setup.
-Whit
On
On Thu, Sep 13, 2012 at 05:02:54PM -0400, li li wrote:
Dear all,
Y, X are bivariate normal with all the parameters known.
I would like to generate numbers from the distribution Y | X c
where c is a constant.
Does there exist an R function generating
random numbers from such
Dear community,
I've check it while working, but just to reassure myself. Let's say we have
2 models:
model1 - lm(vdep ~ log(v1) + v2 + v3 + I(v4^2) , data = mydata)
model2 - lm(vdep ~ log(v1) + v2 + v3 + v4^2, data = mydata)
So in model1 you really square v4; and in model2, v4*^2
On 13.09.2012 23:47, R. Michael Weylandt wrote:
On Thu, Sep 13, 2012 at 6:44 PM, Ray Griner rgri...@sdac.harvard.edu wrote:
Dear R-help list -
Here's my problem. I have some programs I want to share with other users in the
department. I currently tell users to use
On 14.09.2012 09:41, agent dunham wrote:
Dear community,
I've check it while working, but just to reassure myself. Let's say we have
2 models:
model1 - lm(vdep ~ log(v1) + v2 + v3 + I(v4^2) , data = mydata)
model2 - lm(vdep ~ log(v1) + v2 + v3 + v4^2, data = mydata)
So in model1 you
Thank you so much,
Yes I need to look up how to make a reproducible example,
William I will try your advice, I believe this will be my salvation here, once
I get my computer.
Pancho
On Sep 13, 2012, at 23:59, William Dunlap wdun...@tibco.com wrote:
-Original Message-
From:
Hi guys,
How could I specify minimum and maximum x-axis values in logi.hist.plot?
My code is something like this:
plot(mydata$Temperature,mydata$Mortality,
,xlab=Temperature,ylab=Probability of mortality)
curve(predict(temp.glm,data.frame(Temperature=x),type=resp),add=TRUE,
col=red)
Hi,
I have two data-frames which I want to match by a date-variable.
The problem now is, that the date-variable has a different format in this two
data-frames.
So simply matching by this variable wouldn't work.
In the one data-frame the date-variable contains strings that look like this:
Hello,
Try the following.
x - c(14.Aug.2012-16:32, 3.Sep.2012-16:50)
y - c(2012.08.14, 2012.09.3)
as.Date(x, format = %d.%b.%Y-%H:%M) == as.Date(y, format = %Y.%m.%d)
(Or use %in% instead of ==)
But note that the format returned by as.Date uses '-' as separator:
as.Date(x, format =
plot(..., xlim = c(a,b)) to limit the plot to [a, b]
Cheers,
M
On Friday, September 14, 2012, eddie smith wrote:
Hi guys,
How could I specify minimum and maximum x-axis values in logi.hist.plot?
My code is something like this:
plot(mydata$Temperature,mydata$Mortality,
Hello!
I am struggling for quite some time with proper printing of local
characters in pdf plot, via Sweave snippets in Rnw file.
When I am working directly within R, all is fine and I can get local
character properly, like:
pdf('figs/fig-relativeEntropy0.pdf', h=6, w=6, encoding='CP1250')
I'm not sure about R code, but wget (http://www.gnu.org/software/wget/)
should work. Here are some examples:
http://www.editcorp.com/Personal/Lars_Appel/wget/v1/wget_7.html
Cheers!
Tom
On Thu, Sep 13, 2012 at 1:37 PM, ucakmde maris_...@hotmail.com wrote:
Dear R-helpers,
I am trying to
On 14/09/2012 7:07 AM, Petar Milin wrote:
Hello!
I am struggling for quite some time with proper printing of local
characters in pdf plot, via Sweave snippets in Rnw file.
When I am working directly within R, all is fine and I can get local
character properly, like:
Estaré ausente de la oficina desde el 14/09/2012 y no volveré hasta el
01/10/2012.
Responderé a su mensaje cuando regrese.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Well, if nothing else it is probably time to update anyway.
If that does not help, we can hope that some of the gurus have some
suggestions. It really sounded like a local network problem :(
Best of luck
John Kane
Kingston ON Canada
-Original Message-
Hi everyone,
I am trying to do a horizonplot using my own time series
data. I know that there is a horizonplot function in latticeExtra, but on
closer examination i think that the graph itself is slightly wrong (it displays
some regions as triangles and i think they should be trapezoids,
Dear All,
This might be a tiny question but I do not know the reason. On my desktop,
sometimes when I use the mouse or ctrl+A select a piece of R script, and
use ctrl+R to run it,
there is no action. If I redo it, it may work.
On my laptop, it is fine. Both are running windows 7.
Best wishes,
On Sep 14, 2012, at 12:41 AM, agent dunham wrote:
Dear community,
I've check it while working, but just to reassure myself. Let's say we have
2 models:
model1 - lm(vdep ~ log(v1) + v2 + v3 + I(v4^2) , data = mydata)
If you want to create a second degree polynomial for proper
Hello!
Many thanks, but this is not working. I put:
\usepackage[utf8]{inputenc}
and than:
'pu\u010dina (open-sea)'
Blank space is printed in pdf file, after R CMD Sweave and pdflatex.
Please, any other idea?! I am losing my mind here...
Best,
PM
This isn't really an answer to your question,
Dear Sir /Madam
I have a problem in installing quantreg package. I recieve this message
package quantreg is not available (for R version 2.13.1).
Thanks in adavance for the help.
Regards,
Ledile Mankga
-
Ledile Mankga
yes thats right. I want to see all points - so I need a line with dots
--
View this message in context:
http://r.789695.n4.nabble.com/plot-dataframe-with-inconsistently-relations-tp4642766p4643116.html
Sent from the R help mailing list archive at Nabble.com.
HI,
Try this:
vec1-c(14.Aug.2012-16:32,3.Sep.2012-16:50)
vec2-c(c(2012.08.14, 2012.09.3))
vec1DateF-as.Date(vec1,format=%d.%b.%Y-%H:%M)
vec2DateF-as.Date(vec2,format=%Y.%m.%d)
ifelse(difftime(vec1DateF,vec2DateF)==0,TRUE,FALSE)
#[1] TRUE TRUE
A.K.
- Original Message -
From: Martin
Dear R users,
Does anyone knows how to run a glmm with one fixed factor and 2 random
numeric variables (indices)? Is there any way to force in the model a
separate interaction of those random variables with the fixed one?
I hope you can help me.
#eg.
Reserve - rep(c(In,Out), 100)
fReserve -
Dear all,
I want to concatenate the elements of two vectors such as
a-c(a1,a2)
b-c(b1,b2)
and obtain
a1b1, a1b2,a2b1,a2b2
I tried the paste and paste0 functions, but they yielded elementwise
concatenation such as
a1b1,a2b2
I am wondering that is there an efficient way of doing my wish, for
Given my reproducible example
test-structure(list(site = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L,
3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L), .Label = c(A,
Does anyone know the appropriate R code to obtain census tracts (tract
and block level) from addresses by using Zack Almquist's USCensus2000
package?
Thank you,
--
Michael Leitson
michael.leit...@gmail.com
__
R-help@r-project.org mailing list
Hi,
i created a set of graphs and want to find out how to swap the bar colours
(from white to black) and set it to be default for all graph that i will
produce (f.e. for different number of columns). some sugggestions?
thanks in advance.
here is an example:
##
a
On Fri, Sep 14, 2012 at 9:20 AM, Geophagus f...@retposto.net wrote:
yes thats right. I want to see all points - so I need a line with dots
+ geom_point() should work, but I'm not at my computer to test.
Also, as I said before: please quote context
Cheers,
Michael
It looks like they are all corrupted. I tried several other CRAN sites across
the world. How can we notify the package owner?
--
View this message in context:
http://r.789695.n4.nabble.com/Cannot-install-package-xlsx-tp4643054p4643142.html
Sent from the R help mailing list archive at
Hi!
I have a question concerning the rcdk package:
I generated a test sdf-file with 3 molecules in it and tried to perform
clustering by fingerprints and plot the results in a dendogramm.
This is what I did:
mols - load.molecules (molecules.sdf)
fp.list - lapply (mols, get.fingerprint,
76 matches
Mail list logo