You do this by adjusting the margin sizes. Can I suggest you read 'An
Introduction to R', which has a section on the layout of graphics (perhaps
the only area in which it is comprehensive).
On Mon, 16 Jun 2008, MeMooMeM wrote:
Hi R world,
I am such a newbie, with only 4-5 days of R experie
Hello everyone:
I have some quesions about the R code for linear mixed model,hope some one
can give me some hints,thanks a lot~
Say:
we have A B C three factors
(i)A is fixed ,B and C are random,and B is nested in C,so the R code for
this case would be:
case1<-lme(y~A+B+C+..,rando
Hi R friends,
I’m interested in fitting data to the following system of difference
equations:
Change in M = k*(1-M.prev)*(1-F.prev) + R*(1-M.prev)*F.pref
Change in F = k*(1-F.prev)*(1-M.prev) + S*(1-F.prev)*M.pref
Parameters:
k, R, S, are transition rates
Variables:
M = Proportion of Males
Etienne Bellemare Racine gmail.com> writes:
>
> Hi,
>
> I am working with airborne LiDAR data and R. Do you know of any package
> or project designed to achieve that goal ? Largely, if you are using R
> and .las files (or lidar data), I would like to hear about your work
> flow. Otherwise, i
On Tue, 17 Jun 2008, Chua Siang Li wrote:
Hello there. I am new to R.
I want to output the coefficients and their significance value (Pr(>|z|)
into a csv file. I managed to do it for the "coefficients" but not their
significances. Pls help. Thanks.
mylogit<- glm(response~pric
I have data set that I wish to solve with the Matrix package's sparse
matrix functionality. The speed improvements that it has achieved are
amazing, with my dense matrix solutions never taking really long
enough to time in what I've been able to time so far. However, before
I can solve my f
Hello there. I am new to R.
I want to output the coefficients and their significance value (Pr(>|z|)
into a csv file. I managed to do it for the "coefficients" but not their
significances. Pls help. Thanks.
mylogit<- glm(response~price, family=binomial(link="logit"),
Dear useRs:
Is there a way to include a raster image (e.g., .gif, .jpg, .bmp) representing
a company logo, a school logo, etc. into an R graphic? For example, it would be
nice to be able to include the logo of the school into the charts that a
student produces for her dissertation. Similarly, whil
Try max(mydensity$x), max(mydensity$y), quantile(mydensity$x), etc.
See ?density
--- On Tue, 17/6/08, Edward Wijaya <[EMAIL PROTECTED]> wrote:
> From: Edward Wijaya <[EMAIL PROTECTED]>
> Subject: [R] Accessing Max/Min Value of Density Function
> To: r-help@r-project.org
> Received: Tuesday, 17
>> mydensity <- density(x)
>> print(mydensity)
>
> x y
> Min. : -92.14 Min. :0.000e+00
> 1st Qu.: 356.66 1st Qu.:5.530e-09
> Median : 805.45 Median :4.681e-05
> Mean : 805.45 Mean :5.564e-04
> 3rd Qu.:1254.24 3rd Qu.:3.370e-04
> Max. :1703.04 Max.
Dear all,
Currently I have the following output
> mydensity <- density(x)
> print(mydensity)
x y
Min. : -92.14 Min. :0.000e+00
1st Qu.: 356.66 1st Qu.:5.530e-09
Median : 805.45 Median :4.681e-05
Mean : 805.45 Mean :5.564e-04
3rd Qu.:1254.24 3rd Qu.:3
On Jun 16, 2008, at 12:13 PM, Lord Yo wrote:
Hello everyone
I have dataset containing a monetary value (ABS) and two factors (Fct,
Group). I am able to create useful using:
bwplot(ABS~Group|Fct)
and
dotplot(ABS~Group|Fct)
Question: What do I have to do to overlay the dotplot with the
bwpl
You could try overwriting the axis with a white axis line
so it blends into the background.
plot(1:10, axes = FALSE)
axis(1)
axis(1, col = "white", tcl = 0)
On Mon, Jun 16, 2008 at 2:52 PM, Andrew Yee <[EMAIL PROTECTED]> wrote:
> I've been trying to figure out a parameter that will let you separ
Hi,
I have two figures where
I stacked together as one PNG.
Top -> scatter plot
Bottom -> density plot.
However the density curve on the bottom was cut-off
on its tip.
Is there a way I can resize the figure?
So that it can show proportional and
complete plots?
The code I have is the followin
Is this what you tried:
> x <- data.frame(a=1:10, b=1:10)
> x <- data.frame(a=1:3, b=1:3)
> do.call(rbind, list(x,x,x,x,x))
a b
1 1 1
2 2 2
3 3 3
4 1 1
5 2 2
6 3 3
7 1 1
8 2 2
9 3 3
10 1 1
11 2 2
12 3 3
13 1 1
14 2 2
15 3 3
>
On Mon, Jun 16, 2008 at 7:36 PM, Lana Schaffer <[EMAIL PRO
I am in awe of the project you developed but wish to do something simpler.
I have recently decided to give up on Microsoft Office and wish to use
OpenOffice exclusively. In the past I have used the RODBC package to
read data from excel files. Do you know how to read data directly from
OpenOffice ca
Does anyone know how to insert a jpeg image into a widget.
This is a portion of the code I have used to create the widget. I would like
to insert an image in the bottom corner.
Cheers
tt<-tktoplevel()
fontHeading2 <- tkfont.create(family="arial",size=16, slant="italic")
fontHeading3 <- tkfont
Hi
hadley wickham wrote:
> On Mon, Jun 16, 2008 at 7:52 AM, mfrumin <[EMAIL PROTECTED]> wrote:
>> Dear all,
>>
>> With normal plotting, one can size a set of points in a plot using a vector
>> argument to cex in the points() function. This works whether you are using
>> one of the standard R sym
Hello
does anyone know what algorithm is used to produce the hierarchical
clustering in the gplots package using the function heatmap.2? I think it
may be the complete linkage clustering algorithm, but I can't find a source
that seems reliable.
Thank you and sorry if I posted this in the wrong
Jim,
Is it possible to do a do.call(rbind,d.frame) with a list
like using lapply(data,func,args)?
I have not successfully executed this.
Lana
-Original Message-
From: hadley wickham [mailto:[EMAIL PROTECTED]
Sent: Friday, June 13, 2008 2:54 PM
To: jim holtman
Cc: Lana Schaffer; r-help@r-
Hi R world,
I am such a newbie, with only 4-5 days of R experience. I did a search in
forum history but couldn't find a solution to my problem... Sorry if it's
obvious:
I managed to draw a barplot (yey!) with xlabels of 'long' names (filenames,
to be particular). To make them readable, I place
On 16/06/2008 4:57 PM, Wesley Tansey wrote:
If your data is really a bunch of points, not on a grid, then you need
to use plot3d() to plot it as points, or convert it to a surface. The
interp() function in the akima package can do that. (There are lots of
other possibilities too.)
Actually
On 16/06/2008 4:49 PM, Peter Dalgaard wrote:
Andrew Yee wrote:
I've been trying to figure out a parameter that will let you separately
adjust the parameters for the axis line from the tick mark.
In the following example, I would like to suppress the axis line, but keep
the tick marks.
The so
On Mon, Jun 16, 2008 at 4:40 PM, Lord Yo <[EMAIL PROTECTED]> wrote:
>
> Hello everyone
>
> I have dataset containing a monetary value (ABS) and two factors (Fct,
> Group). I am able to create useful using:
>
> bwplot(ABS~Group|Fct)
> and
> dotplot(ABS~Group|Fct)
>
> Question: What do I have to do
On Mon, 2008-06-16 at 14:50 -0700, Taylor, Z Todd wrote:
> I generally do something like
>
> storage.mode(my.matrix) <- "double"
Ahh, that's what I was missing; I was side tracked by working out why I
couldn't write an S3 method for as.double for matrices, and missed the
proper solution in th
thanks for your correction Gavin. i read ?data.matrix and neglected to
pay attention to the last line of the description:
"Factors and ordered factors are replaced by their internal codes".
On Mon, Jun 16, 2008 at 5:30 PM, Gavin Simpson wrote:
On Mon, 2008-06-16 at 15:35 -0500, [EMAIL PR
I generally do something like
storage.mode(my.matrix) <- "double"
--Todd
--
Why do you recite at a play and play at a recital?
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Gavin Simpson
> Sent: Monday, June 16, 2008 2:31 PM
> To: [EMAIL
Dear R Helpers,
At the moment I'm working on the project to implement "optimal binning"
function. It will be primarily used as a tool for logistic regression.
something very similar to
http://www2.sas.com/proceedings/forum2008/153-2008.pdf* *but applied in
diferent problem space...*
*The prob
On Mon, Jun 16, 2008 at 4:22 PM, Alberto Monteiro wrote:
> Why does as.numeric convert matrices and arrays to vectors?
>
Because that's what it's documented to do!
" as.numeric is a generic function, but S3 methods must be written for
as.double. It is identical to as.double (and as.real). "
A
Hello everyone
I have dataset containing a monetary value (ABS) and two factors (Fct,
Group). I am able to create useful using:
bwplot(ABS~Group|Fct)
and
dotplot(ABS~Group|Fct)
Question: What do I have to do to overlay the dotplot with the bwplot (same
data set)?
I've found a couple of posts
On Mon, 2008-06-16 at 15:35 -0500, [EMAIL PROTECTED] wrote:
> hi: you can do below but i don't know if it's worth it ?
>
> newx <- data.matrix(data.frame(x))
> print(newx)
That doesn't work Mark:
> str(data.frame(x))
'data.frame': 2 obs. of 2 variables:
$ X1: Factor w/ 2 levels "1","2": 1 2
I do it a bit differently:
Once you've got myhist, you can plot a line on the existing histogram
plot with something like:
smoo<-spline(myhist$breaks[1:length(myhist$breaks)-1],myhist$counts)
lines(smoo$x,smoo$y,col='green')
You'll have to muck a bit with scale factors to make the spline curve
All,
I'm trying to adapt some code provided by Deepayan Sarkar from a previous
thread (https://stat.ethz.ch/pipermail/r-help/2005-October/081571.html) on
this topic.
## This code produces a graph w/o error bars:
xyplot(Y ~ Hr, data, groups=DRUG,
panel=function(x,y,...){
panel.x
> If your data is really a bunch of points, not on a grid, then you need
> to use plot3d() to plot it as points, or convert it to a surface. The
> interp() function in the akima package can do that. (There are lots of
> other possibilities too.)
Actually, my data is on a grid. The values ar
Thanks for tracking this down.
Andrew
On Mon, Jun 16, 2008 at 4:49 PM, Peter Dalgaard <[EMAIL PROTECTED]>
wrote:
> Andrew Yee wrote:
>
>> I've been trying to figure out a parameter that will let you separately
>> adjust the parameters for the axis line from the tick mark.
>>
>> In the following
Andrew Yee wrote:
I've been trying to figure out a parameter that will let you separately
adjust the parameters for the axis line from the tick mark.
In the following example, I would like to suppress the axis line, but keep
the tick marks.
The source code is the ultimate reference. We have i
On 6/16/2008 4:40 PM, Duncan Murdoch wrote:
On 6/16/2008 4:25 PM, Floyd poole wrote:
I am new to R and what to find out which package would be best to create a
surface plot, 2d or 3d if possible. I have a matrix (depth.dat) which has
over 15k depth measurements from an ultrasonic corrosion mapp
On 6/16/2008 4:25 PM, Floyd poole wrote:
I am new to R and what to find out which package would be best to create a
surface plot, 2d or 3d if possible. I have a matrix (depth.dat) which has
over 15k depth measurements from an ultrasonic corrosion mapping
application. These depth measurements onl
hi: you can do below but i don't know if it's worth it ?
newx <- data.matrix(data.frame(x))
print(newx)
On Mon, Jun 16, 2008 at 4:22 PM, Alberto Monteiro wrote:
Why does as.numeric convert matrices and arrays to vectors?
as.numeric(matrix(c("1", "2", "3", "4"), 2, 2))
[1] 1 2 3 4
I could
On 6/16/2008 1:32 PM, Wesley Tansey wrote:
Hi all!
I'm very new to R, and I'm having trouble figuring out how to go from a file
of points that I have to a 3d surface plot of the data. I typically have
something like this:
X Y Z
0.0050.02334.45
0.0035 0.63
I am new to R and what to find out which package would be best to create a
surface plot, 2d or 3d if possible. I have a matrix (depth.dat) which has
over 15k depth measurements from an ultrasonic corrosion mapping
application. These depth measurements only have one axis (x) but I have used
Octave
I originally sent this to Doug Bates but have received no reply yet so I
thought I would expand to a wider source.
I've been trying to estimate linear mixed effect models in lmer() from the
lme4 package using the weights option. The help and code for lmer()
suggest to me that this is implement
Why does as.numeric convert matrices and arrays to vectors?
as.numeric(matrix(c("1", "2", "3", "4"), 2, 2))
[1] 1 2 3 4
I could only figure out ugly ways to bypass this, like:
x <- matrix(c("1", "2", "3", "4"), 2, 2)
array(as.numeric(x), dim = dim(x), dimnames = dimnames(x))
Alberto Monteiro
_
Jin Wang wrote:
I tried to compare if cch() and coxph() can generate same result for
same case cohort data
Use the standard data in cch(): nwtco
Since in cch contains the cohort size=4028, while ccoh.data size =1154
after selection, but coxph does not contain info of cohort size=4028.
The roug
On Mon, 16-Jun-2008 at 11:35AM -0700, hanen wrote:
|>
|>
|> in order to minimize the police of the title of my graph i tried to use:
|>
|> >title("nombre de fleurs données journellement par 6 cereus peruvianus du
|> 1/07/02 au 31/09/02",font=0)
|>
|> but this message appears:
|>
|> Erreur dan
Hi all
I would be grateful you can help me with my problem.
I try to run an optimization code . in one line I have runif in order to
sample the PDF. I get this error while i run it.
Error in runif(1, f$d[[n.of.u.vars + n.of.o.vars + j]][[2]][1],
f$d[[n.of.u.vars + :
invalid argumen
Perhaps something like this:
x <- sample(10:50, 10, rep=TRUE)
cols <- as.character(cut(x,
breaks = c(0, 10, 20, 30,max(x)),
labels = c("green", "blue", "yellow",
"red")))
barplot(x, col = cols)
On Mon, Jun 16, 2008 at 4:15 PM,
- begin included message -
I tried to compare if cch() and coxph() can generate same result for
same case cohort data
Use the standard data in cch(): nwtco
Since in cch contains the cohort size=4028, while ccoh.data size =1154
after selection, but coxph does not contain info of
Hello,
I want to use R to produce nice meteograms...
For this I would need to color-code my barplot according to the value plotted
say if value x > 30 ... the bar is red ; if x between 20 and 10, the bar
is blue etc...
any ideas how to proceed ?
Thank You a lot,
maria
Marc Bernard yahoo.fr> writes:
>
> Is there a function to calculate the prediction intervals for random
effects in non-linear mixed
> models? I found a way to do it for linear mixed models but not for
non-linear mixed one.
Please do not send HTML mail to the list.
A discussion and a caveat by
I've been trying to figure out a parameter that will let you separately
adjust the parameters for the axis line from the tick mark.
In the following example, I would like to suppress the axis line, but keep
the tick marks.
Thanks,
Andrew
foo <- data.frame(x=1:3, y=4:6)
plot(foo$x, foo$y, type="
Try with mtext:
plot(rnorm(100))
mtext(text = "Left", side = 3, adj = 0)
On Mon, Jun 16, 2008 at 2:39 PM, Carlos Gershenson <[EMAIL PROTECTED]>
wrote:
> Hi,
>
> I am trying to insert a letter in a plot corner outside the plotting area.
> Thus, "legend" and "text" don't seem to work. "title" does
I have a problem where I need to label the vertical axes of a Boxplot
with related, but different quantities (flow & height), which have a
known relationship. Primarily I want to plot the variable as a flow on
the left axis and on the opposing right axis, show the corresponding
height. Is it po
in order to minimize the police of the title of my graph i tried to use:
>title("nombre de fleurs données journellement par 6 cereus peruvianus du
1/07/02 au 31/09/02",font=0)
but this message appears:
Erreur dans title("nombre de fleurs données journellement par 6 cereus
peruvianus de 1/07/02
Hello,
I am attempting to assign significance levels to a UPGMA cluster analysis as
part my doctoral research. The pvclust function works well but doesn't include
the similarity index I need (morisita's) as an option for computing a distance
matrix. Morisita's is available in vegdist in the
Hello list,
I am trying to test a model but for the beginning I want to do this by using
simulated dataset. The model is
Y_t = X_t %*% beta + e
Where Y : (Nx1); X: (Nxp); beta: (0.6,0.3,0.1); e-uncorrelated normally
distributed variates for each t. and later I want to use to use this datas
Hi,
I am trying to insert a letter in a plot corner outside the plotting
area. Thus, "legend" and "text" don't seem to work. "title" does the
trick, but I cannot find a way of moving it from the center to the
left corner... I already tried with a few parameters from par, but
title does no
Thanks your suggestion works -:-)
From: milton ruser [mailto:[EMAIL PROTECTED]
Sent: Thursday, June 12, 2008 2:42 PM
To: Manisha Brahmachary
Cc: r-help@r-project.org
Subject: Re: [R] save workspace while running R on a cluster
Hi Manisha,
How about you i
Dear UseRs,
I wrote the following function to use MLE.
-
mlog <- function(theta, nx = 1, nz = 1, dt){
beta <- matrix(theta[1:(nx+1)], ncol = 1)
delta <- matrix(theta[(nx+2):(nx+nz+1)], ncol = 1)
sigma2 <- theta[nx+nz+2]
gamma <- theta[nx+nz+3]
Hello everyone
I have dataset containing a monetary value (ABS) and two factors (Fct,
Group). I am able to create useful using:
bwplot(ABS~Group|Fct)
and
dotplot(ABS~Group|Fct)
Question: What do I have to do to overlay the dotplot with the bwplot (same
data set)?
I've found a couple of posts
Hi Hadley,
thanks for your tip! It worked fine (though I had to use "tile" instead
of "point"). I see your point to handle legends the way you do, which is
a very convenient feature of ggplot. It's great for points, lines, bars
etc., but in the particular case of error bars, I think it would be
re
OK, it seems like I do not succeed in expressing what I do, or want to do.
So, I give you the example that bring me to this kind of analysis. I wrote
the paper "Chromosomal periodicity of evolutionary conserved gene pairs"
(which you can download at http://www.pnas.org/cgi/reprint/104/25/10559). In
I tried to compare if cch() and coxph() can generate same result for
same case cohort data
Use the standard data in cch(): nwtco
Since in cch contains the cohort size=4028, while ccoh.data size =1154
after selection, but coxph does not contain info of cohort size=4028.
The rough estimate between
Try scatterplot3d function in the package with same name.
On Mon, Jun 16, 2008 at 2:32 PM, Wesley Tansey <[EMAIL PROTECTED]> wrote:
> Hi all!
>
>
>
> I'm very new to R, and I'm having trouble figuring out how to go from a
> file
> of points that I have to a 3d surface plot of the data. I typicall
Hi all!
I'm very new to R, and I'm having trouble figuring out how to go from a file
of points that I have to a 3d surface plot of the data. I typically have
something like this:
X Y Z
0.0050.02334.45
0.0035 0.63 28.48
.
I've tried looking at the persp
Hi,
did you have a look at the manual "An Introduction to R"?
(http://cran.r-project.org/doc/manuals/R-intro.html).
Chapter "Statistical Models in R" should answer (most of) your questions
concerning linear regression.
kayj wrote:
Also each x is a categorical variable with 4 categories so I
John is correct & the error message you got was clear: the error SSP
matrix is of deficient rank.
Possible solutions:
- Regroup the factor levels to fewer and/or
- Combine some of the species to broader classes so you have fewer than
17 response variables
- Collect more samples
hope this help
Dear all,
This is a relist of my previous question. I noticed that some
characters were not displayed in the previous version.
Is there a function to calculate the prediction intervals for random
effects in non-linear mixed models? I found a way to do it for linear mixed
models but not
> thanks for your tip! It worked fine (though I had to use "tile" instead
> of "point"). I see your point to handle legends the way you do, which is
> a very convenient feature of ggplot. It's great for points, lines, bars
> etc., but in the particular case of error bars, I think it would be
> reas
> Thanks for getting back to me so quickly. This mostly works, in that:
>
> - When I ran it just as you sent it, I get this error:
> Error in inherits(x, "factor") : object "fill" not found
Oh yeah, that's a change in the development version.
> - but when I took out "fill = fill" from your s
I think that you have to define some type of criteria on what you
expect the performance to be. There is a cost/benefit tradeoff in
having multiple servers with a load balancer and this all depends on
how much you want to spend on a solution. You might want to prototype
it in R to at least get an
Hi list,
I'm struggling with plotting a 3d data set.
I have a data frame with x,y and z values. They occur in no particular order
(ie x and y are not ordered).
However, I want to visualize the data using a contour plot.
Now, the x,y and z data going into contour() are supposed to be ordered.
Dear UseRs,
I wrote the following function to estimate parameters using MLE.
mlog <- function(theta, nx = 1, nz = 1, dt){
beta <- matrix(theta[1:(nx+1)], ncol = 1)
delta <- matrix(theta[(nx+2):(nx+nz+1)], ncol = 1)
sigma2 <- theta[nx+nz+2]
gamm
The logic behind the round to even rule is that we are trying to represent an
underlying continuous value and if x comes from a truly continuous
distribution, then the probability that x==2.5 is 0 and the 2.5 was probably
already rounded once from any values between 2.45 and 2.54...,
On Mon, Jun 16, 2008 at 3:55 PM, hadley wickham <[EMAIL PROTECTED]> wrote:
> On Mon, Jun 16, 2008 at 7:52 AM, mfrumin <[EMAIL PROTECTED]> wrote:
> >
> > Dear all,
> >
> > With normal plotting, one can size a set of points in a plot using a
> vector
> > argument to cex in the points() function. Th
Everything was read in the same way, and str(junk1) confirms that they are
the same structure. This is very strange.
## original data:
> str(junk1)
'data.frame': 96 obs. of 3 variables:
$ Hour: int 0 3 5 0 3 5 0 3 5 0 ...
$ Drug: Factor w/ 2 levels "D","P": 2 2 2 1 1 1 2 2 2 1 ...
$ Aldo:
Greg Snow wrote:
> It looks like your original data may be tab seperated, if that is the case
> then just use read.delim or use sep='\t' in read.table or scan.
>
I think that was only half the problem. If you do that, you end up with
one column containing both City and State, comma-separated. P
Hi David,
If the delimier is "tab" try this.
my.df<-read.table("my_file.txt", head=T, sep="\t")
Cheers,
Miltinho
Brazil
>
>
>
> > -Original Message-
> > From: [EMAIL PROTECTED]
> > [mailto:[EMAIL PROTECTED] On Behalf Of David Arnold
> > Sent: Friday, June 13, 2008 5:15 PM
> > To: r-hel
To get some sort of frequency which in your case seem to be cycles per
distance? Is a valid use of a fourier transform as long as it is a distance
that is measured in a way that would be analogous to a time series- In
other words if the distance proceeds from an origin in one direction-
geophysic
It looks like your original data may be tab seperated, if that is the case then
just use read.delim or use sep='\t' in read.table or scan.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
> -Original Message-
> From: [EMA
All,
I'm re-running some analysis that has been augmented with additional data.
When I use the exact same code for the augmented data, the behavior of the
aggregate function is very strange, viz., one of the resulting variables is
now coded as a factor while it was coded as numeric for the origi
You can embed plots using the subplot function in the TeachingDemos package.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Alex Rey
On Mon, 16 Jun 2008, Peter Dalgaard wrote:
Prof Brian Ripley wrote:
On Mon, 16 Jun 2008, Peter Dalgaard wrote:
Millo Giovanni wrote:
The only thing I wish I were able to do now is to have the graphics
windows defaults changed to a size fitting the small 9'' screen, as now
I have to reduce i
On Mon, Jun 16, 2008 at 7:52 AM, mfrumin <[EMAIL PROTECTED]> wrote:
>
> Dear all,
>
> With normal plotting, one can size a set of points in a plot using a vector
> argument to cex in the points() function. This works whether you are using
> one of the standard R symbols (i.e. 19+) or some ascii sy
Hi there,
The subset function is available on your Base R. So you need not to load any
package.
May be it is not working because I suppose that your "rownames" is a
sequential one (1...4000). Case it is not true, the command will not work.
Anyway, I think that those line that are working find is
Hi Carsten,
In my eagerness to get every component of the legend working
automatically, I've completed neglected a way to turn bits off of you
don't want them (the philosophy is that every geom that uses an
aesthetic should appear in the legend in some way). The best thing I
can suggest at the mo
Hi All,
I have a data with 15 columns. the first 10 column are dependent(y’s)
variables and the following 5 an columns are the independent variables(x’s)
( MY DATA HAS A HEADER WITH THE NAME OF THE VAIABLES). I need to apply
linear regression y=a+bx where I need x to be one of the independent
va
Dear UseRs,
I wrote the following function to use MLE.
-
mlog <- function(theta, nx = 1, nz = 1, dt){
beta <- matrix(theta[1:(nx+1)], ncol = 1)
delta <- matrix(theta[(nx+2):(nx+nz+1)], ncol = 1)
sigma2 <- theta[nx+nz+2]
gamma <- theta[nx+nz+3]
Thanks for all the replies.
My data frame is actually 4000 rows and 10 columns.
This line,
df <- subset( df, !( rownames(df) %in% 1:25 ) ) Did not work for me. Is
there a certain library I need to load for this to work (I'm assuming there
isn't)? I didn't get an error message, warning, or anyt
Perhaps I'm applying spec.pgram wrong as you said. I will explain what I
want, so you can tell me why I'm wrong and perhaps what I have to do to do
it well.
I have some points in a 1-D space and I want to know if they are spaced at a
certain periodic distance. So, I computed all the distances betwe
Try using Sweave
On Mon, Jun 16, 2008 at 10:37 AM, David Keegan <[EMAIL PROTECTED]>
wrote:
>
> Hi,
>
> I have a non-interactive R script that currently produces
> various graphs in png or pdf format. I need to program the
> script to combine the graphs with various pages of textual
> information,
Hi,
I have a non-interactive R script that currently produces
various graphs in png or pdf format. I need to program the
script to combine the graphs with various pages of textual
information, including some in tabular format, into an
output report in pdf or html format.
What is the recommended
On Mon, 2008-06-16 at 14:38 +0200, Luis San José García wrote:
> Hi,
>
> Somebody knows how to make contrasts if i'm using the function adonis?
What is wrong with the defaults? From ?adonis:
Usage:
adonis(formula, data, permutations = 5, method = "bray",
strata = NULL, contr.un
On Jun 16, 2008, at 6:00 AM, [EMAIL PROTECTED] wrote:
> From: eugen pircalabelu <[EMAIL PROTECTED]>
> Date: June 15, 2008 7:16:09 PM EDT
> To: R-help <[EMAIL PROTECTED]>
> Subject: [R] R vs SAS and HLM on multilevel analysis- basic question
>
>
> Hi R users!
>
> I am trying to learn some multilev
Look here:
http://msenux.redwoods.edu/mathdept/R/CentralLimit.php
D.
On Jun 15, 2008, at 11:25 PM, Gundala Viswanath wrote:
Hi,
Currently I have a density plot generated with this
snippet.
Is there a way I can add a line curve on top of it?
I mean in one figure
__BEGIN__
myhist <- hist
Hi ,
I don't know whether it is a Mac specific problem. I hope it does
not bother you.
When I start R each time I got the following warning , I reinstall for
many times the problem still exists.
R(4053,0xa0619fa0) malloc: *** error for object 0x15630350: double
free
*** set a breakpo
We may want to move this discussion to the R-SIG-Mixed-Models list,
which I have cc:'d on this reply.
On Sun, Jun 15, 2008 at 6:16 PM, eugen pircalabelu
<[EMAIL PROTECTED]> wrote:
> Hi R users!
> I am trying to learn some multilevel analysis, but unfortunately i am now
> very confused. The reaso
Hi,
Somebody knows how to make contrasts if i'm using the function adonis?
Thanks.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provi
Dear all,
With normal plotting, one can size a set of points in a plot using a vector
argument to cex in the points() function. This works whether you are using
one of the standard R symbols (i.e. 19+) or some ascii symbol, such as '/'
eg:
plot(1:10, 1:10, type='n');
points(1:10, 1:10, cex = 1
Prof Brian Ripley wrote:
> On Mon, 16 Jun 2008, Peter Dalgaard wrote:
>
>> Millo Giovanni wrote:
>>>
>>> The only thing I wish I were able to do now is to have the graphics
>>> windows defaults changed to a size fitting the small 9'' screen, as now
>>> I have to reduce it and move it to the right b
1 - 100 of 123 matches
Mail list logo
