Hello,
I was under the impression that in a call to gls, if i specify
gls(.., cor=AR1(form= ~ 1 | subject))
and suppose there are N subjects, then i would have obtained estimates
for N Phi (s) (the AR1 coefficient) - but I only get 1 estimate.
I then checked "Linear Mixed-Effects Models Using
return to the prompt and the
print in telemStats should not appear
On Mon, Mar 16, 2015 at 4:02 PM, David Winsemius dwinsem...@comcast.net wrote:
On Mar 16, 2015, at 3:08 PM, Saptarshi Guha wrote:
Hello,
I would like a function X to return to the place that called the
function XParent
Hello,
I would like a function X to return to the place that called the
function XParent that called this function X.
Y calls XParent
Y = function(){
XParent()
print(hello)
}
XParent calls X
XParent = function(){
X()
print(H)
}
X returns to the point just after the call to XParent.
Hello,
Couldn't think of a better subject line. Rather than a matrix like
x,y
..,..
.,..
I have a matrix like
x,y,n,
..,..,..,
..,..,..
and so on. Also, sum(n) is roughly few hundred million. The number of rows
is 1MM
Are they routines to fit a 2d kde estimate to data provided in this form?
I
hello,
This
rbindlist(list(list(a=NA,b=NA),list(a=20,b=FALSE)))
returns
a b
1: NANA
2: TRUE FALSE
as per the documentation ?rbindlist
is there a way to specify the column class of 'a' to be numeric?
In actual usage, i wont be able to re-order the 2nd list entry to be the
Hello,
I find Map to be nice interface to mapply. However
Map calls mapply which in turn calls mapply via .Internal.
Is there a parallel version of mapply (like mcapply) or do I need to write
this myself?
Regards
Saptarshi
[[alternative HTML version deleted]]
two statements (a sequential
and a parallel statement) -- if a parallel backend is running, it will use
that, otherwise it'll execute in sequential mode.
--j
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] on
behalf of Saptarshi Guha
Hello,
Consider the following transcript
x=NULL
x$date=10
x
$date
[1] 10
or
x=NULL
x[10]=10
x
[1] NA NA NA NA NA NA NA NA NA 10
Wouldn't one expect that x remain NULL despite the further additions (i.e
the x$date, x[10]) etc? Is coercing appropriate here?
Cheers
Saptarshi
Hello,
I was looking at the R (installed on RHEL6) shell script and saw
R_HOME_DIR=/usr/lib64/R. Nowhere (and I could have got it wrong) does
it read in the environment value R_HOME_DIR. I have the need to rsync
the entire folder below /usr/lib64/R to another computer into another
directory
Hello,
If i were to embed R (in C), i could over ride
ptr_R_ShowMessage = Re_ShowMessage;
ptr_R_WriteConsoleEx =Re_WriteConsoleEx;
where Re_ShowMessage,Re_WriteConsoleEx are my own functions.
How would i do this in a package?
I would like the flexibility of having my function be called
Hello,
Is there an R function that tells me the packages that come with R e.g.
c(base,boot,methods,mgcv,...)
Thank you
Saptarshi
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
utils base
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 Saptarshi Guha
Sent: Tuesday, March 19, 2013 11:37 AM
To: R-help@r-project.org
Subject: [R] List
thanks much.
Saptarshi
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Hello,
How would i use bquote to create names for a list.
e.g. to create list(a=10)
bquote(list(.(X)=10), list(X=a))
does not work.
The best i could come up with is
bquote({ a=list(10);names(a)=.(X); a}, list(X=a))
which is quite ugly.
Is there an elegant way to solve this?
Regards
Hello,
[ Sorry, I sent the last email as HTML, this time it's in plain text ]
Suppose my variable,S, (time for something to start) is a skew-normal
distribution [1]. Can glm and mgcv handle this type of distribution for the
dependent variable?
Regards
Saptarshi
[1]
I'm trying to install a package on OS X (Snow Leopard) using the following
command, but instead of respecting the arch directive it install i386. How
can this be resolved?
R --arch=x86_64 CMD INSTALL
* installing to library
Hello,
I have a function 'ewrap' (see below for definition).
It takes 3 expressions and returns another expression e.g.
map - ewrap({
len - length(r$addon)
rhcollect(len,1)
})
becomes:
expression({
NULL
result - mapply(function(.index, k, r) {
{
len -
Hello,
(Sorry for the repost, i am resending in plain text)
I have a function 'ewrap' (see below for definition).
It takes 3 expressions and returns another expression e.g.
map - ewrap({
len - length(r$addon)
rhcollect(len,1)
})
becomes:
expression({
NULL
result -
Hello,
I am lattice scatterplot that has 2 panels (could be a few more). Both
panels have a y-axis label that is different on the
left and right sides. However the right hand side axis labelling depends on
which panel i'm in .
(I am plotting two curves in one panel, the left y-axis has the scale
Hello,
Is there a way to add ylab on the right hand side also (in lattice)?
Different from the left hand side?
Cheers
Saptarshi
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Hello,
I have the following code
asTheEconomist(bwplot(ver~quans,data=z,cex=0.5,subset=abs(quans)1.1,scales=list(x=list(tick.num=20,cex=0.55,alternating=3))
,xscale.components = xscale.components.CF
,main=list('FOOO',cex=0.8)))
The xscale.components displays the exponent of the
Hello,
I would like to write a function where substitute operates on the
parameter, but ...
Expression = function(o,l) substitute(o, l)
Expression({x=.(FOO)}, list(FOO=2))
o
How do i get substitute to work on the contents of o.
Regards
Saptarshi
Hello,
I'm trying to use rjson to parse a JSON object. The object can be found here
http://pastebin.com/np0s5hgM
(you'll probably need to add quotes around the content)
fromJSON returns the error
glibc detected home/sguha/lib64/R/bin/exec/R: realloc(): invalid next
size: 0x0209e210
Hello,
I notice that e.g /home/sguha/lib64 is hard coded into the /bin/R file .
I nstalled R as ./configure --prefix=$HOME ...
What i need to do is ship the entire R distribution to remote nodes,
and run R. These are shipped to ephemeral directories
so I dont know the path ahead of time.
R_HOME
Hello,
Given the following display
library(maps)
library(mapproj)
m - map('world',plot=FALSE)
map('world',proj='mollweide',bg=bgcolor[2],col='white')
map.grid(col=2,lim=c(-175,-175,-180,180),label=FALSE,lty=2)
Is there a way to color the region outside the dotted red border?
Thank you
Hello,
I should see how to mirror a google group (if indeed possible).
There are few things stopping from having RHIPE work on the remote
while the client sends queries from a Windows client.
But as of now, it does not work.
Regards
Saptarshi
On Fri, Mar 18, 2011 at 7:50 AM, XiaoboGu
. scan(textConnection(test),
what=character(0))
Hope this helps
Allan
On 16/03/11 22:18, Saptarshi Guha wrote:
Hello R users,
I have this regex see [1] for apache log lines. I tried using R to parse
some data (only because I wanted to stay in R).
A sample line is [2]
(a) I saved the line
Hello R users,
I have this regex see [1] for apache log lines. I tried using R to parse
some data (only because I wanted to stay in R).
A sample line is [2]
(a) I saved the line in [1] into ~/tmp/a.txt and [2] into /tmp/a.txt
pat - readLines(~/tmp/a.txt)
test - readLines(/tmp/a.txt)
test
Hello,
Consider this little function (not concerned about efficiency here)
m=function(r) { a-serialize(r,NULL); zz-rawConnection(raw(0),r+);
writeBin(length(a),zz); c(rawConnectionValue(zz),a)}
a=m(10)
p=rawConnection(a,r+); u=readBin(p,int)
unserialize(p)
Error in unserialize(p) : unknown
Hello,
I have this script which will be invoked as
Rscript a.r a.r
a.r follows
#!/usr/bin/Rscript --vanilla
f=file(stdin)
while(TRUE){
y=readLines(f,n=1,warn=TRUE)
if(length(y)==0) break else print(y)
}
But it only reads one line from a.r
How can I read line by line from standard
Hello,
Would one use gsub to perform pythons translate?
in python
table = string.maketrans(abc,xyz)
abc.translate(table)
xyz
how does one do this in R?
Thank you
Saptarshi
__
R-help@r-project.org mailing list
Yes, thanks much.
On Thu, Oct 14, 2010 at 9:36 AM, jim holtman jholt...@gmail.com wrote:
?chartr
Is this what you are looking for?
On Thu, Oct 14, 2010 at 12:10 PM, Saptarshi Guha
saptarshi.g...@gmail.com wrote:
Hello,
Would one use gsub to perform pythons translate?
in python
table
Dirk, thanks.
On Thu, Oct 14, 2010 at 9:42 AM, Dirk Eddelbuettel e...@debian.org wrote:
On 14 October 2010 at 11:21, Dirk Eddelbuettel wrote:
| On 14 October 2010 at 08:49, Saptarshi Guha wrote:
| | Hello,
| | I have this script which will be invoked as
| |
| | Rscript a.r a.r
Hello,
I have question about the overhead in lapply.
x is a list of 3000 lists. Each of the i (1=i=3000) list elements is
pair of two elements: a string vector and a data frame
x is roughly 235MB.
gc()
##
z - system.time(y - lapply(x,function(r){
system.time(serialize(r,NULL))['elapsed']
Hello,
I take this back, wrong implementation of the experiment.
Sorry for the noise.
Regards
Saptarshi
On Sun, Aug 15, 2010 at 12:22 PM, Saptarshi Guha
saptarshi.g...@gmail.com wrote:
Hello,
I have question about the overhead in lapply.
x is a list of 3000 lists. Each of the i (1=i=3000
Hello,
I have the following data frame, using lattice and qqmath, how do I display the
quantile plot (e.g distribution = qunif), but highlight the f==TRUE and f==FALSE
differently?
library(lattice)
pp - function(n) ((1:n)-0.5)/n
x=data.frame(x=runif(20),f=sapply(runif(20),function(r) r0.5))
I
Hello,
Do you need /all/ the data in memory at one time? Is your goal to
divide the data (e.g according to some factor /or/ some function of
the columns of data set ) and then analyze the divisions? And then,
possibly, combine the results ?
If so, you might consider using Rhipe. We have analyzed
Dear R Users,
Using codetools I obtained the text representation of the parse tree for this
snippet
z=quote({x[1]-2})
showTree(z)
({ (- ([ x 1) 2)) (A)
If I understand correctly, x[1]-2 ought to be [-(x,1,2), so shouldn't i see
({ ( [- x 1 2 ) )
If indeed the parse tree in (A) is
Hello,
I would like to compose a series of functions each taking one
argument. Let us assume the values returned by the functions are
contained within the domains.
So
b-function(r) r[-1]
compose(b,b)
b(b(x))
This is my solution, is there something succinct ?
Regards
Saptarshi
compose -
Hello,
I wrote a package, which in the NAMESPACE file exports functions like this:
exportPattern(^\\rh)
On R-2.8 (Linux, 64), upon loading the package I have the rh functions present.
On R-2.10, Mac OS X, (32 bit), it builds, loads, but the functions are not
loaded, i.e the only function is
As a followup, i can do things Package:::rh*, there just no exported.
On Tue, Dec 1, 2009 at 11:18 AM, Saptarshi Guha
saptarshi.g...@gmail.com wrote:
Hello,
I wrote a package, which in the NAMESPACE file exports functions like this:
exportPattern(^\\rh)
On R-2.8 (Linux, 64), upon loading
wrote:
You probably just want : exportPattern(^rh)
What did you expect the backslash to do ?
On 12/01/2009 05:18 PM, Saptarshi Guha wrote:
Hello,
I wrote a package, which in the NAMESPACE file exports functions like
this:
exportPattern(^\\rh)
On R-2.8 (Linux, 64), upon loading the package
x - foo(10)
eval(x)
Error in eval(expr, envir, enclos) : object 's' not found
Q: How do I force a substitution so that the returned expression has the value
of 's' 'hardcoded' in it e.g
foo(10)
returns
expression({
v - 10
print(v)
})
Regards
Saptarshi
--
saptarshi guha | http
The substitute() function can do it, but the bquote() function is the most
convenient:
foo - function(s) {
+ bquote(expression({
+ v - .(s)
+ print(v)
+ }))
+ }
foo(10)
expression({
v - 10
print(v)
})
Duncan Murdoch
Thanks, after sending this email, i used
I have some beginner's questions regarding local, in the docs, it says that
local evaluates an expression in a local environment.
Q1: why is B different from A? In B, is a-a+1 getting evaluated
before eval proceeds?
#A
a=0
eval(quote(a-a+1),new.env())
a # 0
#B
a=0
eval(a-a+1,new.env())
a # 1
: it is possible to
extend this to serialize functions, expressions, environments and
several other objects.
However that is some time in the future.
Regards
Saptarshi
[1] http://code.google.com/apis/protocolbuffers/docs/overview.html
Download:
http://ml.stat.purdue.edu/rproto_1.0.tar.gz
Saptarshi Guha
Hello,
(Previous message was scrubbed)
I apologize if this is not the list for package announcements. Having
said that I've released a package called rprotos which implements a
simple
serialization using Googles protocol buffers[1].
The package also includes some miscellaneous functions for
Hello,
I have used socketConnection to connect to a TCP server. I havent
figured out a way to do the same with a UDP server.
i.e I have a server listening on 9000, communicating via UDP. I would
like to , from R, send packets to this server,
This does not work
u -
Hello,
I tried this pseudo-generator style experiment (see below). The -
operator assigns to in the calling environment which would be the
environment of getN.
Yet when the function incr is returned, isn't this environment lost?
Also the print displays GlobalEnv, yet the globalenv does not have
-project.org] On
Behalf Of Saptarshi Guha
Sent: Friday, 17 July 2009 5:16 PM
To: R-help@r-project.org
Subject: [R] Remembering a value in multiple calls of a function
Hello,
I tried this pseudo-generator style experiment (see below). The -
operator assigns to in the calling environment
Hello,
On Fedora 8, yum install R , installs R-2.6.
I dont have much experience with yum, how do change the repos(?) to
install R-2.9?
Thank you
Saptarshi
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Hello,
I am using EC2 to launch several instances. On each instance, I
perform the following commands
yum -y install R (installs R-2.6)
R CMD javareconf
sleep 10
wget rpackage.tgz
R CMD INSTALL rpackage.tgz
Now, rpackage.tgz needs to be built with jni libraries. When the
instance is fully up, i
Hello,
I'm using R-2.9 and have a question about integer ranges. On C, (OS X
platform),
LONG_MIN and LONG_MAX(from limits.h) correspond to -2147483648(-2^31)
and 2147483647 (2^31-1)
In Java(1.4), Integer.MIN_VALUE and Integer.MAX_VALUE have the same values.
On R, as.integer(2^31-1) is the
You've hit on the reason -- representation of NA
src/main/arithmetic.c:157: R_NaInt = INT_MIN;
Aha! Yes, I missed that.
(R-devel seems more appropriate for this question).
Yes, i was not sure.
Thank you
Saptarshi
__
R-help@r-project.org
://ml.stat.purdue.edu/rhipe
There are some downsides to RHIPE which are described at
http://ml.stat.purdue.edu/rhipe/install.html#sec-5
Regards
Saptarshi Guha
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
, I
haven't
had the chance to compare with rmpi and snow.
TODO(will appear soon ): another function: rhmr - mapreduce using R.
Regards
Saptarshi Guha
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch
, I
don't get any such error messages.
Is there a possible reason why?
Regards
Saptarshi
Saptarshi Guha
__
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
the baseline and same height, if necessary the paper width and height
can be adjusted
( i tried setting the paper width and height to no avail).
Is the way to base align the two figures? Do I have to get down the grid level?
Regards
Saptarshi Guha
__
R-help@r
understand it. What does the
option do and why did I initially run into an error ?
From man R
--max-nsize
Set max number of cons cells to N
Why should there be a maximum?
Thanks for your time
Saptarshi Guha
__
R-help@r-project.org mailing list
https
information and examples please visit this page: http://www.stat.purdue.edu/~sguha/rhipe
.
If anyone would like to contribute to this project, please email me
directly - any help is welcome.
Regards
Saptarshi Guha
__
R-help@r-project.org mailing
to this instead of iterating over all objects E
and assigning into .GlobalEnv?
Thank you
Saptarshi
--
Saptarshi Guha - saptarshi.g...@gmail.com
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
),con=conn,endian=big) #conn opened in binary
But the response on the server side is wrong, eg. 1 is received as
4294967297
How do I then write long across the network from R?
Thank you
Saptarshi
Saptarshi Guha | saptarshi.g...@gmail.com | http
]
*** REXPFactory unable to interpret
[EMAIL PROTECTED] *
** REXPFactory unable to interpret
[EMAIL PROTECTED]
How can I in general,force the unserialize to return a list?
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
Marijuana is nature's way of saying, Hi
Saptarshi
R.version = 2.8
On Nov 30, 2008, at 5:53 PM, Saptarshi Guha wrote:
Hello,
I have some code which generates lattice objects. The function
recieves serialized forms of the lattice objects which it then
unserializes and then adds to an ArrayListREXP.
REXPRaw
dataset, the y-axis
labels overlap so I would like to split it across several panels(1
panel per page) with about 25-50 per panel.
Is there a way to do this?
Regards
Saptarshi
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML
for each panel based on the pulled in data and
plotting that.
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo
.)
Is there a way to solve this? Also, each level of C appears in exactly
one
combination of A and B (the combinations that appear in the data set,
not the ones produced by lattice)
Kind regards Saptarshi
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative
))
However, I can see the line in the bottom panel, but nothing in the
top. Is there some
clipping going on?
Thanks
Saptarshi
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML version deleted]]
__
R-help@r
Got it. It should be plot1.panel.1.2.off.vp and likewise.
Thanks
Saptarshi
On Sep 19, 2008, at 12:08 PM, Saptarshi Guha wrote:
Hello,
I have a multi-page display each consisting of two-panels above each
other.
I need to draw a line from the top panel to bottom panel. Using
Forgot to send to the list.
Begin forwarded message:
From: Saptarshi Guha [EMAIL PROTECTED]
Date: September 19, 2008 1:43:50 PM GMT-04:00
To: Gabor Grothendieck [EMAIL PROTECTED]
Subject: Re: [R] Lines between panels in lattice
Certainly.
x=cbind(runif(10),rnorm(10),c(rep(1,5),rep(2,5
)
but how can I get around it? Should I specify a width to the device?
Thank you
Saptarshi
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https
the bottom of the curve and the lower horizontal axis and
between the maximum values of the curve and the upper horizontal axis.
How can i increase the padding - top,bottom,left and right?
Thanks in advance for your time.
Regards
Saptarshi
Saptarshi Guha | [EMAIL
Hello,
Is there a way to add minor tick marks to the Y-axis of a lattice
plots?
Thank you
Saptarshi
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML version deleted]]
__
R
Hello,
I need to plot a line graph using lattice so that the endpoints are
flush with the left and right vertical axes.
So the scale end points and horizontal axes length should coincide.
This works
u=newdata[1:100,] #u[,1] is from 1700-1795
and
the
vertical axes.
How can I remove this padding?
A reproducible example would have helped, but perhaps you want
scales = list(x = list(relation = free, axs = i))
-Deepayan
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
in prepanel but what?
Is there any function i can call to set the actual length of the axes?
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r
scales.
I suppose i need to do something in prepanel but what?
Is there any function i can call to set the actual length of the axes?
Thank you
Saptarsho
Saptarshi Guha | [EMAIL PROTECTED] | http://www.stat.purdue.edu/~sguha
[[alternative HTML version
Thank you.
On Feb 17, 2008, at 1:18 PM, Deepayan Sarkar wrote:
See ?print.trellis.
p - xyplot(1 ~ 1, aspect = 0.5)
p
plot(p, panel.width = list(2, inches), panel.height = list(1,
inches))
This overrides 'aspect=' though, that is, you will need to specify
both width and height if you
Hello,
I have an optimization problem which I would like to solve in R but
am not sure how to approach this.
I have a cost function C which is quadratic in rx and ry (if it
helps, the cost function is separable in rx and ry) and is always
positive for all rx,ry.
I
79 matches
Mail list logo