Hi Sacha,
I was referring to the BRugs package, where I usually proceed as follows:
result = BRugsFit (...)
param_samples = samplesSample(param)(1)
But (1) takes a long time for big chains. So my question was specific to
BRugs, if there an alternative way to get the samples.
openbugs()
Hi,
I am trying to use the system2 function to execute external
applications(actually it is psql) inside R 2.14.0 on X64 Windows, but
the psql command has escape character inside it's full command line,
can you help to figure out the correct parameters for the system2
function,
The working
Hi robin,
I'm not sure is what you need, but that's an esthetically nice
solution (one single line without any loop :) )
matrix(apply(log(cbind(as.numeric(a),as.numeric(b),as.numeric(c),as.numeric(d))),1,sd),3)
hope it could help,
PF
On Mon, Dec 12, 2011 at 5:15 PM, Robin Cura
Emma,
If, as you say, each unit is the same you can just repeat the units to
obtain the required number of units. For example,
unit_size - 10
n_units - 10
unit_id - rep(1:n_units, each=unit_size)
pid - rep(1:unit_size, n_units)
senior - ifelse(pid = 2, 1, 0)
pop -
arunkumar akpbond007 at gmail.com writes:
Please any one help in finding the predicted value for lmer function
model- lmer(formula =formula,data=data,verbose=TRUE,family = gaussian)
I need to get predicted value for this model.
Please see the code at http://glmm.wikidot.com/faq ,
Have a left-censored dataset, attempting to use a Tobit model and am
working
with the censReg package. I like how easy it is to move from glm
models to
predictions with 'predict' and wanted to ask if there was a way to do
so
with objects output from 'censReg'.
You can fit tobin
Dear all,
I am using the ISOdatetime function like this:
test.info$TradeTime = with(test.info, mapply(FUN = ISOdatetime, Year, Month,
Day, Hour, Minute, 60, EST))
Where Year, Month etc are all numeric.
I think ISOdatetime should return a POSIXct object. However, the result I
obtained from
Hi, I am looking for options to plot the following type of matrices:
A B C D
A A C C
A A A C
as a image like this: http://www.phaget4.org/R/image002.jpg
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Dear Jean,
That is a simple and effective solution. Thank you!
Best,
Shige
On Tue, Dec 13, 2011 at 4:16 PM, Jean V Adams jvad...@usgs.gov wrote:
Shige Song wrote on 12/12/2011 07:49:39 PM:
Dear All,
I am fitting a simple GAM model using the gam function in the mgcv
package. The only
I may be confused, but this seems ill-defined since, in the example,
the color corresponds to a numerical value: characters are only there
as labels. Do you intend to map the characters onto some numerical
spectrum? If so, how?
As a side note, one of the best places to get visualization ideas in
On 14/12/2011 16:24, Alex Zhang wrote:
Dear all,
I am using the ISOdatetime function like this:
test.info$TradeTime = with(test.info, mapply(FUN = ISOdatetime, Year, Month, Day, Hour,
Minute, 60, EST))
Where Year, Month etc are all numeric.
I think ISOdatetime should return a POSIXct
Thank you very much, Prof Ripley. The problem is solved and my understanding is
improved.
Happy holidays!
- Alex
From: Prof Brian Ripley rip...@stats.ox.ac.uk
To: Alex Zhang alex.zh...@ymail.com
Cc: r-help@R-project.org r-help@r-project.org
Sent:
Forgot to cc the list.
M
On Wed, Dec 14, 2011 at 12:19 PM, R. Michael Weylandt
michael.weyla...@gmail.com wrote:
It's a bit of a hack-job, but this might get you started:
d - structure(c(A, A, A, B, A, A, C, C, A, D,
C, C), .Dim = 3:4)
plot(0, type = n, xlim = c(0,ncol(d)), ylim = c(0,
Dear Jan,
Thanks for your reply.
The first solution works well for my needs for now, but I have a question about
the second. If I run your code and then call the function:
generate_unit(10)
I get an error that
Error in unit$size : $ operator is invalid for atomic vectors
Did you experience
Actually, scratch that, sorry!
I put the second part of your second solution code into a function and get the
right data frame in the end. So:
generate_unit- function(unit) {
pid- 1:unit$size
senior- rep(0, unit$size)
senior[sample(unit$size, 2)] - 1
return(data.frame(unit_id=unit$id, pid=pid,
Emma,
That is because generate_unit expects a data.frame with one row and
columns id and size:
generate_unit(data.frame(id=1, size=10))
Jan
Emma Thomas thomas...@yahoo.com schreef:
Dear Jan,
Thanks for your reply.
The first solution works well for my needs for now, but I have a
Dear List,
I am unsure if this is the correct list to post to, if it isn't I apologise.
I am using SSH to access a Linux version of R on a remote computer as it
offers more memory and processing power. The model will take 1-2 days to
run, I am accessing R through Putty and when I close the
Thank you,
I guess it didn't work for me, maybe is not possible?
I've tried:
con- odbcDriverConnect(Driver=SQL Server;
Server=...\\...;Database=...;Uid=...;Pwd=... ;)
v1=sqlQuery(con, select v1 from sqltable where v3 =cte and v2 in (select
v2 from R_dataframe) order by (select v2 from
Sorry for some typo last mail.
I corrected it,and resent.Sorry for it.
Sir:
I find out that for 2 level factor, if I set it to factor, then I'll get
error reply.
For the instance last mail, if I use:
result1-glm(y ~ factor(gender),family = binomial);#gender has 2 levels
Hi R users,
I'm using the 'bam' function in mgcv to examine trends in a remotely sensed
vegetation index. I have one random effect variable, 'cons' (with six
levels) which identifies different subjects within this analysis.
My model is specified as follows:
Hi
I am trying to download and unzip a zip file from the following location.
http://www.nseindia.com/content/historical/DERIVATIVES/2011/DEC/fo08DEC2011bhav.csv.zip
I initially there was issue with useragent so i ran the following code
options(HTTPUserAgent=Firefox/8.0.1)
u =
I have one equation, two unknowns, so I am trying to build the solution set
by running through possible values for one unknown, and then using uniroot
to solve for the accompanying second solution, then graphing the two
vectors.
p0 = .36
f = function(x) 0.29 * exp(5.66*(x - p0))
f.integral =
Dear all,
here is a simple problem that surely you already come across, but is giving
me a big headache...
I have a dataframe like this:
set.seed(3)
mydata - data.frame(var = rnorm(100,20,1),
temp = sin(sort(rep(c(1:10),10))),
subj =
Good Morning ,
I want to sample the following time series for every third hour. For example at
00:00,03:00,06:00,09:00 etc.
2011-01-01 00:00:00 0.00e+00
2011-01-01 01:00:00 1.471667e+01
2011-01-01 02:00:00 1.576667e+01
2011-01-01 03:00:00 0.00e+00
2011-01-01
All,
Given the following commands:
ag.m35-read.table(m35.txt,header=TRUE,sep=,)
ag.m35.lp-subset(ag.m35, race==lp)
aov.m35.lp=aov(year~time,data=ag.m35.lp)
anova.m35.lp=anova(aov.m35.lp)
anova.m35.lp
Analysis of Variance Table
Response: year
Df Sum Sq Mean Sq F value Pr(F)
time
Duncan,
If you are taking suggestions for expanding the tables package (looks great)
then I would suggest some way to get the tables into MS products. If I create
a full output/report myself then I am happy to work in LaTeX, but much of what
I do is to produce tables and graphs to clients
Colors probably are not the best for so many levels and combinations. Look at
the symbols function (or the my.symbols and subplot functions in the
TeachingDemos package) for ways to add symbols to a map showing multiple
variables.
-Original Message-
From: r-help-boun...@r-project.org
On Wed, 14 Dec 2011, Chris Mcowen wrote:
As a Linux newbie, I was wondering if anybody here knew how to keep R
running and interactive and return to it on a later date?
Chris,
Screen will retain the session and keep the process going. Read 'man
screen'.
Whether this works from a remote
What I did in the past (not with R scripts) is to start my jobs using
at (start the job at a specified time e.g. now) or batch (start the
job when the cpu drops below ?%)
at now R CMD BATCH yourscript.R
or
batch R CMD BATCH yourscript.R
something like that, you'll have to look at the man
On Wed, Dec 14, 2011 at 4:35 AM, Chris Mcowen chrismco...@gmail.com wrote:
Dear List,
I am unsure if this is the correct list to post to, if it isn't I apologise.
I am using SSH to access a Linux version of R on a remote computer as it
offers more memory and processing power. The model
kchkchkch wrote
I have one equation, two unknowns, so I am trying to build the solution
set by running through possible values for one unknown, and then using
uniroot to solve for the accompanying second solution, then graphing the
two vectors.
p0 = .36
f = function(x) 0.29 *
1) Use dput() to submit data.
2) Would this work? (It requires your data are evenly spaced, but I
think that's it) d[seq(1, nrow(d), by = 3), ]
Michael
On Wed, Dec 14, 2011 at 7:17 AM, abcdef ghijk lineh...@yahoo.com wrote:
Good Morning ,
I want to sample the following time series for every
Greg,
Please look at the SWord package. This package integrates MS Word with R
in a manner similar to the SWeave integration of LaTeX with R.
Download SWord from rcom.univie.ac.at
If you have a recent download of RExcel from the RAndFriends installer, then
you will already have SWord on your
On Dec 14, 2011, at 11:43 AM, matteo wrote:
set.seed(3)
mydata - data.frame(var = rnorm(100,20,1),
temp = sin(sort(rep(c(1:10),10))),
subj = as.factor(rep(c(1:10),5)))
and I need to make a scatter plot for each subj, not a problem, but...
what i want is
On Dec 14, 2011, at 10:01 AM, Tony Stocker wrote:
All,
Given the following commands:
ag.m35-read.table(m35.txt,header=TRUE,sep=,)
ag.m35.lp-subset(ag.m35, race==lp)
aov.m35.lp=aov(year~time,data=ag.m35.lp)
anova.m35.lp=anova(aov.m35.lp)
anova.m35.lp
Analysis of Variance Table
Response:
Am trying to install package JGR and, by necessity, JavaGD on a few
Linux platforms. R rev. 2.14.0. Had some success on one platform, but
am running into a problem on a new platform. I installed rJava, JavaGD,
iplots, and JGR. Things went - I thought - without a hitch. Now, when
I try to use
On Wed, 14-Dec-2011 at 08:43AM -0800, matteo wrote:
| Dear all,
|
| here is a simple problem that surely you already come across, but is giving
| me a big headache...
|
| I have a dataframe like this:
|
| set.seed(3)
| mydata - data.frame(var = rnorm(100,20,1),
| temp =
On Wed, 14-Dec-2011 at 09:55AM -0800, Peter Langfelder wrote:
| On Wed, Dec 14, 2011 at 4:35 AM, Chris Mcowen chrismco...@gmail.com wrote:
| Dear List,
|
|
|
| I am unsure if this is the correct list to post to, if it isn't I
apologise.
|
|
|
| I am using SSH to access a Linux version
I have a big SQL Server application that needs some help from R, and this R
needs to load a large workspace in order to service the database request. I'd
like to keep an R process running continually because loading the workspace
takes about 1 minute. The data in the workspace doesn't change
Heh. Yes,
Berend Hasselman wrote
Shouldn't the line with uniroot be (why try as function?)
sol = uniroot(fcn, lower = p1, upper = 1)
Before the line with for(i in 1:n) insert the following
fcn = function(p2) p1*f(p1) + (.20/5.66)*(exp(5.66*(p2 - p0)) -
exp(5.66*(p1 - p0))) + (1 -
On Wed, Dec 14, 2011 at 7:04 AM, agent dunham crossp...@hotmail.com wrote:
Thank you,
I guess it didn't work for me, maybe is not possible?
I've tried:
con- odbcDriverConnect(Driver=SQL Server;
Server=...\\...;Database=...;Uid=...;Pwd=... ;)
v1=sqlQuery(con, select v1 from sqltable where
Do unz / unzip get the job done?
Michael
On Wed, Dec 14, 2011 at 8:19 AM, himanshuy himanshu@gmail.com wrote:
Hi
I am trying to download and unzip a zip file from the following location.
http://www.nseindia.com/content/historical/DERIVATIVES/2011/DEC/fo08DEC2011bhav.csv.zip
I
On 14.12.2011 19:43, Cable, Sam B Civ USAF AFMC AFRL/RVBXI wrote:
Am trying to install package JGR and, by necessity, JavaGD on a few
Linux platforms. R rev. 2.14.0. Had some success on one platform, but
am running into a problem on a new platform. I installed rJava, JavaGD,
iplots, and
Richard,
I have looked at SWord before, but to my knowledge it does not deal directly
with the tabular objects created by the tables package (please correct me if I
am wrong). These objects do have a matrix as the main data, but the attributes
are different from the usual dimnames. There
There is also the problem that SWord's license does not allow for commercial
use. R2wd, write.table, and odfWeave don't have this restriction.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
From: Richard M. Heiberger
On Wed, Dec 14, 2011 at 4:43 PM, Greg Snow greg.s...@imail.org wrote:
There is also the problem that SWords license does not allow for
commercial use. R2wd, write.table, and odfWeave dont have this
restriction.
The free license for SWord does have the commercial use restriction.
There is
Hi all,
I have solved my question
Output- matrix(0, length(x), 1)
for (i in 1:length(x)){
files - paste('KT', i, '.csv', sep = '')
da1 - read.csv(files)
Output[i,1]-da1[1,2]
}
2011/12/14 Trying To learn again tryingtolearnag...@gmail.com
Hi all,
I have 100 csv files always with
First PRoblem solved
read.csv(KT80.csv)
X x
1 1 0.01331361
I ommitt in calling the file name...
2011/12/14 Trying To learn again tryingtolearnag...@gmail.com
Hi all,
I have 100 csv files always with this information (I Attach two example
excels)
KT80.csv contains:
Hello,
I am writing a test function, and would like to have it return an
appropriate 'htest'
object. This is very easy, but the test can be run for several orders, a
frequent situation
in time series tests.
It would be nice to return a data.frame with order, params, test statistic,
p.value.
Thanks David and Patrick,
really appreciate tour help...
m.
Il giorno 14 Dec 2011, alle ore 18:33, David Winsemius ha scritto:
On Dec 14, 2011, at 11:43 AM, matteo wrote:
set.seed(3)
mydata - data.frame(var = rnorm(100,20,1),
temp = sin(sort(rep(c(1:10),10))),
Hi everyone,
I have been following these forums closely for the last few months but this
is my first time posting. Basically I am trying to get an Oracle SQL
Developer Binary Large Object(BLOB/CLOB) file in as an R object in R. For
those not familiar with a BLOB/CLOB file, it is basically lik a
Dear R members,
I would like to store 18 graphics at the end of a loop (as we do for data
frames and arrays). For each iteration, I use the function x11() and I want
to keep my graphs in a single window to re-use them later. Do you know a
function for storing graphics ?
Many thanks
Julie
--
Hi all
a-c(2,Inf)
max(a)
[1] Inf
How can I say...omit Inf?
I cannot find anything...
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
range(a, finite=TRUE)[2] # [1] for the min
[1] 2
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Trying To learn
again
Sent: Wednesday, December 14, 2011 3:55 PM
To:
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of William Dunlap
Sent: Wednesday, December 14, 2011 4:19 PM
To: Trying To learn again; r-help@r-project.org
Subject: Re: [R] Omit Inf
range(a, finite=TRUE)[2] # [1] for the
Please follow the posting guide and provide some example code.
Bottom line, most functions that draw a graph return a graphic
object that, like any other R object, can be given a name in the
global environment, put into a list, etc.
Alternatively, depending on what you mean -- it was not
There is also a windrose (though not with that name) function in the
openair package.
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 12/12/11 1:20 PM, Adrienne Wootten amwoo...@ncsu.edu wrote:
Greetings!
I'm having an
library(openair)
?windRose
--
Clint BowmanINTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET: cl...@math.utah.edu
Department of Ecology VOICE: (360) 407-6815
PO Box 47600FAX:(360) 407-7534
That only puts different colors on different axes. I was wondering
about the (silly) possibility of having a range of colors for the ticks
on a given axis.
quote
From: Hans W Borchers hwborchers_at_googlemail.com
Date: Tue, 13 Dec 2011 14:21:31 +
Carl Witthoft carl at witthoft.com
Bert,
The OP may be looking for recordPlot/replayPlot.
We were not told the OS (x11 exists even on Windows). The standard
screen devices on Windows (windows() ) and Mac OS X (quartz()) have plot
histories that can be used to flick between plots. And 'flick' is the
operative word here: it
On 2011-12-14 18:32, Carl Witthoft wrote:
That only puts different colors on different axes. I was wondering
about the (silly) possibility of having a range of colors for the ticks
on a given axis.
I think you would have to either hack the code or issue
several axis() calls, e.g.
I tried using unzip but it gave the following error Invalid zip name
argument
options(HTTPUserAgent=Firefox/8.0.1)
u =
http://www.nseindia.com/content/historical/DERIVATIVES/2011/DEC/fo08DEC2011bhav.csv.zip;
o = getURLContent(u, verbose = TRUE, useragent = getOption(HTTPUserAgent))
after this
Hello,
I have a data frame with hourly or sub-hourly weather records that span
several years, and from that data frame I'm trying to select only the
records taken closest to noon for each day. Here's what I've done so far:
#Add a column to the data frame showing the difference between noon and
Dear R-users,
I am a relatively new R user and I have a simple question. Is there a
command attached to a R function, which performs matrix pooling?
Thanks in advance,
Jana
--
--
The world is full of mysteries. Life is one. The curious limitations of
finite minds are another.(J.B.S.
I am still interested in this. Is there no way to pay directly online? via
paypal or other?
Thanks,
Xu
--
View this message in context:
http://r.789695.n4.nabble.com/Re-Monetary-support-to-the-R-project-Was-Re-Executable-for-Production-Use-tp1585186p4198369.html
Sent from the R help mailing
Dear R helpers,
I have a question about drawing random numbers from many categorical
distributions.
Consider n individuals, each follows a categorical distribution defined
over k categories.
Consider a simple case in which n=4, k=3 as below
catDisMat -
Hi,
I use the function fpot of packages evd. If I call the fit that I obtain
fit, I want to modify the name of the axis and the main title that is
produced by plot(fit1). To do this, I want to create a new function where I
would have the names modified. The problem is that I can't find the
kchkchkch wrote
Heh. Yes,
Berend Hasselman wrote
Shouldn't the line with uniroot be (why try as function?)
sol = uniroot(fcn, lower = p1, upper = 1)
Before the line with for(i in 1:n) insert the following
fcn = function(p2) p1*f(p1) + (.20/5.66)*(exp(5.66*(p2 - p0)) -
68 matches
Mail list logo
