Hi
what about do inside some function a subset of your whole data frame
fff <- function( data, rows) {
data.1 <- data[1:rows,]
get all necessary stuf on data.1
return what you want
}
You can put a dimension check if you want the function to be more robust
Regards
Petr
r-help-boun...@r-pr
I have recently discovered the "playwith" library, which is great for
creating complex lattice objects.
If you start with a simple lattice plot then modify it using
playwith, you can export the code to produce the spiffed up plot.
I noticed this function at the bottom of the xyplot documentat
Thanks guys, writing smaller pieces did the trick. Still takes about 4 hours
to do but that more manageable than over 30 hours without even writing a
line.
jholtman wrote:
>
> What do you want to export it to and in what format? You will probably
> have
> a problem trying to write that large a
Im trying to export a file from R to excel using the xlsReadWritePro function
and I keep getting the error below. I get the same error when I use the
oledatetime function. My current date format is mm/dd/, which is how I
want it in excel. Also is there anyway to export different dataset
I have a variable that identifies citizen type but some of my rows are blank. I
want to replace the blank entries with the correct value which is conditional
on another variable. In SAS, my conditional statement would be as follow:
if citizenType="" and primarylanguage="English" then citizenType
Dear kind R-experts.
with miximum likelihood method i found the following estimations for my
parameters
vector:
mu
H
A1
A2
.
.
.
A10
D1
D2
.
.
.
D10
but i want to add sum to zero constraints. sum(a1:a10)=0 and sum(d1:d10)=0
how i will do this?
_
I have a data set like this
ID=c("A","A","A","A","A","A","A","B","B","B","B","B","B","B")
s=c(1.1,2.2,1.3,1.1,3.1,4.1,4.2,1.1,2.2,1.3,1.1,3.1,4.1,4.2)
d=c(1,2,3,4,5,6,7,1,2,3,4,5,6,7)
t=c(-3,-1,0,1,2,3,4,-3,-2,-1,0,1,2,3)
mydata<-data.frame(cbind(as.character(ID),as.numeric(s),as.integer(d),as.nu
On 2009-07-01, David M Smith wrote:
> REvolution Computing has just released three new packages for R to
> CRAN (under the open-source Apache 2.0 license): foreach, iterators,
> and doMC. Together, they provide a simple, scalable parallel computing
> framework for R that lets you take advantage o
Hi all,
I am a novice user of the heatmap.2 function in the Bioconductor gplots package.
I am trying to create a heatmap diagram from a hierarchical clustering result
based on the complete linkage and the pearson correlations as the distance
metrics.
After reading the help page for heatmap.2,
Barry Rowlingson wrote:
On Wed, Jul 1, 2009 at 8:41 PM, Michael wrote:
Hi all,
How could I set a timer in R, so that at fixed interval, the R program
will invoke some other functions to run some tasks?
Use timer events in the tcltk package:
z=function(){cat("Hello you!\n");tcl("after",100
Hello,
I would like to make a simple console application using R. Is there a
good library (such as ncurses) for writing console applications in R?
Thanks,
-stephen
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
On Jul 1, 2009, at 9:15 PM, Murray Cooper wrote:
Dear All,
I am new to R and slowly learning how to use the system.
The following code is an exercise I was trying.
The intent is to generate 10 random samples of size 5 from
a vector with integers 1:10 and 2 missing values. I then want
to genera
Rolf Turner auckland.ac.nz> writes:
>
>
> On 2/07/2009, at 12:20 PM, Hsiu-Khuern Tang wrote:
>
> > Is this expected behavior?
> >
> >> z <- 1:5
> >> z[1] <<- 0
> > Error in z[1] <<- 0 : object "z" not found
> >
> > The documentation seems to suggest that z will be found in the global
> > envir
Dear Murray,
Here is one way: create a function that takes k sample()s from any vector
(e.g., x) and then calculates the number of NA values in it. Then replicate
the procedure as many times as you want.
# The function
foo <- function(x, k = 5){
xsample <- sample(x, k)
Dear All,
I am new to R and slowly learning how to use the system.
The following code is an exercise I was trying.
The intent is to generate 10 random samples of size 5 from
a vector with integers 1:10 and 2 missing values. I then want
to generate a matrix, for each sample which shows the freque
An Econometrician may help you with more theoretical insights, but you could
do Monte-Carlo simulations of data and analyze the effects you are
interested in.
Daniel
-
cuncta stricte discussurus
-
-Ursprüngliche Nachricht-
Von: r-help-boun.
On 2/07/2009, at 12:20 PM, Hsiu-Khuern Tang wrote:
Is this expected behavior?
z <- 1:5
z[1] <<- 0
Error in z[1] <<- 0 : object "z" not found
The documentation seems to suggest that z will be found in the global
environment and modified accordingly.
I would agree that the documentation wou
R Colleagues:
I¹m moving toward building my own package, and it occurs to me that it might
be useful to have some method of listing or better, graphically displaying,
which functions call other functions within the package. In other words,
I¹m seeking some means of seeing how the functions relate
> More generally, you can always write a loop. They aren't necesssrily
fast
> or elegant, but they're pretty general. For example, to calculate
the max
> of the previous 50 observations (or fewer near the start of a
vector), you
> could do
>
> x <- ... some vector ...
>
> result <- numeric(le
Dear R users,
Thanks in advance.
I am Deb, Statistician at NSW Department of Commerce, Sydney.
I am using R 2.9.1 on Windows XP.
This has reference to the package âMCMCpackâ. My objective is to
select a better model using various alternatives. I have provided here
an example code from
Is this expected behavior?
> z <- 1:5
> z[1] <<- 0
Error in z[1] <<- 0 : object "z" not found
The documentation seems to suggest that z will be found in the global
environment and modified accordingly.
Best,
Hsiu-Khuern.
__
R-help@r-project.org mailin
Hi,
I just came across the following issue regarding mixed effects models:
In a longitudinal study individuals (variable ind) are observed for some
response variable. One explanatory variable, f, entering the model as
fixed effect, is a (2-level) factor. The expression of that factor is
consta
On Wednesday 01 July 2009, roger koenker wrote:
> url:www.econ.uiuc.edu/~rogerRoger Koenker
> emailrkoen...@uiuc.eduDepartment of Economics
> vox: 217-333-4558University of Illinois
> fax: 217-244-6678Urbana, IL 61801
>
> On
The mean function works in data.frames:
lapply(loadfiles, mean)
[[1]]
Time Pressure
1. 323.
[[2]]
Time Pressure
1. 323.
On Wed, Jul 1, 2009 at 1:50 PM, plpd00 wrote:
>
> Dear all,
>
> I know it is as simple as c <- (a + b)/2 to compute the average
> (element-wise)
On 2/07/2009, at 9:53 AM, Thomas Mang wrote:
Hi,
I am looking for a function which can convert a single factor (a
vector), or a list of multiple factors, into the indicator (dummy)
matrix with 0/1 entries (or other contrasting schemes). I could
cook it
myself, but I am actually sure this th
It's because the levels of factor, don't characters:
f <- factor(c("AAA's", "B(s)", "CCC"))
levels(f)[levels(f) == "AAA's"] <- "AA"
On Wed, Jul 1, 2009 at 4:43 PM, Chris Anderson wrote:
> I have a several character variables that I need to recode, but some of
> them have special characters
Try this:
xtabs(x ~ t + id, data = X)
On Wed, Jul 1, 2009 at 2:35 PM, Young Cho wrote:
> Hi, thanks everyone for any help in advance.
>
> I found myself dealing with a tabular time-series data formatted each row
> like [ time stamp, ID, values]. I made a small examples:
>
> X = data.frame(t
Hi,
I am looking for a function which can convert a single factor (a
vector), or a list of multiple factors, into the indicator (dummy)
matrix with 0/1 entries (or other contrasting schemes). I could cook it
myself, but I am actually sure this thing already exists - but forgot
the name and pa
On Wednesday 01 July 2009, roger koenker wrote:
> It's not clear to me whether you are looking for an exploratory tool
> or something more like formal inference. For the former, it seems
> that estimating a few weighted quantiles would be quite useful. at
> least
> it is rather Tukeyesque. While
Dear R Users,
I have been trying to fit models using Umacs. However I seem to be lost in
specifying random effects. Can anyone advise me (if possible with a small
hierarchical or mixed models data set) on how one can specify the random
effects sampling algorithm or possibly point me to the example
Hello, Yuhan:
If I wanted to get something sensible today, I'd do ordinary least
squares using lm(y~x), the garchFit on the residuals. This will give
you a reasonable answer except that the confidence intervals from "lm"
will not be accurate. I'd want to do normal probability plots of
Thank you so much. I emailed the people who are in charge of maintaining
miRecords and related info. I asked them for an available data transfer
protocol like ftp or similar to avoid manually downloading haindreds of
sequences from their web site.
I got no feedback at all.
Thanks again,
Maura
Dear Young,
I am sorry, my bad, it should have been
with(X, tapply(x, list(id, t), function(y) y ))
to get exactly what you asked for.
HTH,
Jorge
On Wed, Jul 1, 2009 at 4:59 PM, Jorge Ivan Velez
wrote:
> Dear Young,
> Try this:
>
> with(X, tapply(x, list(t,id), function(y) y ))
>
> HTH,
>
>
Dear Young,
Try this:
with(X, tapply(x, list(t,id), function(y) y ))
HTH,
Jorge
On Wed, Jul 1, 2009 at 1:35 PM, Young Cho wrote:
> Hi, thanks everyone for any help in advance.
>
> I found myself dealing with a tabular time-series data formatted each row
> like [ time stamp, ID, values]. I
Is this what you want:
> X = data.frame(t=c(1,1,1,2,2,2,2,3,3,3,4,4,4,5,5),id =
+ c('a','b','c','c','b','d','e','b','a','e','a','b','d','b','c'))
> X$x = rnorm(15)
> X
t id x
1 1 a -0.6264538
2 1 b 0.1836433
3 1 c -0.8356286
4 2 c 1.5952808
5 2 b 0.3295078
6 2 d -0.82046
I have a several character variables that I need to recode, but some of them
have special characters like single quote or ( ). How tell R to ignore these
special characters so I can rename them? Below is the error message I am
getting.
social$FamilySupport[social$FamilySupport=="Mr. XXXs' exte
Can't say whether its any faster but read.zoo in the devel version of
zoo can do this using the split= argument where split=2 in the
example below says to split it into time series defined by the
second column.
Lines <- '"t" "id" "x"
1 "a" -1.71941257904109
1 "b" 1.33629503083329
1 "c" 1.613373720
This may be closer; forgot about the NAs
> cast(X.m, id ~ t, function(x) if (length(x)==0) NA else sum(x))
id 1 2 3 4 5
1 a -0.6264538 NA 0.5757814 1.5117812NA
2 b 0.1836433 0.3295078 0.7383247 0.3898432 -2.214700
3 c -0.835628
Hi all,
I'm using RScalapack library for parallelizing some heavy matrix
operations required by MCMC methods for spatio-temporal models. The
package reference manuals (dated 2005) states that the library needs
LamMPI to work but we have a Linux Cluster with OpenMPI. We have found
(http://cvs.fedor
Dear all,
I know it is as simple as c <- (a + b)/2 to compute the average
(element-wise) of two data vectors. However, I can't work out to compute the
average when you have many data vectors in a directory. I have done this:
setwd("/.../data/")
listfiles <- li
Hi, thanks everyone for any help in advance.
I found myself dealing with a tabular time-series data formatted each row
like [ time stamp, ID, values]. I made a small examples:
X = data.frame(t=c(1,1,1,2,2,2,2,3,3,3,4,4,4,5,5),id =
c('a','b','c','c','b','d','e','b','a','e','a','b','d','b','c'))
By the way, you'll probably have to reinstall some or all of your
packages (and dependencies) if you are using R64.app, probably
downgrading them in the process.
--j
Steve Ellis wrote:
Dear R-helpers,
I am running R version 2.9.1 on a Mac Quad with 32Gb of RAM running
Mac OS X version 10.5.
Steve:
Are you running R64.app? If not, grab it from here:
http://r.research.att.com/R-2.9.0.pkg
(http://r.research.att.com/ under "Leopard build") .
As far as I know (and I actually just tried it this morning), the
standard R 2.9.1 package off the CRAN website is the 32 bit version,
It's not clear to me whether you are looking for an exploratory tool
or something more like formal inference. For the former, it seems
that estimating a few weighted quantiles would be quite useful. at
least
it is rather Tukeyesque. While I'm appealing to authorities, I can't
resist recallin
For another generic approach, you might be interested in the Reduce
function,
rolling <- function( x, window=seq_along(x), f=max){
Reduce(f, x[window])
}
x= c(1:10, 2:10, 15, 1)
rolling(x)
#15
rolling(x, 1:10)
#10
rolling(x, 1:12)
#10
Of course this is only part of the solution to the
On Wed, Jul 1, 2009 at 12:54 PM, Duncan Murdoch wrote:
> On 01/07/2009 1:26 PM, Mark Knecht wrote:
>>
>> On Wed, Jul 1, 2009 at 9:39 AM, Duncan Murdoch
>> wrote:
>>>
>>> On 01/07/2009 11:49 AM, Mark Knecht wrote:
Hi,
I have a data.frame that is date ordered by row number - earliest
On Wed, Jul 1, 2009 at 8:41 PM, Michael wrote:
> Hi all,
>
> How could I set a timer in R, so that at fixed interval, the R program
> will invoke some other functions to run some tasks?
>
Use timer events in the tcltk package:
> z=function(){cat("Hello you!\n");tcl("after",1000,z)}
> tcl("after"
I use Windows. Thank you!
On Wed, Jul 1, 2009 at 12:53 PM, Eduardo Leoni wrote:
> I think you are better off writing the R script and invoke it using a
> OS specific tool. For Unix-like systems there is cron.
>
> hth,
>
> -e
>
> On Wed, Jul 1, 2009 at 3:41 PM, Michael wrote:
>> Hi all,
>>
>> How c
Just wanted to leave a note on this, after I got my new iMac (and
installed R64 from the AT&T site) -- quantreg did run, after topping out
at whopping 12GB of swap space (MacOS X, at least, should theoretically
have as much swap space as there is space on the HD -- it will
dynamically increase
James Allsopp wrote:
No, that's made no difference, sorry.
Sorry I forgot to check the print method for rcorr. If P<.0001 it prints
as 0. To print under your control print the object $P from the list
created by rcorr:
r <- rcorr(. . .)
r$P
Frank
Frank E Harrell Jr wrote:
James Allsopp
I would do this as a lattice plot. Continuing with your data:
tmp <- data.frame(sapply(data, tapply, data[1:2], mean))
tmp$time <- factor(tmp$time)
xyplot(thanaa+thalcho+thalino+ponaa+pocho+poino ~ time,
group=BMIakt, data=tmp, type="l", scales=list(relation="free"),
auto.key=list(ti
On 01/07/2009 1:26 PM, Mark Knecht wrote:
On Wed, Jul 1, 2009 at 9:39 AM, Duncan Murdoch wrote:
On 01/07/2009 11:49 AM, Mark Knecht wrote:
Hi,
I have a data.frame that is date ordered by row number - earliest
date first and most current last. I want to create a couple of new
columns that show
Thanks Roger. Your comments were very helpful. Unfortunately, each of
the 'groups' in this example are derived from the same set of data, two of
which were subsets-- so it is not that unlikely that the weighted medians
were the same in some cases.
This all leads back to an operation attempting
On Wed, Jul 01, 2009 at 01:35:39PM -0400, miguel bernal wrote:
> Emacs or X-emacs with ess (Emacs Speaks Statistics) is great on Linux and
> Mac (can be the console you saw on Mac) for syntax highlight, programming
> and debugging.
Also see
http://vgoulet.act.ulaval.ca/en/ressources/ema
Hi,
I am starting to play around with neural networks and noticed that there are
several packages on the CRAN website for neural networks (AMORE, grnnR,
neural, neuralnet, maybe more if I missed them).
Are any of these packages more well-suited for newbies to neural networks?
Are there any relati
Hi all,
How could I set a timer in R, so that at fixed interval, the R program
will invoke some other functions to run some tasks?
Thank you very much!
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
On Wed, Jul 1, 2009 at 2:10 PM, Sunil
Suchindran wrote:
> #Highlight the text below (without the header)
> # read the data in from clipboard
>
> df <- do.call(data.frame, scan("clipboard", what=list(id=0,
> date="",loctype=0 ,haptype=0)))
>
> # split the data by date, sample 1 observation from each
Hi,
(apologies for initial html posting)
I have a data frame where one column is a list of lists. I would like to
subset the data frame based on membership of the lists in that column and be
able to 'denormalise' the data frame so that a row is duplicated for each of
its list elements. Example co
REvolution Computing has just released three new packages for R to
CRAN (under the open-source Apache 2.0 license): foreach, iterators,
and doMC. Together, they provide a simple, scalable parallel computing
framework for R that lets you take advantage of your multicore or
multiprocessor workstation
Hi,
I have a data frame where one column is a list of lists. I would like to
subset the data frame based on membership of the lists in that column and be
able to 'denormalise' the data frame so that a row is duplicated for each of
its list elements. Example code follows:
# The data is read in in
#Highlight the text below (without the header)
# read the data in from clipboard
df <- do.call(data.frame, scan("clipboard", what=list(id=0,
date="",loctype=0 ,haptype=0)))
# split the data by date, sample 1 observation from each split, and rbind
sampled_df <- do.call(rbind, lapply(split(df,
df$
seeliger.c...@epamail.epa.gov wrote:
There is no IDE for R in the same way that there is for other languages --
something that supports integrated versioning, debugging and testing,
perhaps using Eclipse. Boy howdee, I hope someone knows otherwise.
There is a feature-rich R plug-in for Ecl
Dear Christopher,
Try this:
merge(x,y,all=TRUE)
HTH,
Jorge
On Wed, Jul 1, 2009 at 2:51 PM, Hane, Christopher A <
christopher.h...@ingenixconsulting.com> wrote:
>
> Hello,
>
> I'm trying to vectorize some assignment statements using match(), but
> can't seem to get it correct.
>
> I have 2 da
Hello,
I'm trying to vectorize some assignment statements using match(), but
can't seem to get it correct.
I have 2 data frames each with a key column of unique values. I want to
copy a column from one frame to another where the key values are the
same. The data frames are not the same length
> ... I saw my friend has a R Console window which has automatic syntax
> reminder when he types in the first a few letters of R command. ...
You might be thinking of JGR (Jaguar) at
http://jgr.markushelbig.org/JGR.html . This editor also prompts you with
function argument lists, including for f
2009/7/1 miguel bernal
> I think there is a package to visualize the links between
> functions in a package, but I don't know its name (if anybody knows it, I
> will love to know it).
reminds me of roxygen's callgraph (relies on graphviz), is that what you
meant?
baptiste
[[alternativ
David,
Using this mail I think I found a simple solution for something I
knew I was going to have to learn about.
Thanks,
Mark
On Wed, Jul 1, 2009 at 10:48 AM, David M
Smith wrote:
> I write about R every weekday at http://blog.revolution-computing.com
> . In case you missed them, here are som
Great roundup - thank you David.
On Wed, Jul 1, 2009 at 8:48 PM, David M Smith <
da...@revolution-computing.com> wrote:
> David
--
--
My contact information:
Tal Galili
Phone number: 972-50-3373767
FaceBook: Tal Galili
My Blogs:
http://www.r-sta
Dear all,
A very basic terrain calculated as a matrix from Spatial Points Patterns:
#interpolate using the akima package
library(akima)
terrain=interp(ppoints$x,ppoints$y,ppoints$marks,xo=x0,yo=y0, linear=F)
> class(terrain)
[1] "list"
> class(terrain$x) #these are the x-coord i.e: [1...1000]
[
I write about R every weekday at http://blog.revolution-computing.com
. In case you missed them, here are some articles from last month of
particular interest to R users.
http://bit.ly/tygLz announced the release of the "foreach" and
"iterators" packages on CRAN, for simple scalable parallel progr
If you are coming to useR! next week, then you might want to check the
session on "Workbenches":
http://www.agrocampus-ouest.fr/math/useR-2009/abstracts/schedule.html
Romain
On 07/01/2009 06:58 PM, Michael wrote:
Hi all,
Could anybody point me to some latest productivity tools in R? I am
int
Here's my code
library(sound);
q1<-loadSample(" path to wav");
q2<-loadSample("path to wav");
q3<-loadSample(" path to wav");
m1<-read.table("txt",header=FALSE);
m2<-read.table("txt",header=FALSE);
m3<-read.table("txt",header=FALSE);
layout(matrix(c(1,2,3,4,5,6), 6, 1, byrow = TRUE));
par(mar = c(
Dear R-helpers,
I am running R version 2.9.1 on a Mac Quad with 32Gb of RAM running
Mac OS X version 10.5.6. With over 20Gb of RAM "free" (according to
the Activity Monitor) the following happens.
> x <- matrix(rep(0, 6600^2), ncol = 6600)
# So far so good. But I need 3 m
Emacs or X-emacs with ess (Emacs Speaks Statistics) is great on Linux and
Mac (can be the console you saw on Mac) for syntax highlight, programming
and debugging. I think there is a package to visualize the links between
functions in a package, but I don't know its name (if anybody knows it, I
will
This maps each string to one of the form yearQqtr at which point
you can sort them. Modify the mapping as necessary.
> library(gsubfn)
> dd <- c("2002", "2003H1", "2003H2", "2004", "2005Q1", "2005Q2")
> gsubfn("H.|Q.|$", list(H1 = "Q1", H2 = "Q2", Q2 = "Q2", Q3 = "Q3", Q4 = "Q4",
> "Q1"), dd)
On Wed, Jul 1, 2009 at 9:39 AM, Duncan Murdoch wrote:
> On 01/07/2009 11:49 AM, Mark Knecht wrote:
>>
>> Hi,
>> I have a data.frame that is date ordered by row number - earliest
>> date first and most current last. I want to create a couple of new
>> columns that show the max and min values from
Hi Michael,
Great topic - I hope to see others respond.
For me there are several big "time savers" with using R (on windows XP),
search them on google :
1) tinn-r, for syntax highlighting.
2) "Rexcel" package - for getting data from excel. (BTW, for excel, I also
recommend the ASAP utillities)
3)
How are you creating the plots before adding the line? This sounds like you
may be mixing graphics types (creating the plots with grid graphics (lattice or
ggplot2) then using abline from the base graphics system) or using a plot
function that plays with the graphics settings and leaves them in
It can be done without setting locales using chron:
> library(chron)
> as.Date(chron("1970-Jan-01", format = "Year-Month-Day"))
[1] "1970-01-01"
On Wed, Jul 1, 2009 at 10:09 AM, Ben Bolker wrote:
> Dear r-helpers,
>
> This is a little bit more of a Windows problem than
> an R problem, but ...
>
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
> Behalf Of Steve Jaffe
> Sent: Wednesday, July 01, 2009 9:59 AM
> To: r-help@r-project.org
> Subject: [R] sorting question
>
>
> I've asked about custom sorting before and it appears that -
Hi,
I have a multiplot of 6 rows and 1 column.I need to draw vertical lines in
each plot.However when I use abline(v=locator(1)$x) in some plots the line
only comes for half the box and it goes beyond the box in others.I suspect
this has something to do with the margins.any help?
--
Rajesh.J
I am trying to set up a Grass project and need to set up the region so
that I can view the map. I can look at a map and find the lat/lon,
but the map projection is in UTM NAD38 WGS84 and I need to set the
eastings and northings. Is there a package that will help me
calculate this in R.
thanks
--
I'm going to both UseR! (in Rennes) and DSC (in Copenhagen), and will
be blogging about the talks and other interesting things I learn here:
http://blog.revolution-computing.com/
# David Smith
--
David M Smith
Director of Community, REvolution Computing www.revolution-computing.com
Tel: +1 (20
See the by= argument.
On Wed, Jul 1, 2009 at 11:08 AM, Andriy Fetsun wrote:
> Hi,
>
> I am trying to calculate the volatility on not overlapping basis. Do you
> know functions for not overlapping calculation?
>
> It is like to take first 20 observations and apply st.dev to 20 and then
> take next
I've asked about custom sorting before and it appears that -- in terms of a
user-defined order -- it can only be done either by defining a custom class
or using various tricks with "order"
Just wondering if anyone has a clever way to order "vintages" of the form
2002, 2003H1, 2003H2, 2004, 2005Q
Hi all,
Could anybody point me to some latest productivity tools in R? I am
interested in speeding up my R programming and improving my efficiency
in terms of debugging and developing R programs.
I saw my friend has a R Console window which has automatic syntax
reminder when he types in the first
You are describing a "generalized nonlinear least-squares" estimation procedure.
This is implemented in the gnls() function in "nlme" package.
?gnls
Ravi.
Ravi Varadhan, Ph.D.
Assistant Professor,
Division of Geriatric Medici
On 01/07/2009 11:49 AM, Mark Knecht wrote:
Hi,
I have a data.frame that is date ordered by row number - earliest
date first and most current last. I want to create a couple of new
columns that show the max and min values from other columns *so far* -
not for the whole data.frame.
It seems
Hi,
I am trying to calculate the volatility on not overlapping basis. Do you
know functions for not overlapping calculation?
It is like to take first 20 observations and apply st.dev to 20 and then
take next 20 observations and calculate st. deviation.
I tried with function rollapply(), but it d
Hi all,
I'm using RScalapack library for parallelizing some heavy matrix
operations required by MCMC methods for spatio-temporal models. The
package reference manuals (dated 2005) states that the library needs
LamMPI to work but we have a Linux Cluster with OpenMPI. We have found
(http://cvs.fedor
Hi All,
I would like to do double boostrapping to estimate 95% CI coverage.
So, I can only need to estimate 95% confidence interval from each
bootstrapped sample.
Since we don't have a closed form of 95% CI, in order to get 95% CI for each
sample, we need to use bootstrapping.
For outer bootstra
Dear r-helpers,
This is a little bit more of a Windows problem than
an R problem, but ...
any idea how to query the *available* locales from
within R (or otherwise) on a Windows system? Teaching
in a Spanish-language setting and would like to do
something like
Sys.setlocale("LC_TIME","en_
Does anyone using linux KDE heard about RKWard? Is it good?
Is it better than emacs/ess? Any thought?
Ubirajara Alberton
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide ht
I tried tseries::garch() and was getting a lot of False Convergence so
tried fGarch::grachFit.
Looking a bit further I find that from fGarch::garchFit I get
..@ fit:List of 17
.. ..$ convergence: int 1
.. ..$ message: chr "singular convergence (7)"
for any and all fits. I
Hi,
I have a data.frame that is date ordered by row number - earliest
date first and most current last. I want to create a couple of new
columns that show the max and min values from other columns *so far* -
not for the whole data.frame.
It seems this sort of question is really coming from m
No, that's made no difference, sorry.
Frank E Harrell Jr wrote:
> James Allsopp wrote:
>> Hi,
>> I've just run an rcorr on some data in Spearman's mode and it's just
>> produced the following values;
>> [,1] [,2]
>> [1,] 1.00 -0.55
>> [2,] -0.55 1.00
>>
>> n= 46
>>
>>
>> P
>> [,1] [,
Dear all,
When doing nonlinear regression, we normally use nls if e are iid normal.
i learned that if the form of the variance of e is not completely known,
we can use the IRWLS (Iteratively Reweighted Least Squares )
algorithm:
for example, var e*i =*g0+g1*x*1
1. Start with *w**i = *1
2.
Hi Maura --
mau...@alice.it wrote:
> I deal with a huge amount of Biology data stored in different databases.
> The databases belongig to Bioconductor organization can be accessed through
> Bioconductor packages.
> Unluckily some useful data is stored in databases like, for instance, miRDB,
> mi
Maura,
Try the RCurl package, specifically the functions getURL and getForm.
Greg
mau...@alice.it wrote:
I deal with a huge amount of Biology data stored in different databases.
The databases belongig to Bioconductor organization can be accessed through
Bioconductor packages.
Unluckily some u
I deal with a huge amount of Biology data stored in different databases.
The databases belongig to Bioconductor organization can be accessed through
Bioconductor packages.
Unluckily some useful data is stored in databases like, for instance, miRDB,
miRecords, etc ... which offer just an
interacti
1 - 100 of 141 matches
Mail list logo