Hello,
There is a task to derive a linear trend, for ex. 5 units/month for
values located at the georeferenced global latidute-longitude grid.
There are several grids corresponding to time steps, say for a month
with 6 hour interval.
Sure, one can create a loop and derive a trend value for each
Dear members,
I am performing multivariate analysis on marine benthic populations
using R. At first glance I found ca and VEGANO packages to be the
suitable for the task, but neither has incorporated Detrended Canonical
Correspondence Analysis (DCCA), which is just the method I want to apply
On Tue, 2011-12-20 at 09:31 +0100, Juan Santos wrote:
Dear members,
I am performing multivariate analysis on marine benthic populations
using R. At first glance I found ca and VEGANO packages to be the
That would be the `vegan` package
suitable for the task, but neither has incorporated
Hi
I need some guidance on how to install RServ 0.6-7 on a Redhat machine
R2.14 has been installed on this machine and when I try to install Rserve it
returned the following:
[root@soso abcde]# R CMD INSTALL Rserve_0.6-7.tar.gz
Warning: unknown option â--enable-R-shlibâ
* installing to
Thank you guys for your suggestions! Florent D - that probably is the way of
writing what I had in my head. The suggestion I took up was one of
JHoltman's which he put in an e-mail to me following his initial post but it
works well with the vectors in this question:
require(plyr)
x1.6 -
Hi all,
I am very new to R (only two days of studies). I know a little bit of
statistical learning and looking for an implementation of CART and random
forest and therefore I am now studying R.
I tested with rpart and randomForest package, they are quite good.
However, I need a classification
Dear all,
is there a function similar to extractAIC based on which I can extract the
BIC (Bayesian Information Criterion) of a coxph model?
I found some functions that provide BIC in other packages, but none of them
seems to work with coxph.
Thanks,
Michael
[[alternative HTML version
Le mardi 20 décembre 2011 à 11:11 +0100, Michael Haenlein a écrit :
Dear all,
is there a function similar to extractAIC based on which I can extract the
BIC (Bayesian Information Criterion) of a coxph model?
I found some functions that provide BIC in other packages, but none of them
seems
Albert-Jan Roskam fomcl at yahoo.com writes:
Hi,
I am using maptools to plot air quality data on a map. Each measurement point
is mapped to a postal code area.
This yields pictures with discrete borders, like so:
http://dl.dropbox.com/u/27415200/baincome.png
The problem is that the size
Dear List,
I am using constrOptim to solve the following
fr1 - function(x) {
b0 - x[1]
b1 - x[2]
((1/(1+exp(-b0+b1))+(1/(1+exp(-b0)))+(1/(1+exp(-b0-b1)/3
}
As you can see, my objective function is
((1/(1+exp(-b0+b1))+(1/(1+exp(-b0)))+(1/(1+exp(-b0-b1)/3 and I would
like to
Dear community,
I would like to try prp function, I read here:
http://www.milbo.org/rpart-plot/prp.pdf,
but I cannot find in which package is included, can anybody help?
Thanks in advance,
u...@host.com
--
View this message in context:
Dear Gavin,
Thank you very much for your fast answer and your explanations.
Yes, you are right; it is in the Vegan package and I would like to run a DCCA
to remove the arch effect that I obtain in a CCA. I am working with samples
stations where I have my species identified and environmental
Hi All,
I'm wanting to convert a ragged list of values into a structured matrix for
further analysis later on, i have a solution to this problem (below) but
i'm dealing with datasets upto 1GB in size, (i have 24GB of memory so can
load it) but it takes a LONG time to run the code on a large
Dear Gavin,
Thank you very much for your fast answer and your explanations.
Yes, you are right; it is in the Vegan package and I would like to run a DCCA
to remove the arch effect that I obtain in a CCA. I am working with samples
stations where I have my species identified and environmental
Hi Jim,
Thank you for your answer, finally I resolved the problem with this code:Â
Â
inverse=function(x){
inv-numeric(length=(length(x)))
for (i in 1:length(x)) inv[i]-x[length(x)-(i-1)]
return(inv)
}
plot(Barerror$MoyArea~Barerror$Mois,type=l)
Hello,
Is there a way to find points in SpatialPoints which lie inside a given Polygon?
--
Kind regards,
Antonio Rodriges
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
That very website states clearly that prp() is in the rpart.plot package:
http://www.milbo.org/rpart-plot/
Sarah
On Tue, Dec 20, 2011 at 4:42 AM, agent dunham crossp...@hotmail.com wrote:
Dear community,
I would like to try prp function, I read here:
http://www.milbo.org/rpart-plot/prp.pdf,
On 11-12-20 12:22 AM, Nicole Marie Ford wrote:
Hello, all.
I have come across a problem. Previously, when I recoded my DV, added the
variable to my dataset and ran the multinom, I was just fine. But I am doing
it again and I am getting this error.
R Poland$trust- trust
Error in
Michael Griffiths wrote
Dear List,
I am using constrOptim to solve the following
fr1 - function(x) {
b0 - x[1]
b1 - x[2]
((1/(1+exp(-b0+b1))+(1/(1+exp(-b0)))+(1/(1+exp(-b0-b1)/3
}
As you can see, my objective function is
Hello,
I want to create contingency table. The data which should be presented in
contingency table look like this.(Its a very long list. I am copying small
part of it)
obs1obs2exp1 exp2
3 8 725 875
0 0 5870
3 7
Michael Pearmain wrote on 12/20/2011 05:21:42 AM:
Hi All,
I'm wanting to convert a ragged list of values into a structured matrix
for
further analysis later on, i have a solution to this problem (below) but
i'm dealing with datasets upto 1GB in size, (i have 24GB of memory so
can
load
Hi all,
Anybody has ever used R to automatically search the coordinates in GE? I
searched the website, but cannot find some useful information.
My addresses are not very regular, so the methods that require the strict
address formats may be not good.
Say my address is kingston ontario
On 20.12.2011 06:47, Vikram Bahure wrote:
Dear R users,
I am getting following error while using boot.ci. I have int.inc function
with 2 values. I am generating CI for the sample estimate.
*med- function(x,i) median(x[i])*
*b1- boot(int.inc,med,2)*
*ci.out- boot.ci(b1,conf =
Milan Bouchet-Valat nalimi...@club.fr
On Tue, 2011-12-20 at 12:00 +0100, r-help-requ...@r-project.org wrote:
Dear all,
is there a function similar to extractAIC based on which I can
extract the
BIC (Bayesian Information Criterion) of a coxph model?
I found some functions that provide
khai wrote on 12/19/2011 11:26:55 PM:
Hi,
I'm very new to working with sparse matrices and would like to know how
I
can column permute a sparse matrix. Here is a small example:
M1 -
spMatrix(nrow=5,ncol=6,i=sample(5,15,replace=TRUE),j=sample(6,
Gavin Simpson gavin.simpson at ucl.ac.uk writes:
On Tue, 2011-12-20 at 09:31 +0100, Juan Santos wrote:
Dear members,
I am performing multivariate analysis on marine benthic populations
using R. At first glance I found ca and VEGANO packages to be the
That would be the `vegan`
I am trying to connect to an internal database and use the sqlQuery command to
reduce and retrieve data using the following code:
channel -odbcConnect(some_dsn, uid=, pwd=)
txt-'SELECT Date, Region, Price FROM TableXYZ WHERE Type=Domestic'
sqlQuery(channel, cat(txt,\n),errors=TRUE,)
Hi,
I just get interesting results owith non-parametric kruskalmc, and i want to
get a graph similar to those obtained for tukey.
I know how to do it for tukey (plot(TukeyHSD(model))).
Anibody knows how to get a graph for kruskalmc
Thanks
-
Mario Garrido Escudero
PhD student
Dpto. de
On 20.12.2011 14:55, Mike Pfeiff wrote:
I am trying to connect to an internal database and use the sqlQuery command to
reduce and retrieve data using the following code:
channel-odbcConnect(some_dsn, uid=, pwd=)
txt-'SELECT Date, Region, Price FROM TableXYZ WHERE Type=Domestic'
Berend Hasselman wrote
Michael Griffiths wrote
Dear List,
I am using constrOptim to solve the following
fr1 - function(x) {
b0 - x[1]
b1 - x[2]
((1/(1+exp(-b0+b1))+(1/(1+exp(-b0)))+(1/(1+exp(-b0-b1)/3
}
As you can see, my objective function is
Thanks. But if I do as you suggested...
txt-'SELECT Date, Region, Price FROM TableXYZ WHERE Type=Domestic'
returns...
SELECT TradeDate, Hub, SettlementPrice FROM PowerSettlements WHERE
Contract=\ERN\
Which because of the extra \' is improper SQL form.
Maybe I should be
Bill thanks so much. I left of the as.matrix and it worked! I really
appreciate the help.
--
View this message in context:
http://r.789695.n4.nabble.com/Random-Forest-Reading-N-A-s-I-don-t-see-them-tp4201546p4218240.html
Sent from the R help mailing list archive at Nabble.com.
On 20.12.2011 16:27, Mike Pfeiff wrote:
Thanks. But if I do as you suggested...
txt-'SELECT Date, Region, Price FROM TableXYZ WHERE Type=Domestic'
returns...
SELECT TradeDate, Hub, SettlementPrice FROM PowerSettlements WHERE
Contract=\ERN\
Which because of the extra \'
Or:
require(xlsx)
test - function(x){
+a - data.frame(A=c(1,2),B=c(10,11))
+write.xlsx(a,file=a.xlsx)
+ }
test()
list.files(patt='xlsx')
[1] a.xlsx
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 12/19/11
On Tue, Dec 20, 2011 at 8:20 AM, Jean V Adams jvad...@usgs.gov wrote:
Hi Jean,
khai wrote on 12/19/2011 11:26:55 PM:
Hi,
I'm very new to working with sparse matrices and would like to know how
I
can column permute a sparse matrix. Here is a small example:
M1 -
Try this:
maxreps=1000
my.medians=rep(NA,maxreps)
for(repnum in 1:maxreps){
my.sample=sample(my.distribution,length(my.distribution)/10) #change
sampling ratio if #desired
this.rep.median=median(my.sample)
my.medians[repnum]=this.rep.median
}
#Now, you can get CI two ways
#first, IF you are
This year JSS, at www.jstatsoft.org, published eight volumes V38-V45. Five of
them were special volumes:
V38 - Competing Risks and Multi-State Models (guest editor Putter)
V41 - Statistical Software for State Space Methods (Guest editors Commandeur,
Koopman, Ooms)
V42 - Poltical Methodology
Dear All,
The second issue of the third volume of The R Journal is now available:
http://journal.r-project.org/current.html.
Thanks to everyone involved.
Heather
--
Editor in chief
heather.tur...@r-project.org
___
r-annou...@r-project.org mailing
How do I return a character string with quotes inside string?
For example, what logic do I use if I want to return the following:
Test Score=A
I tried the following
Score-paste(Test Score=,A,sep='')
But it returned a \ inside:
Test Score=\A
Any
On 20/12/2011 2:10 PM, Mike Pfeiff wrote:
How do I return a character string with quotes inside string?
For example, what logic do I use if I want to return the following:
Test Score=A
I tried the following
Score-paste(Test Score=,A,sep='')
But it returned a \ inside:
On 20.12.2011 20:10, Mike Pfeiff wrote:
How do I return a character string with quotes inside string?
For example, what logic do I use if I want to return the following:
Test Score=A
I tried the following
Score-paste(Test Score=,A,sep='')
But it returned a \ inside:
This is my list.
obs1obs2exp1 exp2
3 8 725 875
0 0 5870
3 7 435 525
10 7 754 910
0 1 145 175
and i want result in contingency table as
obs 3
hello gurus,
i have a data frame like this
HTN HTN_FDR Dyslipidemia CAD t1d_ptype[1:25]
1Y YY T1D
2 T1D
3 Ctrl_FDR
4 T1D
5Y
Thank you Uwe. I can't get my sqlQuery to work and thought that it must be the
\ embedded in the text string, I guess not if R doesn't pass the \. Thanks
for your response.
-Original Message-
From: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de]
Sent: Tuesday, December 20, 2011
Thanks very much for responses. The 2 suggestions each do a part of what I
want. Jean's proposal permutes nonzero elements within each row but shifts
the location of the nonzero element around in the j slot of each i row. And
Douglas's idea keeps the second part but permutes the entire x slot
Duncan, Thank you for referring me to Uwe's answer.
-Original Message-
From: Duncan Murdoch [mailto:murdoch.dun...@gmail.com]
Sent: Tuesday, December 20, 2011 1:14 PM
To: Mike Pfeiff
Cc: 'R-help@r-project.org'
Subject: Re: [R] Quotes inside char string
On 20/12/2011 2:10 PM, Mike
On Dec 20, 2011, at 2:31 PM, Mike Pfeiff wrote:
Thank you Uwe. I can't get my sqlQuery to work and thought that it
must be the \ embedded in the text string, I guess not if R
doesn't pass the \. Thanks for your response.
But it did pass the ''. How does that software deal with an
the short answer... which is a guess cause you didn't provide a
reproducible example... is:
your column (i think its called t1d_ptype[1:25]) is a factor and using
factors is dangerous at best.
you can check with ?str.
see ?factor for how to convert back to strings and see if your code works.
khai wrote on 12/20/2011 01:24:10 PM:
Thanks very much for responses. The 2 suggestions each do a part of what
I
want. Jean's proposal permutes nonzero elements within each row but
shifts
the location of the nonzero element around in the j slot of each i row.
And
Douglas's idea keeps the
On Dec 20, 2011, at 12:36 PM, reena wrote:
This is my list.
obs1obs2exp1 exp2
3 8 725 875
0 0 5870
3 7 435 525
10 7 754 910
0 1 145 175
and i want
Re
your column (i think its called t1d_ptype[1:25]) is a factor and using
factors is dangerous at best.
This depends on how you want to define dangerous. If t1d_ptype ought
take values from a certain set of strings then making it a factor gives
you some safety, since it warns you when you go
Fair enough and good point. How about, dangerous when used unknowingly!
On Tue, Dec 20, 2011 at 1:01 PM, William Dunlap wdun...@tibco.com wrote:
Re
your column (i think its called t1d_ptype[1:25]) is a factor and using
factors is dangerous at best.
This depends on how you want to define
On 21/12/11 01:54, Antonio Rodriges wrote:
Hello,
Is there a way to find points in SpatialPoints which lie inside a given Polygon?
You could try:
* load package spatstat
* convert your SpatialPoints object to an object of class ppp
* convert your polygon to an object of class
I just posted this on the R-hpc list. If the r-help list is the wrong list
to use please ignore this message. Thanks!
Hello,
Hopefully I have the correct server list. I'm new'ish to R and very new to
Linux. I'm looking to do work using AWS and the rzmq package is a
dependency to some other
Ben,
You can email me directly if you have problems.
You can see the build status here:
http://cran.r-project.org/web/checks/check_results_rzmq.html
I don't use windows at all, so I'm hoping for some help from the
community for windows compatibility.
If you want to be that person, then
On Tue, Dec 20, 2011 at 2:22 PM, Whit Armstrong
armstrong.w...@gmail.com wrote:
Ben,
You can email me directly if you have problems.
You can see the build status here:
http://cran.r-project.org/web/checks/check_results_rzmq.html
I don't use windows at all, so I'm hoping for some help from
I need to measure kurtosis, skew, and maybe dip test on some
distributions I have. Currently my data is in the form of 2 vectors x
and y. Where x is 10 bins and y is the number of observations found in
that bin. It seems that the measures I want to run require the actual
observations laid out
bindata - 1:5
nobs - c(2, 3, 1, 4, 3)
rep(bindata, times=nobs)
[1] 1 1 2 2 2 3 4 4 4 4 5 5 5
for the R part, and see below:
Sarah
On Tue, Dec 20, 2011 at 5:45 PM, kbrownk kbro...@gmail.com wrote:
I need to measure kurtosis, skew, and maybe dip test on some
distributions I have. Currently
Let me put my vote in for imageJ as well. It's designed to do (gosh)
image analysis stuff, whereas R is not.
Of course you can write whatever algorithms you like in R :-), but for
specialized work it's often better to use a specialized tool.
quote
From: Andrés Aragón armandres_at_gmail.com
Hi all,
I want to write some C code that I will be able to call from R. The only
problem that I have is, I want to be able to say set breakpoints and do
some general debugging without having to do trial error with the code.
Currently I am using Rstudio, but to my knowledge this doesn't exactly
Apologies for cross-posting
We would like to announce 3 statistics courses:
1. Data exploration, linear regression, GLM GAM in R. With
introduction to R. 6 - 10 February 2012. Coimbra, Portugal.
2. Introduction to linear mixed effects modelling in R. 13 - 17 February
2012. Lisbon,
Hi,
I have a list of files in one of my working directories:
chr17.chunk1.dose.fvd
chr17.chunk1.dose.fvi
chr17.chunk1.prob.fvd
chr17.chunk1.prob.fvi
...
.
chr17.chunk10.dose.fvd
chr17.chunk10.dose.fvi
chr17.chunk10.prob.fvd
chr17.chunk10.prob.fvi
And I
Thanks, Henrik! That's given me some food for thought, and evalWithTimeout
certainly seems to be a nicer way of doing this than setTimeLimit.
I should probably explain a little more clearly: I want to stop
getGraphicsEvent() without stopping the program flow (assume we're running the
code
Hi Heiko,
I run into a similar problem recently. I have a custom kernel which have
some small negative eigenvalues, possibly due to numerical error. I first
used ksvm(K, y) to train the model, but ksvm() freeze (no response to
Ctrl+c). I thought it was due to positive definiteness of the
Hi
Let me put my vote in for imageJ as well. It's designed to do (gosh)
image analysis stuff, whereas R is not.
Of course you can write whatever algorithms you like in R :-), but for
specialized work it's often better to use a specialized tool.
I agree. But I do not know any group of
Hi
hello gurus,
i have a data frame like this
HTN HTN_FDR Dyslipidemia CAD t1d_ptype[1:25]
1Y YY T1D
2 T1D
3 Ctrl_FDR
4 T1D
66 matches
Mail list logo
