On 27-Oct-11 04:09:46, Jim Silverton wrote:
I have a 2 x 3 matrix called snp and I want to compute the following
probability:
choose(sum(snp[,1]), snp[1,1]) * choose(sum(snp[,2]), snp[1,2]) *
choose(sum(snp[,3]), snp[1,3])/choose(sum(snp), sum(snp[1,]))
but I keep getting Infs and NaNs.
Dear all,
i try to anderstand the AlgDesign package, so i make some tests to find with
optFederoy some plan I know.
I would like to find back a simple optimal design : 7 factors with 2
attributes each ( the beter plan is 2^(7-4) and needs 8 rows)
library(AlgDesign)
levels.design =
Hi Users,
I want to fit the maximums of a data series with a cubic spline. How do I
go about this in R.
I failed to figure out how I can use the mgcv library to do this.
Thanks
ZABLONE
__
R-help@r-project.org mailing
On 11-10-27 12:09 AM, Jim Silverton wrote:
I have a 2 x 3 matrix called snp and I want to compute the following
probability:
choose(sum(snp[,1]), snp[1,1]) * choose(sum(snp[,2]), snp[1,2]) *
choose(sum(snp[,3]), snp[1,3])/choose(sum(snp), sum(snp[1,]))
but I keep getting Infs and NaNs. Is
y1,y2 were designed as follow,from the original data file :
y1=data[,1:2]
y2=data[,3:4]
...I am a bit confused with what redesign the whole vector again and
specify indices twice actually mean?...could u point it out in the
script?
Thanks so much I really appreciate it!!
Sent from my -DROID-
Mark Newcomb newc1.mac at me.com writes:
I'm hoping to use R for spatial analysis. In working through examples in
Chapt. 4 of Applied Spatial Data
Analysis with R I've come across the following error in trying to plot lines
with the meuse data set. The
text is verbatim from the book.
This
On 10/27/2011 04:24 AM, John Sorkin wrote:
Windows XP
R 2.12
I am trying to understand how I can take a subroutine (i.e. function) I have
written in Fortran, and call the function in R. I need to know (1) where I
should store the Fortran subroutine (2) how to prepare the function for use
On 11-10-27 8:02 AM, Paul Hiemstra wrote:
On 10/27/2011 04:24 AM, John Sorkin wrote:
Windows XP
R 2.12
I am trying to understand how I can take a subroutine (i.e. function) I have
written in Fortran, and call the function in R. I need to know (1) where I
should store the Fortran subroutine
On 10/27/2011 12:13 PM, Duncan Murdoch wrote:
On 11-10-27 8:02 AM, Paul Hiemstra wrote:
On 10/27/2011 04:24 AM, John Sorkin wrote:
Windows XP
R 2.12
I am trying to understand how I can take a subroutine (i.e.
function) I have written in Fortran, and call the function in R. I
need to know
On 10/27/2011 10:46 AM, Zablone Owiti wrote:
Hi Users,
I want to fit the maximums of a data series with a cubic spline. How do I
go about this in R.
I failed to figure out how I can use the mgcv library to do this.
Thanks
ZABLONE
On Oct 27, 2011, at 12:21 AM, andrewH wrote:
Thanks Jeff! I appreciate you sharing your experience.
My data set is survey data, 13,209 records over nine years,
collected by
someone else, converted from SPSS format. It includes missing values,
identified however SPSS does so, and translated
Dear all,
I am still fairly new to R and newer still to linux (opensuse). I
want to install a package called FEAR from
http://www.clemson.edu/economics/faculty/wilson/Software/FEAR/fear-download.html;.
So, I downloaded the tar.gz file called
FEAR-linux-64bit-2.6.16.60-0.21-smp.tar.gz and it is
In your loop, you assign, for example, pv twice
pv=w$p.value # pv is scalar
pv[i]= w[2] # pv is a vector
give an example for the point
pv=1
pv[5]=2
pv
[1] 1 NA NA NA 2
This may not be what you want.
Weidong
On Thu, Oct 27, 2011 at 7:40 AM, Seb seba@gmail.com wrote:
y1,y2 were
Hi everyone
Do you know about any possibility in R to check for consecutive numbers in
vectors? That is, I do not only want to count the number of observations in
total (which can be done table(x)), but I also want to count for instance
how many times that vector contains a certain number
Hi,
I am try to fit an eGARCH model on an expanding basis using the rugarch
package. I have 6 columns of data, and I am trying refit parameters ~6000
for each column. If I run the following code, I get an error in windows on
the 2nd column (this means I am getting all the way through the first
Hi Joshua,
Thanks for the suggestion.
Sharad
--
View this message in context:
http://r.789695.n4.nabble.com/help-with-paste-tp3942818p3943753.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Hello,
I got an error if I try to install a package. For example package emg:
R CMD check emg
Everything is fine.
R CMD INSTALL emg
I got the following errors:
* installing to library 'c:/R_2_13_1/library'
* installing *source* package 'emg' ...
Error in file_path_as_absolute(dfile) :
Datei
Hello R gurus,
I have the code below for which i need help and pointers to make it run in
parallel on a dual core win7 computer with R 2.13.x, using foreach,
iterators,doMC.
library(scatterplot3d) # Loads 3D library.
library(fields)
library(MASS)
library(ROCR)
library(verification)
library(caret)
Hi,
I am using PostgreSQL 9.1 and want to run
loglinear_simple-simple.lm(x$p_overPnot,x$logQuantity,show.ci
=FALSE,conf.level=confidenceLevel)
without showing the graph, however the show.ci=FALSE is ignored.
Thank you
[[alternative HTML version deleted]]
Hi there,
What is the best way to get a time series of daily stock price observations
into R (from excel).
The time series are daily but there are spaces for w/e's and holidays etc.
So I am not sure a ts object will be suitable but I am not sure what I
should use?
What ever package you
On 27/10/2011 8:43 AM, Samir Benzerfa wrote:
Hi everyone
Do you know about any possibility in R to check for consecutive numbers in
vectors? That is, I do not only want to count the number of observations in
total (which can be done table(x)), but I also want to count for instance
how many
Hi,
Yea Linux! I think you are missing a library file (libgfortran.so.1),
you can use locate to try and find the so file, in bash type:
locate libgfortran
Now three things can happen:
- It doesn't find anything, than you need to install an additional
gfortran package. This is often called
Here is a common scenario that I and probably others often face:
- I write a script, test.R, to produce some graphs, viewing them
on-screen and tinkering until they look right
- I want to save image files (e.g., .png), so I wrap each plot in device
driver calls, e.g.,
png(file='test01.png',
On 27/10/2011 10:13 AM, Michael Friendly wrote:
Here is a common scenario that I and probably others often face:
- I write a script, test.R, to produce some graphs, viewing them
on-screen and tinkering until they look right
- I want to save image files (e.g., .png), so I wrap each plot in
On Oct 27, 2011, at 6:58 AM, Elinor Zeev wrote:
Hi,
I am using PostgreSQL 9.1 and want to run
loglinear_simple-simple.lm(x$p_overPnot,x$logQuantity,show.ci
=FALSE,conf.level=confidenceLevel)
without showing the graph, however the show.ci=FALSE is ignored.
My guess (and it is a guess
thanks for the replies so, i remove the indices pv[i] and
fc[i]..but when i run it like this:
===
fc=0
pv=0
for (i in 1:nrow(data))
{
v1= c(y1[i,1], y1[i,2])
v2= c(y2[i,1], y2[1,2])
fc=v1-v2
w=t.test(v1,v2)
pv=w$p.value
}
results =
On Oct 27, 2011, at 9:21 AM, Duncan Murdoch wrote:
On 27/10/2011 8:43 AM, Samir Benzerfa wrote:
Hi everyone
Do you know about any possibility in R to check for consecutive
numbers in
vectors? That is, I do not only want to count the number of
observations in
total (which can be done
Thanks Dimitris,
but I would like to bind the result on the dataframe, so the length
should be equal to nrow(df1).
BTW, sorry for the example, it wasn't very clear, next try:
#
options(stringsAsFactors = FALSE)
hi,
I want any suggestion for a differents packages to manipulate a time series.
thanks for help.
--
View this message in context:
http://r.789695.n4.nabble.com/package-for-time-series-tp3944825p3944825.html
Sent from the R help mailing list archive at Nabble.com.
My issue seems simple but I can't find any method to work...
I've got a matrix with 4 columns and I want to extract the values from each
row to make a new 2x2 matrix for each row. From these resultant matrices, I
want to get the Eigen Values and output them from the function as usable
numbers
Usually it's encouraged to use a csv or txt intermediary between Excel
and R: depending on what you want to do, you can either convert it
once it's in R or maybe use read.zoo() to simplify a few things. Most
people find the raw ts class hard to use and prefer a contributed
class instead: my usual
This might be an easier method:
suppose your data is in MAT
t(apply(MAT, 1, function(v) eigen(matrix(v, 2))$values))
Your problem is that return() automatically returns the output and
ceases function execution as soon as its hit.
Michael
On Thu, Oct 27, 2011 at 10:44 AM, StephenHughes
This is my first excursion into using reshape2 and I want to ensure that
the melt() function call is syntactically correct.
The unmodifed data frame is organized this way:
head(tds.anal)
site sampdate param quant
1UDS-O 2006-12-06 TDS 10800
4 STC-FS 1996-06-14 Cond 280
7
http://lmgtfy.com/?q=CRAN+time+series+analysis
Michael
On Thu, Oct 27, 2011 at 11:08 AM, rgui wa7@gmail.com wrote:
hi,
I want any suggestion for a differents packages to manipulate a time series.
thanks for help.
--
View this message in context:
Looks like you have problems with indexing. In your code, fc is a
vector of length 2, pv is a scalar(it only keeps the last value of the
loop), then you try to cbind different lengths of vectors
(row.names(y1), fc and pv). cbind only works for vectors with the same
length.
I don't know what you
Is there a package that will perform power calculations for mixed model ordinal
logistic regression? I searched and came up with nothing.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
On 27/10/2011 11:08 AM, rgui wrote:
hi,
I want any suggestion for a differents packages to manipulate a time series.
thanks for help.
See the Time Series task view on http://cran.r-project.org/web/views/.
Duncan Murdoch
__
R-help@r-project.org
Note that for factors with NA in the levels, is.na(f)[2] - TRUE
and is.na(f[2])-TRUE give different results:
f - factor(c(A,A,NA), levels=c(NA, A), exclude=NULL)
str(f)
Factor w/ 2 levels NA,A: 2 2 1
is.na(f)
[1] FALSE FALSE FALSE
is.na(f[2]) - TRUE
str(f)
Factor w/ 2
Hello,
How do I inherit two classes using the R.oo package. Below is kind of a
silly example, but I am trying to create class PerDog from classes Dog and
Person. Error at bottom. I've tried a few other ways of using extend(), but
nothing seems to get me what I want.
Example:
Hi guys,
I have 64-bit Windows 7 machine. I need to execute the latex and
pdflatex commands from R to print out tables in pdf.
I use R 2.13 64-bit, and I have GSView, GSScript and TexWorks
installed. I remember i did this once but I forgot the next step.
What shall I do?
Thanks,
--
Ş. Barış
Scott Raynaud scott.raynaud at yahoo.com writes:
Is there a package that will perform power calculations
for mixed model ordinal logistic regression? I
searched and came up with nothing.
Didn't you ask this question, and get an answer, on this list
yesterday ... ??
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Rich Shepard
Sent: Thursday, October 27, 2011 8:26 AM
To: r-help@r-project.org
Subject: [R] Syntax Check: rshape2 melt()
This is my first excursion into using reshape2
Hello
I am using glmmBUGS to fit a multilevel model. Treatments are nested in
Course are nested in Patients. The predicted variable in total EEG duration.
The predictors are:
at the observation level : Medication dose
at the Course level: Weight in KG and Age
at the Patient level: Weight in KG
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA) wrote:
What I think you want is as simple as
test.melted - melt(test)
Dan,
I see the difference this syntax makes; sampdate is the variable
associated with the value (quant) for a specific ID pair of site and param.
wanted - cast(test.melted,
Hi:
A trip to package sos reveals that simple.lm() is in the UsingR
package. Looking at the code of this function, it plots the (x, y)
pairs and the fitted least squares line without an option to suppress
the plot. Here's a slight hack of the function; it adds a new argument
plot with default
You also might be satisfied with the output of acf() or pacf() if you
are just looking at the auto-correlation of one series.
Michael
On Thu, Oct 27, 2011 at 11:18 AM, R. Michael Weylandt
michael.weyla...@gmail.com wrote:
Usually it's encouraged to use a csv or txt intermediary between Excel
Dear Paul,
Thanks for the response. I could resolve the problem. The problem, as
you had rightly pointed out, was the absence of the gfortran package
called libgfortran41. I located the rpm and installed it. This
contains libgfortran.so.1. That does it.
Thanks again.
sambit
On 27 October 2011
Dear R helpers,
I know this is a simple task, but I'm new to R and I'm still havind
difficulties with the language.
I want to create 30 vectors to be used in a simulation, each with 1 columm and
5 lines, of random numbers N(0,1).
What I tried was this:
N=150
u2-rep(1:150,0)
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Rich Shepard
Sent: Thursday, October 27, 2011 10:22 AM
To: r-help@r-project.org
Subject: Re: [R] Syntax Check: rshape2 melt()
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA)
1) Packages to be used-
For smaller datasets
use these
1. CAR Package http://cran.r-project.org/web/packages/car/index.html
2. GVLMA Package http://cran.r-project.org/web/packages/gvlma/index.html
3. ROCR Package http://rocr.bioinf.mpi-sb.mpg.de/
4. Relaimpo Package
5. DAAG
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Nordlund, Dan (DSHS/RDA)
Sent: Thursday, October 27, 2011 10:47 AM
To: r-help@r-project.org
Subject: Re: [R] Syntax Check: rshape2 melt()
-Original Message-
From:
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA) wrote:
I should apologize. I saw rshape2 and just glossed over that and read
rshape.
Dan,
Apparently the 'e' in the subject line fell off. It's 'reshape2'.
So I was using the rshape package rather than rshape2. I don't know the
relationship
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Rich Shepard
Sent: Thursday, October 27, 2011 10:56 AM
To: r-help@r-project.org
Subject: Re: [R] Syntax Check: rshape2 melt()
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA)
This book is also pretty good:
Introductory Time Series with R
Paul S.P. Cowpertwait and Andrew V. Metcalfe
http://www.springer.com/statistics/statistical+theory+and+methods/book/9
78-0-387-88697-8
-Original Message-
From: r-help-boun...@r-project.org
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA) wrote:
You should be able to get what you want using the dcast function (I think).
test.melted - melt(tds.anal)
dcast(test.melted, site + sampdate ~ param)
Dan,
I don't think the data frame melted properly. I tried the dcast() but
sampdate is
you could try this:
u2=matrix(rnorm(150), nrow=5, ncol=30); u2
Seb
On Thu, Oct 27, 2011 at 1:44 PM, Iara Faria iarag...@yahoo.com.br wrote:
Dear R helpers,
I know this is a simple task, but I'm new to R and I'm still havind
difficulties with the language.
I want to create 30 vectors to
sorry for noise
the simulated data should be like this
d=data.frame(replicate(9, rnorm(40)+10),rep(c('y','n'),20))
colnames(d)-c(LEPTIN,SAA,PTH,sEGFR,IGFBP6,MMP2,OPG,IGFBP3,PDGFAABB,group)
--
View this message in context:
if you are a beginer then try this
http://www.stat.pitt.edu/stoffer/tsa2/R_time_series_quick_fix.htm
there are lots of books and their exercises available on the authors
website. also search in google books
hello,
I got a problem with plotting large time series, since I want to store
the results in a .PDF file (I want to store several pages of plots). The
PDF files get too large to be handled ( 10MB, one was even 200MB big).
So I wonder, if there would be a possibilty to either
- reduce the
my modification gives me error
rows- c(1:nrow(mat))
scores - c()
labels -c()
itr-1000
chnksz-ceiling(itr/getDoParWorkers())
smpopt=list(chunkSize=chnksz)
foreach(icount(itr),.combine=cbind,.options.smp=smpopts)%dopar%
+ {
+ train - sample(rows, length(rows)-1)
+ label =0 ; if
Ben,
Thank you for the incredibly helpful suggestions and links. I've been
exploring each over the past few days, and for anyone else's future
reference, here's what I've found.
(1) I was able to use SANN to specify how to choose new candidate
solutions, but I wasn't able to easily use SANN for
Hi Guys,
First, English is not my native language so sorry if the question is
too difficult to understand. I can rephrase it if necessary.
I have 32-bit Windows Xp SP3, i use R version 2.12.0 (2010-10-15), and
the question is about Redundancy canonical analysis plot problem in
3D using
Thanks Michael, this worked great.
Stephen
--
View this message in context:
http://r.789695.n4.nabble.com/Trying-to-create-a-function-to-extract-values-from-a-matrix-tp3944737p3945236.html
Sent from the R help mailing list archive at Nabble.com.
__
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Rich Shepard
Sent: Thursday, October 27, 2011 11:17 AM
To: r-help@r-project.org
Subject: Re: [R] Syntax Check: rshape2 melt()
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA)
I am having trouble loading RMySQL after installing it. I have looked
at archived r-help files and it seems many people have had issues, but I
didn't find my specific issue. I also went to the webpage at Vanderbilt
(http://biostat.mc.vanderbilt.edu/wiki/Main/RMySQL) and followed the
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA) wrote:
What I think you want is as simple as
test.melted - melt(test)
Dan, et al.:
No, it's not that simple. Unfortunately. I've tried various combinations
of id.vars, measure.vars, and variable, but just cannot get the dates to
display in the
On Thu, 27 Oct 2011, Mehmet Suzen wrote:
This book is also pretty good:
Introductory Time Series with R
Paul S.P. Cowpertwait and Andrew V. Metcalfe
http://www.springer.com/statistics/statistical+theory+and+methods/book/9
78-0-387-88697-8
There's also: Time Series Analysis and its
Hello,
Why do I get preceding X. (that is a and X followed by a period) for
negative numbers and an X for positive numbers when I read a csv file? Am
I stuck with this? If so, how do I convert it to normal numbers?
dat=read.csv(file_path)
dat
[1] X0.0 X.0.240432350374
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA) wrote:
I read in the sample data that you put in an earlier post.
tds.anal - read.table('clipboard',header=TRUE)
tds.anal
site sampdate param quant
1UDS-O 2006-12-06 TDS 10800
4 STC-FS 1996-06-14 Cond 280
7UDS-O 2007-10-04
Hi,
I am converting a data.frame to xts. the data.frame is 4 columns and 1000
rows. I get a message that NROW (x) must match length(order.by)
class is data.frame, mode is list
when I run
dim(x) # I get
1000 4 #which is consistent with 1000 rows and 4 columns
NROW (x) # I get
1000
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Rich Shepard
Sent: Thursday, October 27, 2011 12:00 PM
To: r-help@r-project.org
Subject: Re: [R] Syntax Check: rshape2 melt()
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA)
Data frame is list internally so length(df) = ncol(df)
M
On Oct 27, 2011, at 2:44 PM, Muhammad Abuizzah izzah...@yahoo.com wrote:
Hi,
I am converting a data.frame to xts. the data.frame is 4 columns and 1000
rows. I get a message that NROW (x) must match length(order.by)
class is
Figured it out. Solution:
dat=read.csv(file_path, header=F)
dat
V1 V2V3 V4 V5 V6V7
V8V9V10V11V12V13 V14
V15 V16 V17 V18 V19V20
1 0 -0.2404324 0.3554681 -0.21147 1.18128
Did you load the class package before calling lda()?
Dennis
On Thu, Oct 27, 2011 at 10:14 AM, 1Rnwb sbpuro...@gmail.com wrote:
my modification gives me error
rows- c(1:nrow(mat))
scores - c()
labels -c()
itr-1000
chnksz-ceiling(itr/getDoParWorkers())
smpopt=list(chunkSize=chnksz)
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Ben quant
Sent: Thursday, October 27, 2011 12:26 PM
To: r-help@r-project.org
Subject: Re: [R] preceding X. and X
Figured it out. Solution:
dat=read.csv(file_path, header=F)
Try this, based on your small example:
tds.a - read.table(textConnection(
+ site sampdate param quant
+ 1UDS-O 2006-12-06 TDS 10800
+ 4 STC-FS 1996-06-14 Cond 280
+ 7UDS-O 2007-10-04Mg 1620
+ 9UDS-O 2007-10-04 SO4 7580
+ 19 JCM-10B 2007-06-21Ca79
+ 20
(I apologize if this ends up being posted twice, I tried to
post it earlier, and it didn't seem to 'take'.)
I am having trouble loading RMySQL after installing it. I have looked
at archived r-help files and it seems many people have had issues, but I
didn't find my specific
On Thu, Oct 27, 2011 at 1:46 PM, Jeff Laux jeffl...@gmail.com wrote:
I am having trouble loading RMySQL after installing it. I have looked at
archived r-help files and it seems many people have had issues, but I didn't
find my specific issue. I also went to the webpage at Vanderbilt
So I was using the rshape package rather than rshape2. I don't know the
relationship between those two packages and/or how they differ. I am sure
that there are others that can help you out here.
I, too, don't know how the two packages 'reshape, The Orignal' and
'reshape2, Rebooted'
I think it is what I want. The values look OK. I do get a warning. Here is
what you asked for:
dat=read.csv(file_path, header=F)
Warning message:
In read.table(file = file, header = header, sep = sep, quote = quote, :
incomplete final line found by readTableHeader on 'y:\ALL
On 28/10/11 01:13, Duncan Murdoch wrote:
SNIP
I think nowadays many Linux distros are no more ready than Windows for
this. On Linux, you'll likely need to install development libraries
and tools; on Windows, the same. The advantage of Windows is that the
tools have all been collected in one
New user here. My goal is pull daily averages from a long dataset.
I've been working with some code I got from this list from
https://stat.ethz.ch/pipermail/r-help/2009-March/191302.html
The code how I have been using it is as follows:
library(zoo)
library(chron)
yes the library(MASS) was loaded.
--
View this message in context:
http://r.789695.n4.nabble.com/help-with-parallel-processing-code-tp3944303p3945972.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
I did it using the code below thank you.
Will the fact I have used a ts object but the observations are not really
equally spaced invalidate the results?
vols=read.csv(file=C:/Documents and Settings/Hugh/My Documents/PhD/Swaption
vols.csv
+ , header=TRUE, sep=,)
x-ts(vols[,2])
x
acf(x)
--
If i understand correctly you mean to write the line as below:
foreach(icount(itr),.combine=combine,.options.smp=smpopts,.packages='MASS')%dopar%
--
View this message in context:
http://r.789695.n4.nabble.com/help-with-parallel-processing-code-tp3944303p3945954.html
Sent from the R help
On Fri, 28 Oct 2011, Hadley Wickham wrote:
See here: https://github.com/hadley/reshape
Thanks, Hadley.
Rich
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
On Thu, 27 Oct 2011, Dennis Murphy wrote:
# Define param so that all of its levels are represented:
tds.a - within(tds.a, {
param = factor(param, levels = c('TDS', 'Cond', 'Mg', 'Ca',
'Cl', 'Na', 'SO4'))
sampdate = as.Date(sampdate) } )
Dennis,
This produces an error:
On Thu, 27 Oct 2011, Rich Shepard wrote:
So I did the conversion and melting the long way:
tds.a - tds.anal
tds.a$sampdate - factor(tds.a$sampdate)
I didn't copy the melt command:
tds.anal.m - melt(tds.a)
However, ...
c.tds.anal - dcast(tds.a, site + sampdate ~ param, value_var =
I would like to combine multiple plots (one of which is a pairs plot)
in one output display:
While this works well:
layout(matrix(c(1,2,3,3,3,3), ncol=2, byrow=T))
data - matrix(rnorm(1000), ncol=10)
boxplot(data)
hist(data)
plot(data)
This doesn't:
(because the pairs function seems to
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Rich Shepard
Sent: Thursday, October 27, 2011 2:41 PM
To: r-help@r-project.org
Subject: Re: [R] Syntax Check: rshape2 melt()
On Thu, 27 Oct 2011, Dennis Murphy wrote:
#
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA) wrote:
It appears that what you did differently was to dcast the un-melted data.
In addition, you shouldn't need to turn your dates into factors. I assume
that your unaltered data is in tds.anal.
Dan,
It's been that sort of a day when that
On Thu, Oct 27, 2011 at 4:18 PM, Vinny Moriarty vwmoria...@gmail.com wrote:
New user here. My goal is pull daily averages from a long dataset.
I've been working with some code I got from this list from
https://stat.ethz.ch/pipermail/r-help/2009-March/191302.html
The code how I have been
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of Rich Shepard
Sent: Thursday, October 27, 2011 4:29 PM
To: r-help@r-project.org
Subject: Re: [R] Syntax Check: rshape2 melt()
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA)
On Thu, 27 Oct 2011, Nordlund, Dan (DSHS/RDA) wrote:
I will be leaving work shortly, and won't be able to look at this again
until tomorrow some time. However, I went back and reread the help page
on dcast. Here is part of that page.
fun.aggregate
aggregation function needed if variables do
Try another format (tiff, jpg, etc) to see how they look, what the
sizes are for different resolutions. If you have a lot of single
points, PDF files get very large because of the commands used to print
each point. If you want to keep PDF, then find some way of
aggregating the data points so
On Oct 27, 2011, at 6:27 PM, Ulrich wrote:
I would like to combine multiple plots (one of which is a pairs plot)
in one output display:
While this works well:
layout(matrix(c(1,2,3,3,3,3), ncol=2, byrow=T))
data - matrix(rnorm(1000), ncol=10)
boxplot(data)
hist(data)
plot(data)
This
I have had issues with some parallel backends not finding functions
within a namespace for packages listed in the .packages argument or
explicitly loaded in the body of the foreach loop. This has occurred
with MPI but not with multicore. I can get around this to some extent
by calling the
I am using the mgcv package to develop vis.gam plots and having trouble
figuring out how to relabel the z-axis (image attached). It is currently
labeled as linear predictor, but I would like to change it to a different
name. Currently I am using this code:
Hi All,
I'm working with some SAS code to analyze an experiment set up as follows:
66 subjects (colonies) treated with a random treatment (1-8) and measured at
three time points.
The data structure looks like:
input colony tmt y1 y2 y3;
y=y1; date=*1*; output;
I want to use
for (i in 1:time)
but I want the increments to be monthly.
For example, if I'm adding interest to a virtual bank account monthly, for a
total of 'time' in years which the user has entered.
Is there a simple way of doing this please? I'm not familiar with many R
functions yet.
1 - 100 of 107 matches
Mail list logo