Hi:
Since Deducer was mentioned, I'll add that Ian Fellows has recently released
an 'all-in-one' installer for R, JGR and Deducer at
http://ifellows.ucsd.edu/pmwiki/pmwiki.php?n=Main.WindowsInstallation
Look for the Download Installer link (it's kinda hard to miss :)
On my system (64-bit Win 7),
Did you load the package in R? Sounds like you didn't. Try
require(nlme)
Another possibility is that that function is hidden. Just to be sure, are you
should be using it like this
summary(mod)
Not like this
summary.lme(mod)
HTH
Gavin
Sent from my HTC
- Reply message -
From: "Joann
Dear all,
Since version R-2.12.0, I've got problems building an add-on package for R
("imglib" here) for 32-bit Windows.
The following command is called:
C:\Program Files\R\R-2.12.0\bin\i386\Rcmd build --binary imglib
And here is an error I am getting in the end of processing:
** testing if ins
How do I add data to a .rdata file? In my case, I have a time series that
needs to get updated every day.
Thanks,
Jason
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLE
Eduardo de Oliveira Horta gmail.com> writes:
>
> -- Forwarded message --
> From: Eduardo de Oliveira Horta gmail.com>
> Date: Wed, Nov 17, 2010 at 3:59 PM
> Subject: Re: [R] Numerical integration
> To: David Winsemius comcast.net>
>
> It works, however is not very efficient fo
On 2010-11-17 09:26, Silvia Cecere wrote:
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
Hi all,
I am running R.app on Mac OS X 10.4. I am trying to figure
out if there is a .plist file of some sort, or a similar
parameter file, that will let me turn on and off e.g.
line-wrapping in the R.app GUI console.
The only hint I could find online is this:
https://stat.ethz.ch/pipermail
Hi,
note that to expand tilde, there is also path.expand() which should
give less "surprises" than Sys.glob().
I can confirm that this is not only Windows, but also on Linux;
> dir("~/ttt")
[1] "foo.txt"
> unlink("~/ttt", recursive=TRUE)
> dir("~/ttt")
[1] "foo.txt"
> unlink(path.expand("~/ttt")
some comparisons for interfaces to R on desktop-draft from blog post
(JSS paper is submitted for GUI issue but not vetted)
http://decisionstats.com/2010/10/05/interfaces-to-r/
RKward is number 4
have you researched out all R commander E plugins before settling on R Kward
R commander is by J F
Hi all,
Has any folk tested or been using RKward? Please shed me some light whether it
takes data from the spreadsheet, analyzes the data and puts the data back to
the
spreadsheet, similar to RExcel and "R and Calc". Or RKward is only an editor
of
R, working on front-end?
I have tested REx
In my last e-mails, I have asked for help regarding
1. 'defining functions inside loops'
2. 'integrating functions / vector arithmetics'
3. 'vectors out of lists?'
4. 'numerical integration'
Since some of these topics seemed to be relevant (I'm guessing by the # of
replies I got), I'm posting a mo
On Nov 17, 2010, at 5:38 PM, Elliot Joel Bernstein wrote:
I'm trying to make multiple line plots, each with a different color,
using the xyplot command. Specifically, I have an NxK matrix Y and
an Nx1 matrix x. I would like the plot to contain a line for each
(x, Y[,i]), i=1:K. I know some
I have another question about drawing samples from a data frame. This might
sound really tricky. Let me use a data frame I have posted earlier as an
example:
SubIDCSE1 CSE2 CSE3 CSE4 WSE1 WSE2 WSE3 WSE4
1 6 5 6 2 6 22 4
2
I'm trying to make multiple line plots, each with a different color, using the
xyplot command. Specifically, I have an NxK matrix Y and an Nx1 matrix x. I
would like the plot to contain a line for each (x, Y[,i]), i=1:K. I know
something like
xyplot(Y[,1] + Y[,2] + Y[,3] ~ x, type='l')
will w
Basically I'm just looking for a command that can look up the name of the
directory of the script that is running. If I move or copy the script to
another directory, it should be able to read the name of the new directory
without me having to edit the code.
Once I have identified the directory,
Hello, I would like to setup a non-linear constraint for a portfolio using
the portfolioConstraint function. Does somebody now how to do this? thank
you
Felipe Parra
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
http
Hello, I am trying to use portfolioConstraints from the fPortfolio package
and y would like to write a constraint of the form
t(w)*A=z
where w is the weight vector I am optimizing on, A is another vector and a
is a scalar (which can take zero value). Does somebody know how to setup
this constrain
Hi Alex,
Something like this ?
x <- 1:4
y <- list(good=2:5, bad=3:5)
for (yy in y) {
tryCatch( x <- cbind(x, yy),
warning=function(w) cat("problem values: ", yy, "\n")
)
}
Michael
On 18 November 2010 03:19, Alaios wrote:
> Hello when my code executes I receive the message tha
Hello Group,
I am trying to see if there is way to access data that is inside another
namespace.
For e.g. the addATR function in the quantmod package calculates the ATR
using the TTR package and then plots it to the graph.
Now since it has already calculated the info that I need, can I access tha
Echoing an older inquiry
(https://stat.ethz.ch/pipermail/r-help/2008-November/180184.html), does anyone
know of R code that exists to correct for sample selection with panel data as
in:
J.M. Wooldridge (1995), Selection Corrections for Panel Data Models Under
Conditional Mean Independence Ass
-- Forwarded message --
From: Eduardo de Oliveira Horta
Date: Wed, Nov 17, 2010 at 3:59 PM
Subject: Re: [R] Numerical integration
To: David Winsemius
It works, however is not very efficient for the problem I'm working with,
since it is the same as vectorizing the integrand befor
Dear Fellow R Users,
I am experimenting right now the FKF package. I started by working out the
first example included in the package and I am already confused. Would you
offer some kind suggestions? What I want to do is to write down the state
transition equation and the measurement equation
On Nov 17, 2010, at 6:57 PM, Sam Albers wrote:
On Wed, Nov 17, 2010 at 3:49 PM, David Winsemius > wrote:
On Nov 17, 2010, at 6:27 PM, Sam Albers wrote:
Hello all,
I have a fairly simple data manipulation question. Say I have a
dataframe
like this:
dat <- as.data.frame(runif(7, 3, 5))
On Wed, Nov 17, 2010 at 6:00 PM, Alexander Shenkin wrote:
> Hello All,
>
> Searched around, haven't found a decent solution.
>
> I'd like to translate a vector of numbers to a matrix (or to a list of
> vectors) such that the vector values would serve as indicies of the 1's
> in an otherwise-zero-f
On Wed, Nov 17, 2010 at 3:49 PM, David Winsemius wrote:
>
> On Nov 17, 2010, at 6:27 PM, Sam Albers wrote:
>
> Hello all,
>>
>> I have a fairly simple data manipulation question. Say I have a dataframe
>> like this:
>>
>> dat <- as.data.frame(runif(7, 3, 5))
>> dat$cat <- factor(c("1","4","13","1
On 17/11/2010 6:27 PM, Sam Albers wrote:
Hello all,
I have a fairly simple data manipulation question. Say I have a dataframe
like this:
dat<- as.data.frame(runif(7, 3, 5))
dat$cat<- factor(c("1","4","13","1","4","13","13A"))
dat
runif(7, 3, 5) cat
1 3.880020 1
2 4.062800 4
On Nov 17, 2010, at 6:27 PM, Sam Albers wrote:
Hello all,
I have a fairly simple data manipulation question. Say I have a
dataframe
like this:
dat <- as.data.frame(runif(7, 3, 5))
dat$cat <- factor(c("1","4","13","1","4","13","13A"))
dat
runif(7, 3, 5) cat
1 3.880020 1
2 4.
Well, assuming this refers to the histogram function in the lattice package,
looking at the code for print.trellis shows that the default behavior is to
call plot.trellis with the same arguments. So unless you change the default
behavior, there really is no difference between printing and plott
Hello all,
I have a fairly simple data manipulation question. Say I have a dataframe
like this:
dat <- as.data.frame(runif(7, 3, 5))
dat$cat <- factor(c("1","4","13","1","4","13","13A"))
dat
runif(7, 3, 5) cat
1 3.880020 1
2 4.062800 4
3 4.828950 13
4 4.761850 1
David Winsemius wrote:
>
>> Windows binary: not available, see ReadMe
>
> Sorry. Then in answer to your second question:
>
> http://cran.r-project.org/doc/manuals/R-admin.html
>
> Especially:
> http://cran.r-project.org/doc/manuals/R-admin.html#The-Windows-toolset
>
> (At least as I read the t
On Nov 17, 2010, at 6:12 PM, phaaland wrote:
This is what CRAN says about ncdf4:
Windows binary: not available, see ReadMe
Sorry. Then in answer to your second question:
http://cran.r-project.org/doc/manuals/R-admin.html
Especially:
http://cran.r-project.org/doc/manuals/R-admin.html#The-W
On Nov 17, 2010, at 6:00 PM, Alexander Shenkin wrote:
Hello All,
Searched around, haven't found a decent solution.
I'd like to translate a vector of numbers to a matrix (or to a list of
vectors) such that the vector values would serve as indicies of the
1's
in an otherwise-zero-filled matr
This is what CRAN says about ncdf4:
Windows binary: not available, see ReadMe
--
View this message in context:
http://r.789695.n4.nabble.com/ncdf4-for-Windows-R-2-12-0-tp3047807p3047829.html
Sent from the R help mailing list archive at Nabble.com.
__
Alexander -
If I understand your problem, I think this
function will make the matrix you want:
makeyourmatrix = function(vec){
n = length(vec)
mat = matrix(0,n,n)
mat[cbind(1:n,vec)] = 1
mat
}
makeyourmatrix(c(1,3,3,4))
[,1] [,2] [,3] [,4]
[1,]1000
[2,]
On Nov 17, 2010, at 5:56 PM, phaaland wrote:
We need to install the ncdf4 package for R-2.12.0 on Windows. There
are is no
binary available from CRAN
Are you sure? This is the line on CRAN Package Check as of today, 6PM
EST:
ncdf 1.6.3 OK OK OK OK OK OK OK OK OK ERROR OK Brian Ripley
Hello All,
Searched around, haven't found a decent solution.
I'd like to translate a vector of numbers to a matrix (or to a list of
vectors) such that the vector values would serve as indicies of the 1's
in an otherwise-zero-filled matrix (or list of vectors). For example:
> a = c(1,3,3,4)
# pe
We need to install the ncdf4 package for R-2.12.0 on Windows. There are is no
binary available from CRAN so it looks like a manual install will be
required. Is there a guide for doing this?
--
View this message in context:
http://r.789695.n4.nabble.com/ncdf4-for-Windows-R-2-12-0-tp3047807p30478
What do you expect it to look up? How do you want to specify the
directory? Is it supposed to search through some sequence or use ESP?
On Wed, Nov 17, 2010 at 4:08 PM, Cliff Clive wrote:
>
> Hello everyone,
>
> This should be an easy question, I think.
>
> I'd like to write a command in a progr
On Nov 17, 2010, at 5:22 PM, > wrote:
The function is not exported from the package. You have to get
tough with it, e.g.
library(nlme)
summary.lme
Error: object 'summary.lme' not found
But if you insist:
nlme:::summary.lme
function (object, adjustSigma = TRUE, verbose = FALSE, .
Dear R community,
I have slightly unbalanced survey/land elevation data (over 7,000 data
points, grouped according to 60 transects, 30 in each of two
treatments). Observations in both treatments are taken simultaneously,
and the surveying of all 60 transects took 22 days. Comparisons among
t
On Nov 17, 2010, at 5:07 PM, Joanna Geller wrote:
Thanks for the suggestions - I tried all of them, including
require(lme4) and require(nlme) and summary.lme still cannot be found.
I want to use it to get full output of a multi-level model. Maybe it's
a Mac thing?
Any other suggestions would b
The function is not exported from the package. You have to get tough with it,
e.g.
> library(nlme)
> summary.lme
Error: object 'summary.lme' not found
But if you insist:
> nlme:::summary.lme
function (object, adjustSigma = TRUE, verbose = FALSE, ...)
{
fixed <- fixef(object)
...
You ca
Thanks! it worked. this forum is so powerful!
On Tue, Nov 16, 2010 at 12:20 PM, Douglas Bates wrote:
> On Tue, Nov 16, 2010 at 10:30 AM, poko2000
> wrote:
>
> > Hi I am a newbie in R.
> > I have data with dim of 20.
> > How to use lm if i want to do regression with the whole design matrix? My
>
Thanks for the suggestions - I tried all of them, including
require(lme4) and require(nlme) and summary.lme still cannot be found.
I want to use it to get full output of a multi-level model. Maybe it's
a Mac thing?
Any other suggestions would be much appreciated! Thanks, again.
On Wed, Nov 17, 2
Hi all,
after having many images plots on my window, I want to have one single
legend plot.
For that I use image.plot where the position is defined within smallplot:
image.plot(legend.only=T,zlim=c(0,400),
col=rainbow(50),horizontal=T,smallplot=c(.05,.45, .03,.06))
This works fine, bu
Hello everyone,
This should be an easy question, I think.
I'd like to write a command in a program to set the working directory to
whatever directory the file is currently stored in. Suppose I have a file
called "myRscript.r", and it's stored in "C:\Rprojects\myRscript.r", and it
references ot
On Nov 17, 2010, at 4:48 PM, David Winsemius wrote:
On Nov 17, 2010, at 4:39 PM, Joanna Geller wrote:
Hello,
I have installed the nlme package, but every time I try to use the
function "summary.lme,"
How did you "try to use it"?
I get a message that the function cannot be
found.
If i
On Nov 17, 2010, at 4:39 PM, Joanna Geller wrote:
Hello,
I have installed the nlme package, but every time I try to use the
function "summary.lme,"
How did you "try to use it"?
I get a message that the function cannot be
found.
If it says it cannot be found, then you quite possibly forg
That was my thought, and if you are going to be integrating it, you do need to
be concerned with efficiency to some extent, I would imagine.
My experience is that Vecotrize() is theoretically interesting and it is great
for getting you out of a tight spot like this, but if you can avoid it witho
Hello,
I have installed the nlme package, but every time I try to use the
function "summary.lme," I get a message that the function cannot be
found. I've tried to install the package several times and have
re-started R several times, but to no avail. I have also loaded the
lme4 package, in case th
Thanks for the suggestion. The solution below is much better than my
round-about way.
combn(outcomes, 2, list )
I can't do much about the speed of combn() so I wanted to trim the fat
wherever else I could.
C
On 17 November 2010 15:10, Charles C. Berry wrote:
>
> On Wed, 17 Nov 2010, Chris Carl
On Wed, Nov 17, 2010 at 7:41 PM, Partha Sinha wrote:
> 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
On Wed, Nov 17, 2010 at 2:59 PM, johannes rara wrote:
> I have a vector like this:
>
> a <- c("thisIsName", "thisIsAlsoName", "andThisName")
>
> How to break this into pieces and produce a vector with unique parts:
>
> this
> Is
> Name
> Also
> and
> This
Try this:
library(gsubfn)
strapply(a, ".
Great, thanks!
2010/11/17 Henrique Dallazuanna :
> Try this:
>
> Reduce(union, strsplit(gsub("([A-Z])", ";\\1", a), ";"))
>
> On Wed, Nov 17, 2010 at 5:59 PM, johannes rara
> wrote:
>>
>> I have a vector like this:
>>
>> a <- c("thisIsName", "thisIsAlsoName", "andThisName")
>>
>> How to break thi
has anyone accessed Haver from R ?
I have their api for connection to vb , c++ and other softwares...has anyone
tried to use that api with R ?
Any other help in this regard would be much appreciated.
Thanks
Amit
__
R-help@r-project.org mailing list
h
Try this:
Reduce(union, strsplit(gsub("([A-Z])", ";\\1", a), ";"))
On Wed, Nov 17, 2010 at 5:59 PM, johannes rara wrote:
> I have a vector like this:
>
> a <- c("thisIsName", "thisIsAlsoName", "andThisName")
>
> How to break this into pieces and produce a vector with unique parts:
>
> this
> Is
On Wed, 17 Nov 2010, Chris Carleton wrote:
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 r
Hi,
I used sspir for managing non-gaussian State space models but I observed
that for such models only the smoother is gave while the filter is missing.
Why?
--
View this message in context:
http://r.789695.n4.nabble.com/kalman-filter-in-sspir-tp3047486p3047486.html
Sent from the R help mailing
I have a vector like this:
a <- c("thisIsName", "thisIsAlsoName", "andThisName")
How to break this into pieces and produce a vector with unique parts:
this
Is
Name
Also
and
This
-J
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/lis
On Thu, 18 Nov 2010, Paolo Rossi wrote:
> I sent a post to find a clever way to compute a Rolling Average of columns
> in a matrix and I was given the solution below which I am very pleased
> with.
>
>
> RollingAverage <- function(x, RollingObs) {
> cx <- cumsum(x);
> N <- length(x);
> Temp <-
On Nov 17, 2010, at 1:05 PM, Federico Calboli wrote:
Hi All,
I am trying to figure out how to get the position of the knots in a
pspline used in a cox model.
As far as I understand it, the term "knot" should be used with natural
splines, but not for penalized smoothing splines. See p 124
You've tried?
apply(a, 2, RollingAverage, 7)
On Wed, Nov 17, 2010 at 5:27 PM, Paolo Rossi <
statmailingli...@googlemail.com> wrote:
> I sent a post to find a clever way to compute a Rolling Average of columns
> in a matrix and I was given the solution below which I am very pleased
> with.
>
>
>
Dear List,
I am working with a relatively large correlation matrix (~1600*1600), which
I am looking to plot with function cor.plot in package psych.
cor.plot draws a scalebar with as many tick marks/subdivisions as there are
rows in the matrix, which makes the values unreadable even for relativel
I sent a post to find a clever way to compute a Rolling Average of columns
in a matrix and I was given the solution below which I am very pleased
with.
RollingAverage <- function(x, RollingObs) {
cx <- cumsum(x);
N <- length(x);
Temp <- (cx[RollingObs:N] - c(0, cx[1:(N-RollingObs)]))/RollingO
Hello,
Consider this little function (not concerned about efficiency here)
m=function(r) { a<-serialize(r,NULL); zz<-rawConnection(raw(0),"r+");
writeBin(length(a),zz); c(rawConnectionValue(zz),a)}
a=m(10)
p=rawConnection(a,"r+"); u=readBin(p,"int")
unserialize(p)
Error in unserialize(p) : unk
On Wed, Nov 17, 2010 at 08:19:53PM +0200, arturs.onz...@gmail.com wrote:
> 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 :
>
> 1289008
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
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
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
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
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
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
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
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
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
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,
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
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
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
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
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
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
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
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.
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 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))
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
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
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)))
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
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
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
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
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
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
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
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
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(
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
1 - 100 of 166 matches
Mail list logo
