r-sig-...@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
.
--
Marcelino de la Cruz Rot
Depto. de Biología y Geología
Física y Química Inorgánica
Universidad Rey Juan Carlos
Móstoles España
__
R-help@r-project.org mailing list -- To UNSUBSC
Dear all,
I will express my opinion without knowing the details of the posts John would
like to be removed.
In the current state, people posting on this and other servers have no clear
way to go when trying to remove their posts.
It is a likely event that the number of people attempting the rem
Hi all-
New to R, and I've managed to learn a bit about Lavaan package for SEM, and
pscl package for analyzing count data.
I would like to find a way to combine these analyses as I am hoping to run a
model with a multi-trait, multi-reporter factor (2 reporters by 4 traits
apiece, with one o
standard errors for zz provided by predict.coxph, turns out to be too
small.
Thank you for any suggestions.
Omar.
On Thu, Oct 11, 2012 at 1:26 PM, Omar De la Cruz C.
wrote:
> Thank you, Dr. Therneau, that was very helpful.
>
> Best regards,
>
> Omar.
>
>
> On Mon, Oct 8
Thank you, Dr. Therneau, that was very helpful.
Best regards,
Omar.
On Mon, Oct 8, 2012 at 9:58 AM, Terry Therneau wrote:
>
>> I am interested in producing the expected number of events, in a
>> recurring events setting. I am using the Andersen-Gill model, as fit
>> by the function "coxph" in
Hello,
I am interested in producing the expected number of events, in a
recurring events setting. I am using the Andersen-Gill model, as fit
by the function "coxph" in the package "survival."
I need to produce expected numbers of events for a cohort,
cumulatively, at several fixed times. My ultim
Hello,
I'm a newbie with R, and this is my first time to post here. Apologies in
advanced if I should be posting my question on another list.
I'm working on a problem on Arellano and Bond's GMM estimator, and I've chosen
to use plm/pgmm. My issue is, I'm getting different results from pgmm()
3 4 3
Pennsylvania Utah
3 3
HTH,
MArcelino
>
>I would appreciate your help!
>
>Best regards,
>Madeleine
>
Marcelino de la Cruz Rot
D
David Firth wrote:
Dear R users:
It is planned that participants in the R User Conference, useR! 2011
(August 16-18, http://R-project.org/useR-2011) will each receive a
Better :
http://www.R-project.org/useR-2011
__
R-help@r-project.org mailin
Xin Zhang wrote:
Hi all,
I have never worked with this kind of data before, so Please help me out
with it.
I have the following data set, in a csv file, looks like the following:
Jan 27, 2010 16:01:24,000 125 - - -
Jan 27, 2010 16:06:24,000 125 - - -
Jan 27, 2010 16:11:24,000 176 - - -
Jan 2
them, you can
sort them by their scores and enumerate the number of couples out of order you can get.
Hope this help.
José-Marcio
--
---
Jose Marcio MARTINS DA CRUZ http://j-chkmail.ensmp.fr
Ecole des Mines de Paris
sammyny wrote:
I did a fourier transform on a function in time domain to get the following
functions in frequency domain (in latex):
$Y_1[\omega] = \frac{1}{1-\phi_1 e^{-jw}}$
$Y_2[\omega] = \frac{1}{1-(\phi_1 + \phi_2)e^{-jw} +\phi_1\phi_2e^{-2jw}}$
How do I find the spectrum of this functio
Thanks for the hint. Helped a lot
Gabor Grothendieck wrote:
On Sat, Feb 12, 2011 at 7:48 AM, Jose-Marcio Martins da Cruz
Functions which work with ts typically assume that a full cycle is
represented by 1 unit so if a full cycle is a week then a week must be
one unit and a day must be 1
This question is surely trivial, sorry. I'm afraid I'm misunterpreting
the information I got with the documentation, and I'm a little bit
confused. I'm just an engineer with some little skills in statistics.
Well, I have a time series - 600 days long - with some weekly
periodicity inside. So
D=1 and select d using ndiffs and the KPSS test
# for stationarity.
# ... if no model can be found - return NULL
--
---
Jose Marcio MARTINS DA CRUZ http://j-chkmail.ensmp.fr
Ecole des Mines de Paris
60, bd Saint Mi
I think ther's a bug here :
bdp wrote:
Some code I have been playing with to do this follows ...
get.best.arima<- function(x.ts, minord=c(0,0,0,0,0,0),
maxord=c(2,1,1,2,1,1))
{
# function based on 'Introductory Time Series with R'
best.aic<- 1e8 # a big number
n<- leng
ear ALL,
The last part of my thesis analysis is the cross validation. Right now I am
having difficulty using the cross validation of gstat. Below are my commands
with the tsport_ace as the variable:
nfold <- 3
part <- sample(1:nfold, 69, replace = TRUE)
sel <- (part != 1)
m.model <- x2[sel, ]
m
ce) ~ 1, x2, cutoff=100, width=9)
te<- fit.variogram(v,vgm(0.0437, "Exp", 26, 0))
y <- krige(tsport_ace~1, x2, x2_SG, model = ve.fit)
spplot(y, 1, col.regions = bpy.colors(100), sp.layout =
list("sp.lines",as(rsa2,
"SpatialLines"),no.clip = TRUE))
I'm
In any way, can you help me on either converting tileplot's trellis to
shapefile/image or (solving the problem of voronoi polygons and adding another
column to represent the polygons in tripack).
Thanks a lot in advance,
Pearl dela Cruz
___
R-help,
I seem not to get an object saved neither with ".RData" extension nor
output via dput.
Whenever I try to import the above object in another worspace I just get
nothing.
> geoFeatures <- load("geoFeatures.RData")
> geoFeatures
[1] "geoFeatures"
The "geoFeatures.RData" workspace contain
R-help,
I have a data frame:
> mydata
strata nh NhWhfh
1 I 10 26 0.048 0.385
2 II 32 84 0.154 0.381
3 III 16 42 0.077 0.381
4 IV 4 11 0.020 0.364
5 V 10 26 0.048 0.385
7 VII 64 168 0.309 0.381
8VIII 49 129 0.237 0.380
9 IX 22 58 0.107 0.37
R-help,
I have a list whose elements are data frames.
I want to change the colnames attribute in each element of this list but an
error message
comes up:
> lapply(LD_strataNew,function(x) dimnames(x)[[2]][-1]) <- as.roman(1:9)[-6]
Error in lapply(LD_strataNew, function(x) dimnames(x)[[2]][-1]) <
R-help,
I'm just trying to find the ML (maximum likelihood) estimates
of the mean and standard deviation of a set of observations:
>xx=c(2.5,3.5,4,6,6.5,7.5)
fn<-function(params,x=xx)
{
media<-params[1]
st <-params[2]
pdf=-sum(dnorm(log(xx),log(media),st,TRUE))
return(pdf)
}
optim(c(mu,stdev
R-help,
I'm trying to estimate some parameters using the Maximum Likehood method.
The model describes fish growth using a sigmoidal-type of curve:
fn_w <- function(params) {
Winf <- params[1]
k <- params[2]
t0 <- params[3]
b <- params[4]
R-help,
I 'm trying to optimize a model to data using log-likelihoods
but I encounter the following error message:
> l= c(49.4, 57.7,64.8,70.9,78.7,86.6,88.3,91.6,99,115)
> t=3:12
> fn <- function(params, l=l, t=t) {
Linf <- params[1]
k <- params[2]
R-help,
At the moment the R version installed on my machine is 2.8.1. (Ubuntu 9.04)
I wish to upgrade to R 2.9.1.
I did:
$ sudo apt-get upgrade
..but R is not upgraded although the "sources.list" file is updated with:
deb http://cran.ii.uib.no/bin/linux/ubuntu jaunty/
When I run from the term
t on this I would need also help
as there is no example on Oracle drivers (just MySql, PostgreSQLt)
Thanks in advance
----
Luis Ridao Cruz
Faroe Marine Research Institute
Nóatún 1, P.O. Box 3051
FO-110 Tórshavn
Faroe Islands
Tel : (+298
R-help,
I get the following error message when trying to connect to an Oracle database
through R (2.8.1) under Linux (Ubuntu 9.04).
> channel<-odbcConnect("magnus",uid="luisr",pwd="juanayzakarias")
Warning messages:
1: In odbcDriverConnect(st, ...) :
[RODBC] ERROR: state IM002, code 0, message
R-help,
I'm now on the process of installing RODBC on UBUNTU
but an error occured. I copy-paste the output:
> install.packages("RODBC")
Warning in install.packages("RODBC") :
argument 'lib' is missing: using
'/home/luisridaocruz/R/x86_64-pc-linux-gnu-library/2.8'
trying URL 'http://cran.ii.uib
R-help,
I'm trying to install the rgdal package under Ubuntu
but I get the following warning:
> install.packages("rgdal")
Warning in install.packages("rgdal", configure.args =
"--with-gdal-modules=/usr/local/lib") :
argument 'lib' is missing: using
'/home/luisridaocruz/R/x86_64-pc-linux-gnu-
R-help,
I have installed R under Ubuntu and
I'm very new to a Linux distribution.
To open an empty R session I just type R on the Terminal aplication.
But how can I open a saved workspace?
At present I just start R and then load ("my_workspace") but it must be possible
to do it all at oncer
wo populations. Syrjala SE. Ecology. 1996;77(1):7580.
http://dl.getdropbox.com/u/1047321/Syrjala1996.pdf
[2] https://stat.ethz.ch/pipermail/r-sig-geo/2008-February/ thread.html#3137
___
R-sig-Geo mailing list
r-sig-...@stat.math.ethz.ch
https://stat.ethz.ch/mailman/
R-help,
I'm sorry that I said " I don't have idea ".
What I mean t is that I just don't know how to set a legend
when "image" is called.
Thanks in advance
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
R-help,
I have an image plot and I wish to set a legend to it
but I don't have an idea how to do it. Below it's the code line.
image(yyy, col = rev(heat.colors(10)), axes = T)
Thanks in advance.
__
R-help@r-project.org mailing list
https://stat.ethz.c
")
dist
cluster=hclust(dist,method="ward")
cluster
plot(cluster)
rect.hclust(cluster, k=4, border="red")
I extracted 4 clusters from the data. My question is: is it possible to
produce a summary of every mean values for each variable of each of the 4
clusters?
Thanks a lot
existing functions.
Thanks!
Marcelino de la Cruz Rot
Departamento de Biología Vegetal
E.U.T.I. Agrícola
Universidad Politécnica de Madrid
28040-Madrid
Tel.: 91 336 54 35
Fax: 91 336 56 56
marcelino.delac...@upm.es
_
[[alte
R-help,
How can I specify the distance between the axis and axis annotation?
Thanks in advance
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.htm
ight be syntax errors, I just want to clearly describe my problem
Since paste("X$a>=5") gives type of character; whereas the type of
X$a>=5 should be logical.
How can I do this?
All thoughts are greatly appreciated.
Many Thanks,
cruz
__
R
try this:
temp2 <- temp2[order(temp2[,1]),] # PROBLEM IS HERE
dim(temp2) <- c(5,1)
print(temp2)
On Wed, Feb 11, 2009 at 2:23 PM, wrote:
> i'm sorry. i had an error in my previous code because i left out a letter in
> the rownames.
> while fixing that, i also found a solution. so i'm sorry fo
x27;t
even start a R session in Emacs when I have read the ESS manual
several times, but it is getting clearer to me now as I have been
writing R codes by Emacs for a couple of weeks...
It is exactly the same as we code in Windows.
Thanks,
cruz
PS: The R gurus may find this thread annoying cause it i
On Fri, Feb 6, 2009 at 8:12 PM, Harsh wrote:
> Hell R-list,
>
> At the cost of sounding far-fetched and almost incredulous, I would
> like to know if any R user is remotely considering the use of R on
> Mobile devices, and Android in particular.
>
Check out FreeRunner from Openmoko (http://wiki.o
A
> Dpto. Sistemas y Recursos Forestales
> Ctra. La Coruña, Km. 7.5
> E-28040 Madrid (Spain )
> Tel.: +34 91 347 1461
> e-mail: dmartin inia.es /// dmartinbenito gmail.com
>
> <http://dmartinbenito.googlepages.com>http://dmartinbenito.googlepages.com
>
There is an interesting article "An Introduction to Text Mining in R"
by Ingo Feinerer on R News Volume 8/2, October 2008
(http://www.r-project.org/doc/Rnews/Rnews_2008-2.pdf)
Check it out
On Fri, Feb 6, 2009 at 9:16 AM, spiketide wrote:
>
> hi everyone...
>
> i am a newbie to text mining.
(paste("x", 2, sep = "")) :
non-numeric argument to binary operator
>
All comments are appreciated.
Thanks,
cruz
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
R-help,
I'm using the "for" control flow to graph plots continuously:
par(mfrow=c(3,5), mar=c(0.6,.2,1.2,.2),yaxt="n",xaxt="n")
for(j in 1:11)
{
for(i in 1:15)
{
species <- spAldur[spAldur$ar == 1993+i & spAldur$aldur == j,]
plot(spec...@coords[,1], spec...@coords[,2], xlim=c(-10,-3.5),
abntest[abntest$abntesr > 90,]
On Sun, Jan 4, 2009 at 11:55 AM, wrote:
> Hi:
>
> I'm loading in students test scores with:
>
>> abntest <- read.table("scores.txt")
>
> if I type:
>
>>abntest
>
> I get ALL the values. I want to be able to filter it by various things such
> as:
>
>
> if( abntes
involved.
Thanks a lot,
cruz
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
7; and 'y' lengths differ
In addition: Warning messages:
1: pseudoinverse used at 0.52416 0.029047
2: neighborhood radius 3.8766
3: reciprocal condition number 4.3293e-017
4: There are other near singularities as well. 19.538
>
Any idea woul
I need help with these R simulations
___
The generation of pseudo-random numbers with uniform distribution [0.1]
can be done through a method called Congruencial Committee, which consists
the use of a recursive relationship of the type:
xi +1 = (a.xi + c) mod
Hi,
I am residing in Taipei, Taiwan. My windows has English R (2.7.2)
installed. I noticed that the CRAN mirror for Taipei is Taiwan
(Taipeh), Is Taipeh a Deutsch or English name?
Thanks,
cruz
__
R-help@r-project.org mailing list
https://stat.ethz.ch
i tried to run your code, this is how/why you got NaN:
> mu<-0.2
> sig<-0.2
> S0<-100
> j<-0.2
> dt<-1/252
> int<-0.1
> i<-0
> is.nan
function (x) .Primitive("is.nan")
> k<-rnorm(1,0,1)
> k
[1] 0.3214954
> theta<-ifelse((k<(int*dt)),1,0)
> theta
[1] 0
> m<-rnorm(1)
> m
[1] -0.4525731
> ga
is this what you expect?
> ur.df(PIBTUN$V1,type='none',lags=1)
###
# Augmented Dickey-Fuller Test Unit Root / Cointegration Test #
###
The value of the test statistic is: 6.252
Hi,
What formula is appropriate for calculating the estimate of variance
from the parametric fit? I have a linear regression model i.e. lm(a~b)
Thanks,
cruz
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
use the trick from this thread:
https://stat.ethz.ch/pipermail/r-help/2008-November/178985.html
your problem shall be fixed
On Fri, Nov 7, 2008 at 8:21 PM, mentor_ <[EMAIL PROTECTED]> wrote:
>
> Hi,
>
> imagine the following two matrix:
> Matrix A:
> a b c
> 1 2 3
> 4 5 6
> 7 8 9
>
> Matrix B:
hi friend,
this is from your previous posts on Kruskal-Wallis test:)
i came up with this one:
A5 <- read.table('kew.dat' ,header=TRUE)
plot(factor(A5$Month, levels=month.abb), A5$Rain)
is that what you want?
On Sat, Nov 8, 2008 at 7:03 AM, Swanton0822 <[EMAIL PROTECTED]> wrote:
>
> Hi.
> i hav
is this what you want?
> vector1
[1] 65 1 34 100 42 20 79 43 89 10
> vector2
[1] 34 65 47 91 48 32 23 74 92 86
>
> for (i in 1:10) {
+ for (j in 1:10) {
+ if (vector1[i] == vector2[j])
+ show(c(i,j))
+ }
+ }
[1] 1 2
[1] 3 1
On Sat, Nov 8, 2008
them give me 20.090
Thanks,
cruz
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
>
> Does that answer your question?
>
Thanks:)
I received one from Erin:
x <- NULL
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provi
n of this object
should not be fixed because it will grow...
so in MATLAB, i.e. we define x = [ ],
what about in R?
Thanks,
cruz
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://w
Hi,
I am writing this in a wrong way, can someone please correct me?
> A <- matrix()
> length(A) <- 6
> dim(A) <- c(3,2)
> colnames(A) <- c("X","Y")
> A
X Y
[1,] NA NA
[2,] NA NA
[3,] NA NA
> A$X
Error in A$X : $ oper
On Thu, Nov 6, 2008 at 4:23 PM, Peter Jepsen <[EMAIL PROTECTED]> wrote:
>
> Here is an example:
>
> id <- c(rep("a",4),rep("b",2), rep("c",5), rep("d",1))
> start <- c(c(0,6,17,20),c(0,1),c(0,5,10,11,50),c(0))
> stop <- c(c(6,12,20,30),c(1,10),c(3,10,11,30,55),c(6))
> data <- as.data.frame(cbind(id
On Thu, Nov 6, 2008 at 8:53 AM, dinesh kumar <[EMAIL PROTECTED]> wrote:
> Dear R users,
>
> I have a matrix like
>
> A X1
> B Y2
> C Z3
>
> I want to reshape this matrix into this format
>
>X Y Z
> A 1
> B 2
> C 3
>
is this what you
>
> How can I assign df$counts values to dd$counts for the corresponding
> dd$comp==df$comp
>
?factor the $comp then assign accordingly
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
On Wed, Nov 5, 2008 at 6:47 AM, Swanton0822 <[EMAIL PROTECTED]> wrote:
>
> after that i have been ask to do the following:
> 'The regression equation for EXPEND depending on ECAB is
> E(EXPEND)=119.0+1.73xECAB when there is outliners and
> E(EXPEND)=97.8+1.94xECAB when the outliners are omitted fro
r)
> Q
[,1][,2] [,3]
[1,] -0.7071068 -0.70710680
[2,] 0.000 0.000 1
[3,] -0.7071068 0.70710680
Thanks,
cruz
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Hi,
Is it possible to construct a Q from qr() that some of the rows could
be specified (to be fixed values)?
Thanks,
cruz
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R
Hi,
I have an object, Z:
> typeof(Z)
[1] "list"
It looks like:
[,1] [,2] [,3]
[1,] 111
[2,] 111
How can it label the Rows and Columns with Text? so that it looks like:
Y1 Y2 Y3
X1 111
X2 111
R-help,
Whenever I try to save a plot with "savePlot"
the file is not stored in my hard disk with the selected
format. Several formats are set and none of them
works. I just get the file name with missing extension
and it can't be open with programs like Paint and Microsoft Photo
Editor
Th only on
R-help,
I'm trying to remove some packages with the command
"remove.packages" but R 2.4.1 will come out with this:
> remove.packages(pkgs="FLAssess", lib = "C:\\Program
Files\\R\\R-2.4.1\\library\\FLAssess")
Error in .find.package(pkgs, lib) : there is no package called
'FLAssess'
Thanks in adv
R-help,
I have two data frames df1 and df2:
>head(df1)
lgdcm 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 0 0
R-help,
I have a data frame in which I compute maximum by rows
and I wish to find out the column name at which this maximun occurs.
Does anyone know how to do it?
Thanks in advance
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/list
R-help,
I'm using the code below to plot a shaded area graph.
At the same time I want to plot a second series on the y-axis (from
par(new=T) on)
but as the two series have different x-axis range (first 1994:2007 and
second 1996:2007)
the corresponding x's do not match.
How can this be sorted ou
R-help,
I'm using the code below to plot a shaded area graph.
At the same time I want to plot a second series on the y-axis (from
par(new=T) on)
but as the two series have different x-axis range (first 1994:2007 and
second 1996:2007)
the corresponding x's do not match.
How can this be sorted ou
R-help,
I have a data frame with one column containing roman numbers
The data are not sorted as : III III IV VVI VII VIII IX
XXI XII XIII XIV XV
Using data[order(data$Roman),] does not do the job.
How can this be done?
Thanks in advance.
__
R-help,
Sorry for this question (I guess it has been addressed before but
I could not find it in the archives)
but how can I see a function code
when the following comes up:
> svymean
function (x, design, na.rm = FALSE, ...)
{
.svycheck(design)
UseMethod("svymean", design)
}
Thanks i
R-help,
I have a 3-way array:
> dim(bugvinP)
[1] 13 14 3
The array looks something like this (object trimmed for readability)
, , slag = 1
ar
1994 1995 1996 1997 1998
1 NA 0.000 0.000 0.000 0.000
2 0.036 0.059 0.027 0.000 0.000
3 0.276 0.475 0.491 0.510 0.559
10 1.0
This is additional information to my previous email (pasted below)
## I opened an R session ( R 2.6.0 on Windows) and got the following
pop-up window:
## "This application has failed to start because tk84.dll was not
found. Re-installing the application may fix the problem"
## I have now upgr
R-help,
I opened an R session ( R 2.6.0 on Windows) and got the following
pop-up window:
"This application has failed to start because tk84.dll was not found.
Re-installing the application may fix the problem"
I have now upgraded to R 2.6.1 but the same fatal error comes up.
Can someone let me
R-help,
Sorry if this question has been discussed/posted before
but I can't just find it myself.
How can I print the comment character (") ?
Thanks in advance
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
R-help,
Whenever I try to save a plot with extension ".pdf" , "ps" or "eps"
I get the following error/warning message:
plot(rnorm(10))
savePlot("test",type="pdf")
Error in savePlot("test", type = "pdf") : Invalid font type
In addition: Warning messages:
1: In savePlot("test", type = "pdf") :
f
R-Help,
I have a 3 dimensional array with dimensions:
> dim(test)
[1] 241 104 12
I copy-paste an example of the object
> test[1:5,1:5,1:2]
, , ar = 1996
slagnr
puntar12 34 5
I1 60.30303 36.50610 0 0.0 0
I10 58.19512 35.48209 0 75.0 0
I11 0.
81 matches
Mail list logo