Duncan Murdoch wrote:
> On 15/04/2008 7:20 PM, Hongshu Chen wrote:
>> Hi,
>>
>> I am building a package which contains some Java class files. I am planning
>> to create a subdirectory "inst/Java/" in the source directory and put Java
>> files in it, therefore after installing the package, those fil
The issue is that the X11() device is the default device only in
interactive use, and further than the event loop is only run in
interactive use.
You need to
1) make use of R 2.7.0 RC
2) call R with R --vanilla --interactive
3) arrange for your script to wait for something so the plot remains u
Hi,
as already mentioned, sorting could be a pain.
My solution to that is to write my own "order" routine for a given
language.
The idea is to transform the UTF-8 string into ASCII in such a way
that the built-in order routine outputs the desired result. But this
could be a very stony way.
Hi,
I would like to have at least 2 x11 windows open and make many graphs into
them. I wonder how to draw the graph in a certain device so that the new
graph will be in the place of the old one.
For example I have open x11 with device numbers 1 and 2. I want to make
plot to the device 1 without d
Hi, sorry for jumping in here, but to me your description of why you want to
have only the needed data rows remains ambiguous.
If you just want to select the data you indicate then you do:
selected.data=data[ , needed=="yes"]
where "data" is the name of your long dataset (13 obs). As I see it, y
Here is the original dataset I am working on:
> myTime
Time Series:
Start = 0
End = 1066
Frequency = 0.0769230769230769
[1]0 13 26 39 52 65 78 91 104 117 130 143 156
169 182
[16] 195 208 221 234 247 260 273 286 299 312 325 338 351
364 377
[31] 390 403
Patrick,
my intention was, to perform a one-to-one exact match, which pairs each treated
unit with ONE control unit (without replacement), using my two confounders
(age, school) for matching.
Patrick Connolly schrieb:
On Mon, 14-Apr-2008 at 08:37AM +0200, Udo wrote:
|> Zitat von Peter Alspach <[
Hi use Rs,
I have a csv file:
"1989-90","1990-91"
Barley,23,34
Oats,15,16
Which I want to turn into:
year, Barley, Oats
1 "1989-90", 23, 15
2 "1990-91",34,16
Transpose doesn't quite do it, is there a standard way?
Cheers,
Geoff Russell
__
R-he
On 16-Apr-08 10:21:00, Tommi Viitanen wrote:
> Hi,
> I would like to have at least 2 x11 windows open and make
> many graphs into them. I wonder how to draw the graph in
> a certain device so that the new graph will be in the place
> of the old one.
>
> For example I have open x11 with device numb
?dev.set (on the same page as dev.off).
On Wed, 16 Apr 2008, Tommi Viitanen wrote:
> Hi,
>
> I would like to have at least 2 x11 windows open and make many graphs into
> them. I wonder how to draw the graph in a certain device so that the new
> graph will be in the place of the old one.
>
> For e
On Wed, 16 Apr 2008, Guohui Ding wrote:
> Dear All,
>
> I have written a function in c (in a file named 'a.c') as follows,
*PLEASE* don't claim credit for the work of others (in this case it is me
that you owe the apology and retraction to). That is from 'S Programming'
p. 126, and 'Writing R
Hello everybody,
I want to save a graphic using the "pdf environment", everything in order,
but I can not
to save the graph in "horizontal" way. I have read the manual (to plot,
par), the FAQ, google, but nothing I found that can help me. Anyone can say
me a hint?
Thanks a lot,
--
Josue
Dear All,
I have written a function in c (in a file named 'a.c') as follows,
void convolve(double *a, int *na, double *b, int *nb, double *ab)
{
int i, j, nab = *na + *nb - 1;
for(i = 0; i < nab; i++)
ab[i] = 0.0;
for(i = 0; i < *na; i++)
for(j = 0; j < *nb; j++)
ab[i + j] += a
On 16 Apr 2008, at 12:21, Tommi Viitanen wrote:
> For example I have open x11 with device numbers 1 and 2. I want to
> make
> plot to the device 1 without doing anythin to the 2 and not making a
> new
> x11. Something like ?:
Do you mean something like dev.set(DEVICENUMBER) ?
Have a look at
Hi
I would like to know whether it is possible to implement the concept of "weak
disposability" in FEAR package for Data envelopment analysis. Weak
disposability means equality constraint
with warm regards
Ramelan
___
Yahoo! For Go
Georg Ehret wrote:
> Dear R community,I am creating large graphs with hundreds of
> thousands of datapoints. My usual way for output was pdf, but now I am
> getting file sizes of >30Mb that do not open well (or at all) in Adobe. Is
> there a way to reduce the resolution or get rid of overla
Ted Harding wrote:
>
> One is that "NA" is not a value. Its logical status is,
> in effect, "value not known". Therefore, when 'y' is "NA",
> "x==y" cannot have a definite resolution, since it is
> possible for the unkown value of 'y' to be equal to the
> value of 'x'; and equally possible that
Thanks Hans!
Hans-Joerg Bibiko wrote:
> Hi,
>
> as already mentioned, sorting could be a pain.
>
> My solution to that is to write my own "order" routine for a given
> language.
> The idea is to transform the UTF-8 string into ASCII in such a way
> that the built-in order routine outputs the des
Try this:
x <- read.table(textConnection("\"1989-90\",\"1990-91\"
Barley,23,34
Oats,15,16"), sep = ",", header = T, check.names = FALSE)
t(x)
or if you want year as column and not rownames:
data.frame(Year = names(x), t(unname(x)))
On Wed, Apr 16, 2008 at 7:27 AM, Geoff Russell <[EMAIL PROTEC
See in ?pdf function, 'paper' argument.
On Wed, Apr 16, 2008 at 7:56 AM, Josué Polanco <[EMAIL PROTECTED]> wrote:
> Hello everybody,
>
> I want to save a graphic using the "pdf environment", everything in
> order,
> but I can not
> to save the graph in "horizontal" way. I have read the manual (
Many thanks -- the "unname" makes all the difference.
Geoff.
On 4/16/08, Henrique Dallazuanna <[EMAIL PROTECTED]> wrote:
> Try this:
>
> x <- read.table(textConnection("\"1989-90\",\"1990-91\"
> Barley,23,34
> Oats,15,16"), sep = ",", header = T, check.names = FALSE)
>
> t(x)
>
> or if you want y
Hi,
I have a R script that loads an image R.data, does some operations, then
save to the R.data again. Suppose I have done some computation before
loading the R script, will all the objects before the R script execution be
saved to R.data ? If yes, how can I specify save.image to save only those
o
See ?save
On Wed, Apr 16, 2008 at 8:46 AM, Ng Stanley <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I have a R script that loads an image R.data, does some operations, then
> save to the R.data again. Suppose I have done some computation before
> loading the R script, will all the objects before the R scr
Read and reread, can't make out. Will try an experiment later
On Wed, Apr 16, 2008 at 7:51 PM, Henrique Dallazuanna <[EMAIL PROTECTED]>
wrote:
> See ?save
>
> On Wed, Apr 16, 2008 at 8:46 AM, Ng Stanley <[EMAIL PROTECTED]>
> wrote:
>
> > Hi,
> >
> > I have a R script that loads an image R.
Thank you for your response. I tried searching the archives but did not
know that it was called X11() device.
Can this be done in version 2.6.2 or do I need to have version 2.7.0
running? Also, when you say wait in my script, I need to put a sleep or
something in the script or else it will delet
:) Thanks.
I owe you an apology and retract it, as I copied the codes from others and
did not check the origin of these codes.
2008/4/16, Prof Brian Ripley <[EMAIL PROTECTED]>:
>
> On Wed, 16 Apr 2008, Guohui Ding wrote:
>
> Dear All,
> >
> > I have written a function in c (in a file named 'a.c')
Thanks Mark!
Mark Leeds wrote:
> Hi Ricardo: you can look at ?regexpr but I'm pretty sure that below returns
> a vector of TRUES and FALSES depending on whether the condition ( is Jalapa
> contained in unidad ) is TRUE. So, wherever the vector is TRUE is where
> Jalapa was contained in unidad. So,
The Brown-Forsyth test is included in the HH library in function hov().
hov is "homogeneity of variance".
Use
hov(y ~ x, data=mydata)
and
plot.hov(y ~ x, data=mydata)
Rich
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/
Charles C. Berry wrote:
>
> Like this:
>
> conaguaMexicoSub <- subset(conagua,
> regexpr("Jalapa", as.character( unidad ) ) != -1
> select = c(equipo,X101:X309))
>
> But it looks like you have to sort out some issues concerning dueling
> locales first.
>
> Chuck
>
>
>
Far clear! Thanks
Mark Leeds wrote:
> Hi Ricardo: I'm glad it helped but check
>
> regexpr("Jalapa", as.character( unidad ) ) != -1
>
> because I think it should be
>
> regexpr("Jalapa", as.character( DF$unidad ) ) != -1
>
> if unidad is a column of DF.
Checked! attach(DF) did the trick! :-)
Greetings,
Ricardo
I have problems with emf files because OpenOffice does
a poor job at importing figures in this format. On the other
hand, imported eps figures are not displayed, just printed.
Is there any R graphic gui able to export as SVG (or
other vector format)? I've
tried JGR, iplots and svGUI.
Thanks!
Agus
Using xyplot, I am plotting observed versus predicted values conditional to the
values of a third variable.
I would thus like to have an aspect ratio of 1 and the same axis range.
Since the range of the values differs substantially between panels, the axis
ranges shall be different for each pane
Hans-Joerg Bibiko wrote:
> Hola,
>
> Muchas gracias!
> This is new to me. I learnt Spanish a bit - well - 20 years ago ;)
> But this simplifies it.
This change happens just 14 years ago! You you are not guilty!
>
>
> Recuerdos
>
> Hans
>
Saludos cordiales! "Read" you in Spanish whenever you want!
On Wed, 16 Apr 2008, Pologruto, Thomas wrote:
> Thank you for your response. I tried searching the archives but did not
> know that it was called X11() device.
>
> Can this be done in version 2.6.2 or do I need to have version 2.7.0
Flag --interactive is only in 2.7.0. There is no simple way to
On Wed, 16 Apr 2008 13:20:35 +0200 Agustin Lobo wrote:
AL>I have problems with emf files because OpenOffice does
AL>a poor job at importing figures in this format. On the other
AL>hand, imported eps figures are not displayed, just printed.
AL>Is there any R graphic gui able to export as SVG (or
AL>
Thank you very much. Most appreciated.
Please follow the attached hyperlink to an important disclosure:
http://www.credit-suisse.com/legal/marketcommentary
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: Wednesday, April 16, 2008 8:52 AM
To: Pologruto, T
R 2.7.0 RC has an svg() device on Unix-alikes, and there are RSvgDevice
and RSVGTipsDevice packages on CRAN.
I've seen so many problems with SVG renderers (especially related to
fonts) that I would prefer not to rely on it.
If this is Windows, the preferred route seems to be to use EPS with a
Prof Brian Ripley wrote:
>
> Yes, but only for the same character set. I believe R.app assumes
> UTF-8, and I would not expect to be able to change charset on a
> running console.
>
> Please do use R-sig-mac for MacOS-specific issues.
I've not a clear idea yet about the relationship between loc
Thanks,
I know about saving to file via the graphic devices, and I
do use the cairoSVG for example. But I'm looking
for a minimum interactivity, such in the windows Gui where
you resize your graphic and when you like it
you click on "save as..." The problem is that you only
have emf as vector optio
I've just tried the following in 2.6.2
auk% cat x11.R
X11()
interactive()
plot(1:10)
locator()
q()
explicitly launch a device. It seems to be responsive enough (at least on
my system) and will stay up until you middle- or right-click on the device
surface. In earlier threads others have repor
Hi Thomas,
> I have looked through several R books and searched the web to find
> answers to my questions with no results. I have an ensemble of time
> series data (what are essentially Monte Carlo simulations) which I would
> like to summarize as a time series of boxplots, by date/time at 6-hr
Hello list,
I am performing a sensitivity analysis using the package ROCR. I am using
the class "prediction" in this aim. My question is, could anyone tell me
what the vector "cutoffs" represent in the result?
Thank you all,
Eleni
[[alternative HTML version deleted]]
___
All,
A new Special Interest Group devoted to the use or R in ecology has been
created. The list came about through a very thorough and scientific
survey that identified a need (i.e. I chatted with Sarah Goslee at a
recent conference and she agreed that this might be a good idea).
Myself and Sara
Eleni Christodoulou gmail.com> writes:
> I am performing a sensitivity analysis using the package ROCR. I am using
> the class "prediction" in this aim. My question is, could anyone tell me
> what the vector "cutoffs" represent in the result?
>
As the docs say:
cutoffs:
A list in which each
Dear Fehmi,
Curiously, I check the Rcmdr package on an up-to-date Ubuntu 7.10
system running R 2.6.2. I've just checked again after updating both
Ubuntu and all R packages, and I can't duplicate your problem -- the
"Tools -> Load packages" menu works fine for me, as does everything
else that I che
On 16-Apr-08 11:13:40, Alberto Monteiro wrote:
>
> Ted Harding wrote:
>>
>> One is that "NA" is not a value. Its logical status is,
>> in effect, "value not known". Therefore, when 'y' is "NA",
>> "x==y" cannot have a definite resolution, since it is
>> possible for the unkown value of 'y' to b
> mypal <- trellis.par.get('superpose.line')$col[1:2]
> mycol <- mypal[1]
> trellis.par.set(list(dot.symbol = list(pch = 16, cex = 1.5)))
> dotplot(1:6 ~ audRating, data = txtA.ci,
+ scales = list(y = list(at = 1:6, labels = txtA.ci$pitchAud, col =
'black', cex = 1.1)),
+ ylab = '
John Fox wrote:
> Dear Fehmi,
>
> Curiously, I check the Rcmdr package on an up-to-date Ubuntu 7.10
> system running R 2.6.2. I've just checked again after updating both
> Ubuntu and all R packages, and I can't duplicate your problem -- the
> "Tools -> Load packages" menu works fine for me, as doe
On 4/16/08, Andreas Krause <[EMAIL PROTECTED]> wrote:
>
> Using xyplot, I am plotting observed versus predicted values conditional to
> the values of a third variable.
> I would thus like to have an aspect ratio of 1 and the same axis range.
> Since the range of the values differs substantially
Using the "ergoStool" data cited in Mixed-Effects Models in S and
S-PLUS by Pinheiro and Bates as an example, we have
> library(nlme)
> fm <- lme(effort~Type-1, data=ergoStool, random=~1|Subject)
> summary(fm)
Linear mixed-effects model fit by REML
Data: ergoStool
AIC BIC
Does R have a function for doing an ARAR analysis of a time series?
Tom
--
View this message in context:
http://www.nabble.com/ARAR-Time-Series-Analysis-in-R-tp16722577p16722577.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@
Hi
I have a general statistics question on calculating confidence interval of
log transformed data.
I log transformed both x and y, regressed the transformed y on transformed
x: lm(log(y)~log(x)), and I get the following relationship:
log(y) = alpha + beta * log(x) with se as the standard error
A reasonable question here might be, what do you intend to do with the
SVG file, once you have it. OO.org does not natively import SVG files
and the only import filter that I am aware of for it:
http://wiki.services.openoffice.org/wiki/SVG_Import_Filter
is an incomplete implementation. This me
At 18:03 14/04/2008, Ryan Lauritsen wrote:
>Hi all. I'm brand new to R.
>
>My dataset (stored in MySQL) is a list of weather stations in rows by
>year with various weather variables in columns, for example:
>STNID YEAR TEMP DEWP
>station11990 54 50
>station11991 23 10
As alluded, I've got several off-line replies and suggestions,
notably Tobias Verbeke mentioning Matlab code {as an example of "pseudo
code" :-) } by Burkhardt on the web,
and Barry Rowlinson providing code similar to Peter Dalgaard's,
see below.
I've found myself python code on the web, unfortunat
tom soyer wrote:
> Hi
>
> I have a general statistics question on calculating confidence interval of
> log transformed data.
>
> I log transformed both x and y, regressed the transformed y on transformed
> x: lm(log(y)~log(x)), and I get the following relationship:
>
> log(y) = alpha + beta * l
See comments inline:
Agustin Lobo wrote:
> Thanks for your answers, let me summarize the situation:
>
> 1. I've tested the SVG import for OO, it's not great
> but it is much better than what you get for EMF. I know
> nothing about SVG, though, perhaps it is not the
> best graphics format.
In the
No. The relationship U = V doesn't have to hold for positive-semidefinite
matrices, just as it doesn't for an indefinite matrix (i.e. a matrix with
both positive and negative eigenvalues), since you can have u_i = +/- (v_i
)^T corresponding to the zero eigenvalue.
You may refer to, for example,
save.image() is a wrapper for save() (Type 'save.image' (without quotes) at the
prompt to see the code for the function.)
Try something like this (not tested):
> load("oldvars.RData")
> old.vars <- ls(all=TRUE)
> ... computation ...
> save(list=setdiff(ls(all=TRUE), old.vars), file="newvars.RDat
Rubén Roa-Ureta wrote:
> tom soyer wrote:
>> Hi
>>
>> I have a general statistics question on calculating confidence interval of
>> log transformed data.
>>
>> I log transformed both x and y, regressed the transformed y on transformed
>> x: lm(log(y)~log(x)), and I get the following relationship:
Richard M. Heiberger wrote:
> The Brown-Forsyth test is included in the HH library in function hov().
> hov is "homogeneity of variance".
> Use
>hov(y ~ x, data=mydata)
> and
>plot.hov(y ~ x, data=mydata)
>
>
Are you guys talking about the same thing? One seems to be talking about
a mod
Hello,
I asked this question a little while ago (
https://stat.ethz.ch/pipermail/r-help/2008-April/158761.html ) but got no
response. Can anybody explain to me the difference between /, :, and %in% in
the definition of random effects in lmer, such as:
(1|A/B), (1|A:B), (1|B %in% A)?
My underst
R may not be giving you what you want, but it is doing the right thing.
You can change what the base category is through contrasts but you can't
get the marginal effects for every level of all factors because this
creates a linear dependence in the model matrix.
> -Original Message-
> Fro
Thanks for your answers, let me summarize the situation:
1. I've tested the SVG import for OO, it's not great
but it is much better than what you get for EMF. I know
nothing about SVG, though, perhaps it is not the
best graphics format.
2. I need the figure displayed on the screen, I'm using
Impre
Hi all,
I've read the R for windows FAQ and am a little confused re:
memory.limit and memory.size
to start using R 2.6.2 on WinXP, 2GB RAM, I have the command line "sdi
--max-mem-size=2047M"
Once the Rgui is open, memory.limit() returns 2047, memory.size()
returns 11.315, and memory.size(max=T)
Survfit returns an "Error in temp[, 1] : incorrect number of
dimensions" error message when I attempt to use newdata to validate my
models. The syntax I am using is "sfit <- survfit(mod1,
newdata=testd, individual=T)". The dataframe testd has two columns
and 39 rows (see below). The rows are th
Hi all,
I've read the R for windows FAQ and am a little confused re:
memory.limit and memory.size
to start using R 2.6.2 on WinXP, 2GB RAM, I have the command line "sdi
--max-mem-size=2047M"
Once the Rgui is open, memory.limit() returns 2047, memory.size()
returns 11.315, and memory.size(max=T)
Hi All,
I am a relative newbie at R and I am trying to plot some intraday price data
over a long period of time. I would like to use a weekday(only) scale for
x-axis. Can anyone tell me how to go about this?
When I plot just the prices it looks fine, but when I add in the Dates/Times
vs prices
At the risk of annoying Doug Bates, I'll point out that "glht" in
the multcomp package works with lmer objects. In fact, you can
supply your own degrees of freedom value via the "df" argument
(a fact which is not immediately obvious in the glht help page).
If you don't supply a df value, it will
Hello all,
I should probably know this by now... Anyway:
I have a large matrix (dim(data) is 3000 18000). In each element are
one of the following character strings "0/0", "1/1", "1/2", "2/2". I
wanted to replace "0/0" with NA and the other three with 0,1,2
respectively. To accomplish just the f
Hi All,
I am a relative newbie at R and I am trying to plot some intraday price data
over a long period of time. I would like to use a weekday(only) scale for
x-axis. Can anyone tell me how to go about this?
When I plot just the prices it looks fine, but when I add in the Dates/Times
vs prices
On Wed, 16 Apr 2008, Doran, Harold wrote:
> R may not be giving you what you want, but it is doing the right thing.
> You can change what the base category is through contrasts but you can't
> get the marginal effects for every level of all factors because this
> creates a linear dependence in the
You might find the width and height arguments useful also.
HTH
On Wed, 16-Apr-2008 at 08:19AM -0300, Henrique Dallazuanna wrote:
|> See in ?pdf function, 'paper' argument.
|>
|> On Wed, Apr 16, 2008 at 7:56 AM, Josu? Polanco <[EMAIL PROTECTED]> wrote:
|>
|> > Hello everybody,
|> >
|> > I want
On 17/04/2008, at 7:52 AM, Matthew Keller wrote:
> Hello all,
>
> I should probably know this by now... Anyway:
>
> I have a large matrix (dim(data) is 3000 18000). In each element are
> one of the following character strings "0/0", "1/1", "1/2", "2/2". I
> wanted to replace "0/0" with NA and th
On Thu, 17 Apr 2008, Rolf Turner wrote:
>
> On 17/04/2008, at 7:52 AM, Matthew Keller wrote:
>
>> Hello all,
>>
>> I should probably know this by now... Anyway:
>>
>> I have a large matrix (dim(data) is 3000 18000). In each element are
>> one of the following character strings "0/0", "1/1", "1/2"
On 17/04/2008, at 9:33 AM, Charles C. Berry wrote:
> I'll lay odds that Matthew's 'matrix' is actually a data.frame, and
> I'll not be surprised if the columns are factors.
I suspect that you're right.
***Why*** can't people distinguish between data frames and matrices?
If
I think the clue is that the message you quote comes from gstat, which
does not use R's memory allocator. It is gstat and not R that has failed
to allocate memory.
Try re-reading the help page for memory.size. 'max=T' does not indicate
the limit (that is the job of memory.limit()), but the ma
Yes Chuck, you're right.
Thanks for the help. It was a data.frame not a matrix (I had called
as.matrix() in my script much earlier but that line of code didn't run
because I misnamed the object!). My bad. Thanks for the help. And I'm
VERY relieved R isn't that inefficient...
Matt
On Wed, Apr 16
If z is your zoo object then you plot it like this:
plot(zoo(coredata(z)), xaxt = "n")
and then add a custom axis using axis or Axis. Also the quantmod
package has this type of graphics built in.
On Wed, Apr 16, 2008 at 3:57 PM, rsan <[EMAIL PROTECTED]> wrote:
>
> Hi All,
>
> I am a relative ne
I think any geostatistical program/R package would have trouble handling
12000 observations on a PC. The empirical variogram would be built with
the combinations of 12000 over 2 pairs, nearly 72 millions pairs, and
during kriging, if you didn't restrict the search neighbourhood,
interpolation w
Hi folks,
when trying to update packages in version 2.6.1 I recuieve the following
error messages:
Warning: unable to access index for repository
http://www.sourcekeg.co.uk/cran/bin/windows/contrib/2.6
Warning: unable to access index for repository
http://www.stats.ox.ac.uk/pub/RWin/bin/windows
Hi all
I am trying to calculate the instantaneous hazard (rather than
cumulative hazard) for an event "bird death".
Some sample data is as follows:
year<-c(rep(2003,10),rep(2004,10))
time<-c(5,10,5,15,20,10,15,15,25,25,5,15,15,5,25,10,25,15,35,25)
event<-c(1,0,1,1,0,1,1,0,1,0,1,0,0,1,0,0,0,0,1,0
You're right-- I think we're not talking about the same thing. It
looks to me like the hov() function in the HH package is a test for
whether or not the variances are different. What I was looking for
was a test (also by Brown-Forsyth, and referred to as the F* statistic
in my textbook) for testi
Oddly enough, the variogram modelling is rather quick in Surfer, but
one cannot compute the standard errors. I restricted the search to the
approximate the range of the variogram model (2000m). I can get R to
compute with 12079 observations, but 13453 and I run into the gstat
error message.
Hello,
I am trying to use R2WinBUGS to conduct a mixed treatment comparison (MTC)
analysis. On the surface, it seems to me that I am following the correct
steps: (1) reading the data into R, (2) specifying initial values for the
parameters in the model and (3) fitting the model
I have several old encapsulated postscript graphics generated by R that
I want to include in a LaTeX Beamer presentation so I want them in pdf.
Under Ubuntu linux typing
convert z.eps z.pdf (to use ImageMagick)
results in a nice looking graphic but the bounding box is lost. Doing
convert z.
Frank E Harrell Jr wrote:
> I have several old encapsulated postscript graphics generated by R that
> I want to include in a LaTeX Beamer presentation so I want them in pdf.
> Under Ubuntu linux typing
>
> convert z.eps z.pdf (to use ImageMagick)
>
> results in a nice looking graphic but the
ghement.ca> writes:
>
>
>Hello,
>
>I am trying to use R2WinBUGS to conduct a mixed treatment comparison (MTC)
>analysis.
>Thank you in advance for any help you may be able to provide!
>
Can you please provide your model.bug file?
Ben Bolker
__
This maybe a silly question. I'm trying to figure out a way to draw a line from
a data set which contain NA. Say, I have a set of data as:
x <- c(1.1 2.2 NA 4.4 5.5) ; y <- c(1:5) # as x,y of point a, b, c, d, and e.
I would like to plot this to a line by using dot-line to connect the two
adjac
Hi Folks... No code to troubleshoot here. I need some suggestions about the
right terminology to use in further searching, and any suggestions about R
pkgs that might be appropriate.
I am in the planning stages of a project in which IR, NMR and other spectra
(I'm a chemist) would be collected on
Marc Schwartz wrote:
> Frank E Harrell Jr wrote:
>> I have several old encapsulated postscript graphics generated by R
>> that I want to include in a LaTeX Beamer presentation so I want them
>> in pdf. Under Ubuntu linux typing
>>
>> convert z.eps z.pdf (to use ImageMagick)
>>
>> results in a n
On 17/04/2008, at 12:56 PM, DinoDragon wrote:
> This maybe a silly question. I'm trying to figure out a way to draw
> a line from a data set which contain NA. Say, I have a set of data as:
>
> x <- c(1.1 2.2 NA 4.4 5.5) ; y <- c(1:5) # as x,y of point a, b,
> c, d, and e.
>
> I would like to
Try glht in package multcomp.
Simon.
On Wed, 2008-04-16 at 12:00 -0400, Gang Chen wrote:
> = 1), c(TypeT2 = 1))) :
> Only defined for lm,glm objec
--
Simon Blomberg, BSc (Hons), PhD, MAppStat.
Lecturer and Consultant Statistician
Faculty of Biological and Chemical Sciences
The University of
Dear Mr,
I saw by the web your request about the code in Language C of the method
Nelder-Mead with constraints,
I would like if you got the program, it's posible you can help us to get the
code in C,
I'm grateful for your support,
Regards,
Jorge Angarita
[[alternative HTML version de
Dear R helpers,
I was doing a genetic project with two datasets X and Y. There are
some IDs in both data sets, and others in either data set. I used
"merge(x,y,by="ID",all=TRUE)". The data set Y contains a variable (a
genotype) which is also in data X. When I merge X with Y, these two
variables we
Try this which interpolates the NAs using na.approx from zoo
drawing the interpolated series with lty=3 and then overwrites
solid lines with the gaps using lines.
library(zoo)
plot(y ~ x, as.data.frame(na.approx(zoo(cbind(x, y, lty = 3, type = "l")
lines(y ~ x)
On Wed, Apr 16, 2008 at 8:56 P
Hi,
Design isn't strictly an R base package, but maybe someone can explain
the following.
When lrm is called within a function, it can't find the dataset dd:
> library(Design)
> age <- rnorm(30, 50, 10)
> cholesterol <- rnorm(30, 200, 25)
> ch <- cut2(cholesterol, g=5, levels.mean=TRUE)
>
On Wed, 16 Apr 2008, Lu, Jiang wrote:
> Dear R helpers,
>
> I was doing a genetic project with two datasets X and Y. There are
> some IDs in both data sets, and others in either data set. I used
> "merge(x,y,by="ID",all=TRUE)". The data set Y contains a variable (a
> genotype) which is also in dat
Greetings all,
I am trying to use Deal to learn a Bayesian Network for discrete data. I
came across the following problem with jointprior function.
>library(deal)
> a <- read.csv("prepared.cluster1.csv")
> fit <- network(a)
> fit.prior <- jointprior(fit)
Error in array(1, Dim) : 'dim' specifies t
Hello,
I've got some problems. I hope someone can help me.
First question:
I am trying to get grid on scatterplot3d (from scatterplot3d package). It
seems that scatterplot3d draw on grip on X and Z side. Is it possible to get
Grid on the whole Box?
Second question:
Is it possible to use the st
1 - 100 of 103 matches
Mail list logo
