__
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.
Hi,
I have a series of id's with multiple visits and questionnaire scores. This
is a clinical trial that will be analyzed using the last observation carried
forward method. In other words, in order to comply with intent to treat
analysis when many subjects withdraw, data points for the last visit
On 09/10/11 10:39, Duncan Murdoch wrote:
On 11-10-08 5:32 PM, Rolf Turner wrote:
On 09/10/11 00:18, Duncan Murdoch wrote:
On 11-10-07 5:26 PM, Carl Witthoft wrote:
Just wondering here -- I tested and found to my delight that
% round(325.4,-2)
[1] 300
gave me exactly what I would have expect
On 09/10/11 00:18, Duncan Murdoch wrote:
On 11-10-07 5:26 PM, Carl Witthoft wrote:
Just wondering here -- I tested and found to my delight that
% round(325.4,-2)
[1] 300
gave me exactly what I would have expected (and wanted). Since it's not
explicitly mentioned in the documentation that neg
I evaluated the Bayes factor in the k=2 allele case with a "triangular"
prior under the null as in the example in the help file:
HWETriangBF2(nvec=c(88,10,2))
[1] 0.4580336
When I swap the n11 entry and n22 entry of nvec, I received totally
different Bayes factor:
>
> HWETriangBF2(nvec=c(2,
By going through the code, I came to a similar conclusion also: it
seems the match.2.list function must be modified in the following way
to make it work:
match.2.list <- function(l1,l2){
if (class(l1)=="element"){
l1 <- l1[[1]]
l2 <- l2[[1]]
}
length(l1)==length(l2)
}
because depending wh
I'm not quite familiar with R. Could you recommend some relative material.
Actually, I'm not intend to loop from 1 to 1. I'm just trying to see if it
works. It's supposed to be some len>1.
I'm trying to solve the following example:
there are 3 types of variables(x,w,t), they are not independen
Yeah thanks guys. I can utilize this update function to help clean things
up. It's not as simple as I was hoping though. Also, Bert, by updating the
data argument, don't I still have to update the model anyway to take
advantage of the extra covariates?
Justin
On Wed, Oct 5, 2011 at 12:55 PM, Be
On 10/8/11 6:11 PM, (Ted Harding) wrote:
Carl Witthoft's serendipitous discovery is a nice example
of how secrets can be guessed by wondering "what if ... ?".
So probably you don;t need to tell the secrets.
Taking the "negative digits" to their logical extreme:
round(654.321,2)
# [1] 6
One last comment: add sep="," to the read.csv assignment, to avoid "Error in
scan(file ... : scan() expected 'a real', got" -type errors:
x<-matrix(scan("/home/victoria/R/5.3.2.2.csv",n=5*12,sep=","),5,12,byrow=TRUE)
--
View this message in context:
http://r.789695.n4.nabble.com/Entering-data-
Hi Francy,
A bootstrap test would likely be sufficient for this problem, but a
one-sample t-test isn't advisable or necessary in my opinion. If you use a
t-test multiple times you are making assumptions about the distribution of
your data; more importantly, your probability of Type 1 error will b
The Benchmarking package is generally working well for me but one of its
functions, dea.plot.transform, produces results that look odd to me.
Under Ubuntu 10.04 LTS, I am running R version 2.13.2 and Benchmarking
version 0.20 ($Revision: 118 $). In that environment, I tried executing
the follo
On 08-Oct-11 21:39:07, Duncan Murdoch wrote:
> On 11-10-08 5:32 PM, Rolf Turner wrote:
>> On 09/10/11 00:18, Duncan Murdoch wrote:
>>> On 11-10-07 5:26 PM, Carl Witthoft wrote:
Just wondering here -- I tested and found to my delight that
% round(325.4,-2)
[1] 300
gave
On 11-10-08 5:32 PM, Rolf Turner wrote:
On 09/10/11 00:18, Duncan Murdoch wrote:
On 11-10-07 5:26 PM, Carl Witthoft wrote:
Just wondering here -- I tested and found to my delight that
% round(325.4,-2)
[1] 300
gave me exactly what I would have expected (and wanted). Since it's not
explicitl
On 11-10-08 11:04 AM, Ben qant wrote:
Thank you!
Sorry, I have a couple more questions:
1) How to I turn off the box shading completely? I figured out how to
lighten it up to a grey color with col=c("white"... . I looked in the
package pdf, but I didn't see anything.
2) I read about the zunit i
Update:
1. Here is a better illustration. I assumed that the gene expression was
the same for the three replicates for each gene, in the respective Time and
Treatment groups, so that I could follow (verify) what was going on.
2. To better-view the following, copy and paste into a plain-text edit
Thanks, David.
By going through the code, I came to a similar conclusion also: it seems the
match.2.list function must be modified in the following way to make it work:
match.2.list <- function(l1,l2){
if (class(l1)=="element"){
l1 <- l1[[1]]
l2 <- l2[[1]]
}
length(l1)==length(l2)
}
because d
Hi everybody,
I got a solution for my problem with Eric Paniagua (many thanks!)! Thank you
a lot everybody!
Paniagua said explained that I was losing row.names when I merged a
data.frame with a factor object (data2$color), because a factor objects
don't have row.names. Instead of merging data1 to
*Dear All*
*
*
*I am able to compile the simple program given in the site
http://dirk.eddelbuettel.com/code/rinside.html, but when i execute it a
dialog appear saying "Fatal Error: unable to load base package". Please
help.*
*
*
*With regards*
*
*
*Sudhabrata Majumder*
*
*
This is the source code
Hi Duana,
I am looking for the same information, I am wondering if you were able to
change the line type. If you did, I will apreciate if you can give me a hand
with it.
Thanks,
Juan
--
View this message in context:
http://r.789695.n4.nabble.com/line-type-in-plots-of-the-ade4-package-tp3312465p3
Hi,
I am having trouble understanding how to approach a simulation:
I have a sample of n=250 from a population of N=2,000 individuals, and I
would like to use either permutation test or bootstrap to test whether this
particular sample is significantly different from the values of any other
rando
Thanks!
Yes, that produces tables that are formatted in the same way as the gam
output. I'm hoping to have publication-ready tables that have standard
errors in parentheses listed below coefficients. Do you know of a method to
do that?
David
--
View this message in context:
http://r.789695.n4
Hallo Mohammad,
Under Linux (OpenSuse 11.3 x64) i'v just installed Mesa + developer and
source packages for Mesa:
zypper install Mesa
zypper install Mesa-32bit
zypper install Mesa-debuginfo
zypper install Mesa-debuginfo-32bit
zypper install Mesa-devel
zypper install Mesa-devel-32bit
I would s
"Simple" question: how can I do that?
I think that I need to find a possible way to save my googlevis motion plot
into a SWF file, but how?
Thanks for the help,
Roberto.
--
View this message in context:
http://r.789695.n4.nabble.com/R-GoogleVis-I-need-to-import-my-chart-into-powerpoint-file-tp3
Thank you,
you are right my execution was completly wrong. it was
plot(0,type="n"); rect(4,0,6,11,col=5); plot(x,y).
Thanks for your help :)
H.
On Sat, Oct 8, 2011 at 12:43 PM, Robert Baer wrote:
> This is what you tried? What doesn't work?
> x=1:10
> y=1:10
>
> plot(x,y, type='n')
> rect(4,0,
To run a cmd.exe built-in function use 'cmd /c function arg1 ...', as in
> status <- system("cmd /C echo foo")
foo
> status # 0 means all went well
[1] 0
as opposed to
> status <- system("echo foo")
> status
[1] 127
You get a little more information about that 127 return value
DEL is built in to CMD. Use the R "unlink" function.
---
Jeff Newmiller The . . Go Live...
DCN: Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteries O.O#. #.O#. with
/Softwar
I'd like to keep all observations for y1...
>>
>> I'm trying to plot connected time series of two variables in a
>> lattice plot:
>> xyplot(y1 + y2 ~ t, data=size, type="b")
>>
>> y2 has missing data for some of the observations and some points are
>> therefore not connected. It would make theore
On Oct 8, 2011, at 4:59 AM, Allan Sikk wrote:
Hello,
I'm trying to plot connected time series of two variables in a
lattice plot:
xyplot(y1 + y2 ~ t, data=size, type="b")
y2 has missing data for some of the observations and some points are
therefore not connected. It would make theoretica
This is what you tried? What doesn't work?
x=1:10
y=1:10
plot(x,y, type='n')
rect(4,0,6,11, col=5)
points(x,y)
___
Hello everyone,
if I have:
x=1:10
y=1:10
plot(x,y)
and I plot a rectangle
rect(4,0,6,11, col=5)
it covers the points of the graph.
Is there a way to draw the rectangle
I measured nitrate concentration and primary production (PP) biweekly for
23 months in one headwater stream. I would like to use linear regression
to determine if PP is related to nitrate concentration. My dataframe is
called "data" and consists of the vectors Rdate, PP, and nitrate. Rdate
So I've been puzzling on this for a week or so and don't see any
elegant way to do it. I'd imagine there's a nifty way to do it using
ddply() in conjunction with some melt() tools to get it all in one
data frame (or just possibly mapply() after converting the "iN"
factors to characters to make sure
It's all in the context of calling commands via the system command in R.
I just found out that you can run system("pdflatex file") and it works fine.
So R finds the command "pdflatex" but not the "del" command which I can use
to delete files.
I am wondering how I have to call this command to be a
I'd heartily suggest you read some basic R references before trying to
engage this question: a cursory glance over your code suggests it can
(and should) be entirely rewritten without for loops using
vectorization and possibly a product function like this one that I'll
provide you with:
product <-
What does this have to do with the subject line of your request or,
for that matter, with R?
Michael
On Sat, Oct 8, 2011 at 1:08 PM, syrvn wrote:
> Under MacOS I use the following script to process my latex documents:
>
> pdflatex -halt-on-error document.tex
> bibtex document.aux
> pdflatex -hal
Under MacOS I use the following script to process my latex documents:
pdflatex -halt-on-error document.tex
bibtex document.aux
pdflatex -halt-on-error document.tex
pdflatex -halt-on-error document.tex
rm *.aux *.bbl *.blg *.log *.nlo *.thm
How can I realise this under Windows 7?
--
View this m
My guess is that your intuition was correct, but you probably executed
it incorrectly. Plot (usually) calls a new plotting window; if you
want to add to a window that's already open, you need to use a
specialized command like points, lines, etc.
Try this
x <- y <- 1:10
plot(x,y,type="n"); rect(4
Hello everyone,
if I have:
x=1:10
y=1:10
plot(x,y)
and I plot a rectangle
rect(4,0,6,11, col=5)
it covers the points of the graph.
Is there a way to draw the rectangle under the points?
I was thinking the a solution could be draw an empty plot
then draw the rectangle and after the points,
but it
__
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.
Can you provide an example of where using the zoo methods directly on
the xts object fails to work? You'll note that
R> X = xts(1:5, Sys.Date()+(1:5)/2)
R> class(X)
[1] "xts" "zoo"
so the zoo methods should work automatically. This gives you access to
all sorts of options for handling NAs.
Mich
Thank you!
Sorry, I have a couple more questions:
1) How to I turn off the box shading completely? I figured out how to
lighten it up to a grey color with col=c("white"... . I looked in the
package pdf, but I didn't see anything.
2) I read about the zunit in the package pdf, but no matter what I
Hi,
I have a bunch of irregularly spaced xts time series (with a POSIX index),
and I'm trying to write a function that fillls the missing days. Using a
solution suggested by Gabor Grothendieck for zoo, I wrote the following:
# FD: Fill missing days
FD<-function(ser) {rng<-range(time(ser))
> tem
Dear Johnny,
this is a bug in the hashing-code of sets. Use:
sets_options("hash", FALSE)
lset <- cset(s, matchfun = matchfun(match.2.list))
which will work.
Thanks for pointing this out!
David
Hi,
I tried to use the sets package yesterday, which seems to be a very
powerf
Hello,
I'm trying to plot connected time series of two variables in a lattice plot:
xyplot(y1 + y2 ~ t, data=size, type="b")
y2 has missing data for some of the observations and some points are
therefore not connected. It would make theoretical sense to connect the
points - is there a way of d
May I have a question about dlmForecast() function in the package 'dlm'?
This function 'dlmForecast()' currently only deals with constant models. May
anyone suggest on how to predict using non-constant model? Thanks a lot!
Regards,
Hong Yu
On 11-10-07 5:26 PM, Carl Witthoft wrote:
Just wondering here -- I tested and found to my delight that
% round(325.4,-2)
[1] 300
gave me exactly what I would have expected (and wanted). Since it's not
explicitly mentioned in the documentation that negative 'digits' is
allowed, I just wanted t
On 11-10-07 2:32 PM, Ben qant wrote:
> Hello,
>
> I'm using the rgl package and plotting a plot with it. I'd like to
have all
> the axes values auto-hide, but I want to plot a series of characters
instead
> of the values of the measurement for 2 of the axes. So in the end I will
> have one axis
Hi,
I usually use xtable package to generate LaTeX tables. It also works for
HTML format.
A reproducible example:
-
library(gam)
library(xtable)
model <- gam(Kyphosis ~ s(Age,4)
48 matches
Mail list logo
