would be hugely welcome.
Many thanks,
Hugo
Hugo Pedder
Senior Research Associate
Population Health Sciences
University of Bristol
[[alternative HTML version deleted]]
___
R-packages mailing list
r-packa...@r-project.org
https://stat.ethz.ch/mai
elihood function is a handful but it does return a reasonable
looking vector of values. The "lnwagevec>-200" etc is an inelegant way of
preventing values of Inf and -Inf from entering the sum, the actual values
rarely go as high as 8 or low as -5.
Thank you in advance to anyon
I was looking function Distribution into stats package: it does not exist!
I am working with:
version
_
platform x86_64-pc-linux-gnu
arch x86_64
os linux-gnu
system x86_64, linux-gnu
status
major 3
minor 3.1
year 2016
How can a smacofIndDiff result be evaluated by a jackknife? If not now,
can it be expected?
The alternative seems to be a trip to Monte Carlo many gigaflops away,
hence a third question, is there experience and a model for parallel
running multiple smacofIndDiffs, using all available threads rath
Yes, it also occurs with WriteXLS version 3.2.1.
This test on several computers always leads to the same error.
Hugo Varet
2013/8/17 Rainer Hurling
> Am 13.08.2013 19:40, schrieb Hugo Varet:
> > Dear R users,
> >
> > I've just updated the WriteXLS package (on
y/3.0/WriteXLS/Perl/WriteXLS.pl"
--CSVPath "C:\Users\varet\AppData\Local\Temp\RtmpEzqFNz/WriteXLS" --verbose
FALSE --AdjWidth FALSE --AutoFilter FALSE --BoldHeaderRow FALSE --FreezeRow
0 --FreezeCol 0 --Encoding latin1 "C:\Users\varet\Desktop\irislatin1.xls"'
renvoie u
Dear Martin,
Thank you for your answer. I got the dist.binary() function in the ade4
package. It gives several distances for binary data. I agree with you, some
individuals probably share the same profile in my data.
Regards,
Hugo
[[alternative HTML version deleted
Dear Martin,
Thank you for your answer. Here is the exact call to agnes():
setwd("E:/Hugo")
library(cluster)
load("mydata.rda")
tableauTani<-dist.binary(mydata, method = 4, diag = FALSE, upper = FALSE)
resAgnes.Tani<-agnes(tableauTani, diss = inherits(tableauTa
74 0229
2 02350
3 120 0 15
So, I was wondering if it was normal that the C and FORTRAN codes give
different results?
Best regards,
Hugo Varet
[[alternative HTML version deleted
*(fit.ph$loglik[2]-fit.me$loglik[2]),3) # 0.9421373
1-pchisq(-2*(fit.ph$loglik[2]-fit.me$loglik[3]),3) # 0.5929387
If I'm not wrong, the likelihood ratio tests above indicate adding the
random component is not necessary, which fits well with the way I simulated
the data.
Thank you again,
H
d DevVariance
idp Vmat.1 0.1000 0.0100
Vmat.2 0.02236068 0.0005
Vmat.3 0.02449490 0.0006
The variances are round figures, so I have the feeling that the algorithm
didn't succeed in fitting the model.
Has anyone ever faced to this problem?
Tha
On Sun, 29 Apr 2012 11:27:06 -0700 (PDT)
arun wrote:
> I tried to get xyplot shrink fit for my mixed model._ But, the mixed model
> line is not seen in the graph.
> I would like to know if there is anything wrong with my code.
Typo? Try panel.abline(fixef(fm1),... or panel.abline(ff), ...
__
Dear list members,
what is the reason that one obviously can't do arithmetic operations on
zoo members with different index positions?
> require(zoo)
> z <- zoo(c(1,1,1),order.by=c(1,2,3))
> z
1 2 3
1 1 1
> z[1] + z[1]
1
2
> z[1:2] + z[1:2]
1 2
2 2
> z[1] +
alled. It installs and works seamlessly here on Linux, and also
works with Windows Vista, but I remember that I had to shuffle some dlls
from the postgresql directory. So if you provide some more information,
maybe I'll remember what I did to get it running.
Best regards
Hugo
Am 28.06.
the dataframe to lmer by using data=mydata; mainly
because passing a bunch of variables directly from the environment
hightens the risk of inadvertently using stale copies, but also because
str(mydata) would then produce a well arranged summary of all your
variables.
Kind regards
Hugo
On Sunday
algorithms.
Best regards
Hugo Mildenberger
On Saturday 25 June 2011 15:55:51 Duncan Murdoch wrote:
> On 11-06-25 11:33 AM, Mike Miller wrote:
> > I'm curious about what would cause this (see below), if it isn't a joke.
> > Is it possible that it didn't look ridiculo
Lusk,
if I were in your business, I would tackle such a problem with
help from a little formatting ...
e1=x1-mean(x1)
e2=x2-mean(x2)
n1=length(x1)
n2=length(x2)
N=(n1 + n2)
nu2=sum( c( ( x1 -mean(x1) )^2 , ( x2-mean(x2) )^2 ) )/N
ss=c(e1,e2)
b3=N*sum(ss^4)/(sum( ss^2)^2)
example still triggering the access violation, one could
try to run it under control of a memory guard like valgrind.
Best
Hugo
On Friday 24 June 2011 17:12:26 Saptarshi Guha wrote:
> Hello,
>
> I'm trying to use rjson to parse a JSON object. The object can be found here
>
>
mailing list for linear -
mixed -effect
models:
r-sig-mixed-mod...@r-project.org
But with numeric problems like this one you really should provide a self
contained
example (i.e. with data), at least the output of str(yourdata), if ever
possible.
Kind regards
Hugo
On Friday 24 June 2011
One possibility involves applying a regular expression via gsub:
as.numeric(gsub("^\\[1\\] ","",fvals))
On Friday 24 June 2011 14:39:31 Ravi Varadhan wrote:
> Hi,
>
> I would like to know how to capture the console output from running an
> algorithm for further analysis. I can capture thi
ations; solutions for larger
lambdas returned
So there seem to be many problems worth of a bug report to the author(s) of
glmnet.
Kind regards
Hugo
--- cut here ---
diff -ur a/glmnet/configure.in b/glmnet/configure.in
--- a/glmnet/configure.in 2011-06-15 20:05:30.0 +0200
+++
Please, can you help me with the following? I want to include a Greek letter
as part of a factor label that I need to use as factor on a xyplot. Probably
this was already solved so can I get the place where I can look for?
Thanks,
Hugo
[[alternative HTML version deleted
_0.5-9.tar.gz
patch -p1 < rmpi-check-dlopen-status.patch
tar -cvzf Rmpi_0.5-9a.tar.gz
R CMD INSTALL --preclean Rmpi_0.5-9a.tar.gz
Best
Hugo
cut
diff -ru a/Rmpi/src/Rmpi.c b/Rmpi/src/Rmpi.c
--- a/Rmpi/src/Rmpi.c 2010-11-30 20:38:01.0 +0100
+++ b/Rmpi/src/Rm
same CPU. So it seems that the
> processes are not distributed across the 8 CPUs available.
sfSetMaxCPUs(number=8) should fix at least the thread distribution. However, on
my
notebook all cores are kept busy with 100% load after sfInit(), with nothing
obvious to
compute. I don't kno
/mtest. If the link process failed, then the openmpi
installation was
inconsistent.
Kind regards
Hugo
On Wednesday 15 June 2011 17:45:41 Unger, Kristian, Dr. wrote:
> Dear Hugo
>
> I ran the command with the verbose switch and get the following output:
>
> > mpicc mtest
ried to run install.packages("Rmpi").
Kind regards
Hugo
On Wednesday 15 June 2011 16:22:07 Unger, Kristian, Dr. wrote:
> Thanks Hugo.
>
> I am pretty sure openmpi is installed:
>
> # zypper se openmpi
> Loading repository data...
> Reading installed packages..
n CRAN is Rmpi_0.5-9.tar.gz.
Best
Hugo
R CMD INSTALL --configure-args="--with-Rmpi-include=/usr/include
--with-Rmpi-libpath=/usr/lib64/openmpi --with-Rmpi-type=OPENMPI"
Rmpi_0.5-9.tar.gz
* installing to library ‘/home/hm/R/x86_64-pc-linux-gnu-library/2.13’
* installing *s
Hmm,
looks like there was a trailing blank after the backslash and before end of
line,
resulting in --with-Rmpi-libpath possibly not recognised:
> \--with-Rmpi-libpath=/usr/lib64/mpi/gcc/openmpi/lib64 \
I also doubt there is real need to escape newlines within a string. But another
possib
Andreas,
the jpg files you linked below do not exist, but if all you need
for the moment is a predictive model and graphical displays of
the fitted model and the calibrated sample data, then the R -
package "calib" will do it very well. Usage is very simple.
Best,
Hugo
On Tuesd
Ogbos,
I don't use Ubunto too, but X-11 headers should be installable running
apt-get xorg-dev
and/or (possibly)
apt-get libx11-dev xlibs-dev
Best
On Thursday 26 May 2011 16:13:36 Duncan Murdoch wrote:
> On 26/05/2011 10:05 AM, ogbos okike wrote:
> > Hi Dunc
With PostgreSQL at least, R can also be used as implementation
language for stored procedures. Hence data transfers between
processes can be avoided alltogether.
http://www.joeconway.com/plr/
Implemention of such a procedure in R appears to be straighforward:
CREATE OR REPLACE FUNCTION
Thomas,
your example works with R-2.13.0 and Apache running on
Gentoo. You may try to use "http://127.0.0.1:5984/path/fn.R";,
just to exclude that this problem has anything to do with
name resolution.
Best
Hugo
> Dear List,
>
> I have problems with the function source() u
William,
num_obs obviously isn't a vector, therefore length(num_obs) will
evaluate to one. Hence your for loop control part will expand to
for (i in 1:1)
while it should probably read:
for (i in 1:num_obs)
Best
Hugo
On Wednesday 18 May 2011 17:18:15 armstrwa wrote:
Hello Esther,
you left out the ellipsis argument (...) to the panel function.
That argument serves a placeholder for the rest of the parameters
you did not want to name and pass explicitely. If you don't pass the
ellipsis to panel.Dotplot, all unnamed parameters will get default
values. You can
This should be much faster:
idx <- c(1,2,3)
res <- dbGetQuery(con,paste("select * from db where id in
(",paste(idx,collapse=","),")",sep=""));
print(res);
On Wednesday 09 February 2011 14:13:17 Hugo Mildenberger wrote:
> Arne,
>
&g
operator:
"select * from db where id in (1,2,3);"
but I can't figure it out at the moment.
Best
Hugo
On Wednesday 09 February 2011 10:09:34 Arne H. Schulz wrote:
> Dear list,
> I'm working with a MySQL-Table and would like to extract selected ids that
> are s
d in
http://www.cs.uwaterloo.ca/~shai/TALKS/NIPS07_prob_wkshp.pdf
But this is an uneducated guess only.
Best
Hugo
On Saturday 05 February 2011 00:21:01 DB1984 wrote:
>
> Greg, Dennis - thanks for your input, I really appreciate the feedback, as it
> is not easy to source.
>
> I
Hello Chris,
You may also use the R-package "calib".
Hugo
On Tuesday 01 February 2011 17:08:13 Christopher Anderson wrote:
> Hello,
>
> I am trying to fit my Elisa results (absorbance readings) to a standard
> curve. To create the standard curve model, I performed a
On Sunday 23 January 2011 20:52:21 nandan amar wrote:
> I am trying to use function garsim from package gsarima.
> I can download gsarima package from
> http://cran.fyxm.net/web/packages/gsarima/index.html for linux.
> But how I can add it to R.
>
For example by entering
R CMD INSTALL /pat
Ajay,
you may try to run R as root while updating packages. Some
packages insist to install in directories writeable only by root.
On Saturday 22 January 2011 07:50:10 Ajay Ohri wrote:
> Dear List
>
> I use synaptic to download R on my Ubuntu 10.10. It seems latest version of
> R on Ubuntu i
x27;s bugs database for long.
Maybe getGraphicsEvent still runs on Windows.
Best
Hugo
On Friday 21 January 2011 13:26:26 christiaan pauw wrote:
> HI Everybody
>
> Does anyone know of documentation about different ways of obtaining user
> input in R. I have used readline() but I
ise" Linux after
all, so calling in Redhat support at that point and ask them to do
something for their money should be an option.
Best
Hugo
On Wednesday 19 January 2011 12:54:38 Mauricio Zambrano wrote:
> Dear Hugo,
>
> I tried your memory test program (without further mo
ut if you are still worried simply
change "void main" into "int main()" and add the statement " return 0; "
just before the last curly bracket (without quotation marks, of course):
int main() {
[...]
return 0;
}
Best
Hugo
ALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[37] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[49] FALSE
>length(z[time(z) >= 2 && time(z) <= 7:55])
>[1] 0
Your while loop was therefore equivalent to
>m=vector("numeric",0);
>while
Altay, simply run your tests under control of an exception handler:
help(try)
help(tryCatch)
On Monday 17 January 2011 22:05:07 Altay wrote:
> Hi, everybody.
>
> I am working processing EEG data from 1000 pacients. I have a specific
> syntax to perform the Spectral Analysis and a loop to
t;,size);
}else {
fprintf(stderr,"allocation of %lu bytes failed:%m\n",size);
}
}
put this into a file named, say, "tmalloc.c" and compile it using
gcc tmalloc.c -o tmalloc
Hugo
On Monday 17 January 2011 16:42:43 Mauricio Zambrano wrote:
&g
How about
ls(pattern="NY")
#rm(list=ls(pattern="NY"))
On Sunday 16 January 2011 22:08:16 Erin Hodgess wrote:
> Dear R People:
>
> I have the following:
>
> > ls(pattern="NY*")
> [1] "CRAN_df" "CRAN_df0" "CRAN_df1" "CRAN_mat" "CRAN_sp"
> [6] "CRAN_spdf1" "CRAN_spdf2
Unique certainly works well with Poisson distributed data, but is
essentially a no-op when confronted with continous floating-point
numbers.
Best
Hugo
> prepanel.rootogram
function (x, y = table(x), dfun = NULL, transformation = sqrt,
hang = TRUE, ...)
{
plot.line <- t
Using lattice and the rainfall$Time series as proposed below by Dennis
gives also a nice result:
rainfall$Time <- seq(from = as.Date('1993-01-01'),
to = as.Date('2007-12-01'), by = 'month')
xyplot(rainfall~Time,data=rainfall,type=c("g","p","l","smooth"))
On
Using R-2.12.1 and latticeExtra-0.6-14, I would like to understand
why a rootogram displaying samples from the Poisson distribution looks like I
expected it, whereas a rootogram using the normal distribution does not:
library(latticeExtra)
rootogram(~rpois(1000, lambda = 50), dfun = function(x)
interaction with the random factor, coefficients are 0
multinom(formula = mc$c ~(1 | mc$mom %in% mc$date) + (1 | mc$date) + mc$lma
+ mc$poid + (mc$pop) + mc$male %in% (mc$pop), family = "multinomial")
work well, but is it true ?
hugo mathé hubertEmail : hug...@hotmail.fr
app
5.1 322322 322322 322
5.2 558558 558558 558
5.4 217217 217217 217
5.5 217217 217217 217
Date: Fri, 5 Feb 2010 05:36:12 -0800
Subject: Re: [R] mlogit
From: djmu...@gmail.com
To: hug...@hotmail.fr
Hi:
On Fri, Feb 5, 2010 at 5:23 A
elMode, shape =
"long", alt.var = "mode", method = "nr", print.level = 0)
Coefficients:
alttrainaltbusaltcar time timeair
-3.6776 -4.2339 -4.5232 -0.5963 -2.7048
work prefectly
so, I don't understant any thing, ...
Is anybody understand w
53 matches
Mail list logo
