On Feb 10, 2013, at 11:36 PM, zhu yao wrote:
Dear R-users
I want to generate a growth chart like this:
http://php.med.unsw.edu.au/embryology/images/9/9a/WHO-XX_length_birth_to_2_years.png
can you give some suggestions about the functions I should use?
Looking first at:
Thank you very much for the reference about methodology.
I will try the GAMLESS package.
*Yao Zhu*
*Department of Urology
Fudan University Shanghai Cancer Center
Shanghai, China*
2013/2/11 David Winsemius dwinsem...@comcast.net
On Feb 10, 2013, at 11:36 PM, zhu yao wrote:
Dear R-users
I
Dear list
members,
I would
like to create two series of plotted mean values and error bars, yet with
different locations along the x-axis.
Plotting of
first series using plotCI with the standard arguments goes without any problem.
However I do not succeed to add the second series in the same
Hello,
I am drawing a boxplot using bwplot in lattice.
I would like to learn how to draw a thicker median line instead of the
default setting.
Thank you.
Code
bwplot(,
pch=rep(|,2))
Elaine
[[alternative HTML version deleted]]
__
On 10.02.2013 22:43, David Winsemius wrote:
On Feb 10, 2013, at 12:28 PM, Uwe Ligges wrote:
On 10.02.2013 21:18, Paul Ossenbruggen wrote:
I installed XCode, XQuartz and OpenBUGS on my Mac using WINE.
When I tried to test my installation I received the following error.
If
Hello,
excuse me if this is trivial. I have some survey data with the following
design
MyDesign-(id=ident, weights=~fwgt, fcp=~nval1+nval2, data=hh_data1)
I would like to run simple frequency tables such as those provided by the
freq() function of rgrs to check the impact of the sampling design
On 08/02/2013, e-letter inp...@gmail.com wrote:
Readers,
For a 12 hour time stamp:
testtime-(2013-01-01 01:00:01 PM)
testtime
[1] 2013-01-01 01:00:01 PM
testtime24hour-strftime(testtime,'%H:%M:%S')
testtime24hour
[1] 01:00:01
Have realised that the object is of type 'character'. When
Le lundi 11 février 2013 à 10:37 +0100, Virgile Capo-Chichi a écrit :
Hello,
excuse me if this is trivial. I have some survey data with the following
design
MyDesign-(id=ident, weights=~fwgt, fcp=~nval1+nval2, data=hh_data1)
I would like to run simple frequency tables such as those
Hi Rosalinde,
If you set the argument err equal to x, you'll get horizontal error bars.
Hope this helps,
José
José Iparraguirre
Chief Economist
Age UK
T 020 303 31482
E jose.iparragui...@ageuk.org.uk
Twitter @jose.iparraguirre@ageuk
Tavis House, 1- 6 Tavistock Square
London, WC1H 9NB
Hello everyone,
I was going through the optimization techniques implemented in R. I was
able to find most of the techniques that have been ported in R leaving a
few. Can anyone tell me if a package on quadratic Active set technique has
been developed in R as this was one of the requirements of my
On 02/11/2013 08:01 PM, linde wrote:
Dear list
members,
I would
like to create two series of plotted mean values and error bars, yet with
different locations along the x-axis.
Plotting of
first series using plotCI with the standard arguments goes without any problem.
However I do not succeed
Hello,
Em 11-02-2013 01:26, Pete Brecknock escreveu:
malaka wrote
Hi,
I want to assign the ar1 , ma 1 and the intercept estimated by the
following code to three variables a, b and c respectively.
Can anyone help me with this please?
code:
a0 = 0.05; a1 = 0.1; b1 = 0.85
nu = rnorm(2500)
epsi
On 02/09/2013 11:52 PM, Eike Marie Thaysen wrote:
Hi,
I have run into this problem a couple of times now and hope you can help!
If I want to plot mulitiple series with differing x-axis values (however, all
in the same range) in 1 plot with one common axis R obstruses the plots.
E.g. if I
Dear R list,
I have a problem in assigning values to existing data frames in R. I have a
vector x containing the names of the data frames, which I create to store
the results for each variable (a1,a2,a3) obtained in time series moving
regressions. Thus, say,
x=c('a1','a2','a3')
Moreover, b is a
On 13-02-11 6:58 AM, ivan wrote:
Dear R list,
I have a problem in assigning values to existing data frames in R. I have a
vector x containing the names of the data frames, which I create to store
the results for each variable (a1,a2,a3) obtained in time series moving
regressions. Thus, say,
Readers,
Have since tried to plot converted 24 hour data:
testtimedataset
V1
1 13:01:41
2 13:02:10
3 13:02:38
4 13:03:05
5 13:03259
testdata-seq(1:5)
plot(testdata~testtimedataset)
Error in function (formula, data = NULL, subset = NULL, na.action = na.fail, :
invalid type (list)
Hello,
You're trying to plot the df, not the column. Try
plot(testdata ~ testtimedataset$V1)
Hope this helps,
Rui Barradas
Em 11-02-2013 13:30, e-letter escreveu:
Readers,
Have since tried to plot converted 24 hour data:
testtimedataset
V1
1 13:01:41
2 13:02:10
3 13:02:38
4
Hi
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of e-letter
Sent: Monday, February 11, 2013 2:31 PM
To: r-help@r-project.org
Subject: Re: [R] convert 12 time stamp to 24 hour
Readers,
Have since tried to plot converted
Dear Uwe,
thanks for the response.
I knew my matrix was almost singular, but is there any way to implement the
algorithm is solve(base) with a spam-matrix?
As workaround, I might add a ridge penalty:
ApP - A+10^-1*diag.spam(nrow(A))
solve(ApP)
but this would completely jeopardize
On 11/02/2013, PIKAL Petr petr.pi...@precheza.cz wrote:
str(testtimedataset)
With my real data set, revealed that after conversion, the time series
data consists of 'chr' which I guess means character type data,
reporting the following error:
Error in plot.window(...) : need finite 'xlim'
Berend,
I had pointed out this solution to Axel already, but he aparently wants to
solve more general problems like this, so he wanted a solution from constrained
optimimzers. But, if the more general problems have the similarly discrete set
of candidate parameter values, using constrained
On 11.02.2013 15:00, Camarda, Carlo Giovanni wrote:
Dear Uwe,
thanks for the response.
I knew my matrix was almost singular, but is there any way to implement the
algorithm is solve(base) with a spam-matrix?
Not sure, I think you want another threshold for the error. But I do not
know if
Thank you very much, indeed I have mixing up the x and y values.
Best, Rosalinde
Â
De : Jim Lemon j...@bitwrit.com.au
Cc : r-help@r-project.org r-help@r-project.org
Envoyé le : Lundi 11 février 2013 11h41
Objet : Re: [R] How to plot doubles series
James Jong ribonucleico at gmail.com writes:
I have a binary classification problem where the fraction of positives is
very low, e.g. 20 positives in 10,000 examples (0.2%)
What is an appropriate cross validation scheme for training a classifier
with very few positives?
[snip]
Hi
-Original Message-
From: e-letter [mailto:inp...@gmail.com]
Sent: Monday, February 11, 2013 3:11 PM
To: PIKAL Petr
Cc: r-help@r-project.org
Subject: Re: [R] convert 12 time stamp to 24 hour
On 11/02/2013, PIKAL Petr petr.pi...@precheza.cz wrote:
str(testtimedataset)
I once knew the answer to this, but can't find where I found it.
Say I have a vignette that I don't want R to automatically rebuild, perhaps it uses some
data I can't share or takes too long to run or whatever. Say it is
mypackage/vignettes/charlie.Rnw
It can create the pdf from an R
On 11-02-2013, at 15:00, Camarda, Carlo Giovanni cama...@demogr.mpg.de
wrote:
Dear Uwe,
thanks for the response.
I knew my matrix was almost singular, but is there any way to implement the
algorithm is solve(base) with a spam-matrix?
As workaround, I might add a ridge penalty:
On 11.02.2013 15:44, Terry Therneau wrote:
I once knew the answer to this, but can't find where I found it.
Say I have a vignette that I don't want R to automatically rebuild,
perhaps it uses some data I can't share or takes too long to run or
whatever. Say it is
On Feb 11, 2013, at 1:23 AM, Elaine Kuo wrote:
Hello,
I am drawing a boxplot using bwplot in lattice.
I would like to learn how to draw a thicker median line instead of the
default setting.
Thank you.
Code
bwplot(,
pch=rep(|,2))
At the moment panel.bwplot has code at the end
You can use the functions in lattice to compute the shingle and display the
results, but they do not seem to be stored:
set.seed(42)
vals - runif(100)*100
ints - co.intervals(vals, number=5, overlap=.25)
cts - shingle(vals, ints)
cts contains has the original data and the range of each
Dear Elisha,
It's amazing how informative posting a reproducible example can be:
The problem is that your model has aliased coefficients, and effect() can't
handle that. Coincidentally, the development version of the effects package on
R-Forge *can* handle models with aliased coefficients, so
To get the counts in a variable:
library(lattice)
set.seed(42)
vals - runif(100)*100
ints - co.intervals(vals, number=5, overlap=.25)
shgl - shingle(vals, ints)
summary(shgl)
Intervals:
min max count
1 0.02362834 25.5431425
2 20.76563596 47.4999725
3 39.84827982
Dear R users,[IF THE FORMAT OF MY EMAIL IS NOT CLEAR, I HAVE ATTACHED A TEXT
FILE FOR A CLEAR VIEW]
I would like to use the R output file in Fortran. my file Is exactly in the
following format.
ELISA/BOTTO wATER INN
FROM 1900 11 1 TO 1996 12 31
1901.11. 1 447.000
1901.11. 2
Dear all,
I am running an Ecological Niche Factor Analysis in R but I am stuck with a
problem. As soon as I ask the software to compute the ENFA, I get the
following error message:
Error in eigen(Se) : infinite or missing values in 'x'
Does anyone know what could be wrong. I do not have missing
Hi,
I would like to visualize a tree extracted from a random forest using
getTree {randomForest}.
I'm wondering if there is any way to do it directly or to convert my tree
to any other class of tree repesentation, which then can be plotted?
Regards,
Umang Rathi
Sorry about that. Attached is the relevant data and the following is the
code and error.
Help-read.csv(R-help-effect.csv, sep=,, header=T)
Help$WeeksOut-as.factor(Help$WeeksOut)
Help$Run-as.factor(Help$Run)
GLM-glm(NumberSpecies~Treatment*WeeksOut+Run+Tree., family=poisson
(link=log), data=Help)
Brian thanks for answer me,
The idea is to make a hierarchical cluster of 138037 pixels of 1 km^2
from a study area of colombian Andes as a biogeographical regionalization
. I have distributions models for 89 species so i have a matrix with the
pixels in the rows and is full with
Dear All,
I have been using RExcel for some time and find two features especially useful:
i) The ability to select part of a worksheet and put it into R
ii) The way of viewing the code, the output and the alerts in separate panes.
The rest of R Excel I don't really use.
The major benefit for me
Hello again,
My query may look quite generic, however at this point of time I just
want explain my problem. I am hopeful that somebody can help me.
Let say I have a script file 'Script.R'. Now I want R to run this
script file each day at a specific time, and store the result (the
calculation for
man cron
or something more robust:
http://jenkins-ci.org/
-Whit
On Mon, Feb 11, 2013 at 1:51 PM, Christofer Bogaso
bogaso.christo...@gmail.com wrote:
Hello again,
My query may look quite generic, however at this point of time I just
want explain my problem. I am hopeful that somebody can
Hi Elisa,
It seems also that there are spaces between the second '.' and the number
(1901.11. 1).
You could do something like:
Lines1-readLines(textConnection(1901.11. 1 447.000
1901.11. 2 445.000
1901.11. 3 445.000
1924.11. 4 445.000
1924.11. 5 449.000
1924.11. 6 442.000
I've been searching for information on using deques in R for the past couple of
days. What little I've found has pointed me towards either the Rcpp package or
the cpp project on R-forge (https://r-forge.r-project.org/projects/cpp/), but
I've had trouble finding working examples. Are these the
Depends on system you are working.
Mac : Automator
Windows : save the script in the folder where python is installed and then run
it through command prompt.
Linux and unix specific : cron jobs.
As Mac is unix underneath so cron will work there also but never tried it.
Best Regards,
Dear R-helpers,
The vectors in my biplot are completely obscured by the ~1400 labels R is
printing on my biplot. I don't really care about the labels. How can I make
the biplot without the annoying labels?
See attached, if that helps you see my problem.
Many thanks,
Mark Na
?biplot
set.seed(42)
pcout - princomp(matrix(rnorm(1000), 100, 10))
biplot(pcout)
biplot(pcout, xlabs=rep(., dim(pcout$scores)[1])) # small symbol
biplot(pcout, xlabs=rep(, dim(pcout$scores)[1])) # no symbol
--
David L Carlson
Associate Professor of
cron does work just fine on the Mac, and will work the same way as in
Linux, except (almost certainly) for the path to R itself.
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 2/11/13 11:16 AM, Bhupendrasinh Thakre
This doesn't answer your question, but you may find some of my functions
useful:
https://github.com/geneorama/geneorama/blob/master/R/clipper.R
https://github.com/geneorama/geneorama/blob/master/R/clipped.R
You can quickly get stuff from Excel using the clipboard using clipped().
I also added
You say, Š use the R output file in Fortran Š
I guess that means that R is writing an output file which will then be
used as an input file for a fortran program.
In that case, you need to go back to how R is writing the output file,
find out why it is writing blank lines, and correct it. As far
On Mon, Feb 11, 2013 at 5:20 PM, Michael Holt mh...@terapeak.com wrote:
I've been searching for information on using deques in R for the past couple
of days. What little I've found has pointed me towards either the Rcpp
package or the cpp project on R-forge
On Mon, Feb 11, 2013 at 9:42 PM, Barry Rowlingson
b.rowling...@lancaster.ac.uk wrote:
but I can't find it on CRAN now...
http://cran.r-project.org/web/packages/Containers/
so you could hook the code out of the archives and find out if it still
works...
I just looked at the DESCRIPTION
Dear Macqueen,Although this question might be out of context, but just for
clarification of my own, i got fortran codes from my colleague who is using MAC
OS. Whereas, i am using XP and windows 7 with CVF 6.6. The codes run
efficiently on his system but for me it prompts up with an error that
Hi
I had a similar problem - see
http://tolstoy.newcastle.edu.au/R/e9/help/10/03/7129.html
DF -
data.frame(site = factor(rep(1:5, each = 20)),
height = rnorm(100))
bwplot(height~ site,DF,
pch = |,
panel = function(x, y, ..., horizontal) {
Thanks David. ?biplot is not exactly verbose on the issue, but in hindsight
I can see how I might have figured that out on my own. In any case, I
appreciate the time you took to show me how to do that.
Best wishes, Mark Na
On Mon, Feb 11, 2013 at 2:56 PM, David L Carlson dcarl...@tamu.edu
Hello,
I am using lattice to draw boxplots.
I would like to make the whisker look shorter than the box width.
Please kindly advise how to modify the code.
Thank you.
Elaine
code
bwplot(x,y,
pch=rep(|,2),
box.umbrella = list(col = black, lty = 1))
[[alternative HTML version deleted]]
Hi there,
I had an error when using RSQLite. My code and the error message are
the following:
dbGetQuery(sql.dispersion,
+create table monthly_data_temp as
+ select a.*, b.industry
+ from monthly_data as a left join
+siccd_industry_ff49
Dear R help;
I'll preface this by saying that the example I've provided below is pretty
long, turgid, and otherwise a deep dive into a series of functions I wrote
for a simulation study. It is, however, reproducible and self-contained.
I'm trying to do my first simulation study that's quite big,
Hi,
I've just read your intro, not the full code. But one way first divide your
input data into parts and make a list where each element of the list is one
part of your input data. Rest of code put in a function, which will work on
a part of your data. Build the function such that it takes list
Hi Ben,
I appreciate that the example is reproducible, and I understand why
you shared the entire project. At the same time, that is an awful lot
of code for volunteers to go through and try to optimize.
I would try to think of the structure of your task, see what the
individual pieces
Thanks. I'm more familiar with Java than I am C++ so that might be worth
checking out. I also found this code with looks like it could do what I need:
http://finzi.psych.upenn.edu/R/Rhelp02/archive/15541.html
I've read some things about possible issues with it but it might work for my
use
I have data that looks like below and I would like to re-order the values
within the date column. I would like to have each year organized like so:
Sep-71
Oct-71
Nov-71
Dec-71
Jan-71
Feb-71
Mar-71
Apr-71
May-71
Jun-71
Jul-71
Aug-71
Sep-72
Oct-72
etc...
Is there any way I can order the column in
Dear help list - I have light data with 5-min time-stamps. I would like to
insert four 1-min time-stamps between each row and interpolate the light data
on each new row. To do this I have come up with the following code:
lightdata - read.table(Test_light_data.csv, header = TRUE, sep = ,) #
Hi,
I am trying to write a function which defines some arguments, then uses those
arguments as arguments of other function calls. It's a bit tricky to explain,
so a simple example will have to suffice. I imagine this has a simple solution,
but perusing through environments and other help lists
62 matches
Mail list logo