Dear R forum,
I am trying to execute following code (Page no 259 - VGAM.pdf)
#
.
library(VGAM)
set.seed(123)
fdata - data.frame(y1 = rfrechet(nn - 1000, shape = 2 + exp(1)))
Hi
We have a few questions regarding the use of the isoMDS function.
When we run isoMDS function using 60,000 x 60,000 data matrix,
we get the following error message:
cmdscale(d, k) : invalid value of 'n'
Calls: isoMDS - cmdscale
Hello,
Can I get help? I have the attached kind of material and I should create a
time-dependant covariate from the variate Optime? I have tried using
survsplit function but I don't know exactly how to do it. Or I can't get it
to work. How can I split every id's follow-up time into two rows
On 06/11/2014 06:04, Katherine Gobin wrote:
Dear R forum,
I am trying to execute following code (Page no 259 - VGAM.pdf)
#
.
library(VGAM)
set.seed(123)
fdata -
Dear Mr Michael,
Thanks a lot for your guidance. The pdf file describing VGAM package has
mentioned 'frechet' in the example, so I got the error.
Regards
Katherine
On Thursday, 6 November 2014 2:54 PM, Michael Dewey i...@aghmed.fsnet.co.uk
wrote:
On 06/11/2014 06:04, Katherine Gobin
Dear Ms. Gobin,
There would appear to be a typo in the package manual that appears on
CRAN. Doing ?Frechet (it would be nice to have this aliased also to
frechet) points you immediately to frechet. With frechet2
substituted for frechet your code, everything works.
cheers,
Rolf Turner
You avoid the call to cmdscale() by supplying your own starting configuration
(see the manual page for the y= argument). You could still hit other barriers
within isoMDS() or insufficient memory on your computer.
-
David L Carlson
Department of Anthropology
I wish to simulate the following stochastic process, for i = 1...N
individuals and t=1...T periods:
y_{i,t} = y_0 + lambda Ey_{t-1} + epsilon_{i,t}
where Ey_{t-1} is the average of y over the N individuals computed at time
t-1.
My solution (below) works but is incredibly slow. Is there a faster
Matteo,
I tried your example code using R 3.1.1 on an iMac (24-inch, Early 2009), 3.06
GHz Intel Core 2 Duo, 8 GB 1333 MHz DDR3, NVIDIA GeForce GT 130 512 MB
running Mac OS X 10.10 (Yosemite).
After entering your code, the elapsed time from the time I hit return to
when the graphics appeared was
On 2014-11-05 14:50, Therneau, Terry M., Ph.D. wrote:
This is fixed in version 2.37-8 of the survival package, which has been
in my send to CRAN real-soon-now queue for 6 months. Your note is a
prod to get it done. I've been updating and adding vignettes.
Is your fixed code publicly
Matteo,
Ah — OK, N=20, I did not catch that. You have nested for loops, which R is
known to be exceedingly slow at handling — if you can reorganize the code
to eliminate the loops, your performance will increase significantly.
Tom
On Thu, Nov 6, 2014 at 7:47 AM, Matteo Richiardi
I find that representing the simulated data as a T row by N column matrix
allows for a clearer and faster simulation function. E.g., compare the
output of the following two functions, the first of which uses your code
and the second a matrix representation (which I convert to a data.frame at
the
Hi,
I am working R on windows 2012 R2 platform, I have updated the latest
hotfixes for time zone information(Microsoft KB 2981580.
http://support.microsoft.com/kb/2981580).
But still R is not populating correct date time values(Standard and
Daylight saving).
Could you please advise me, how to
Hi,
I would like to know, once the operating system timezone information is
updated, what step should be carried out at R for reflecting the operating
system(LINUX X64 SOLARIS X64 and WINDOWS X64) update.
Thanks
Vasanth
[[alternative HTML version deleted]]
Log out and log back in again. For many people this may happen in the normal
course of daily use.
---
Jeff NewmillerThe . . Go Live...
DCN:jdnew...@dcn.davis.ca.usBasics:
I am log out and login many times. but no luck.
Thanks
Vasanth
On Fri, Nov 7, 2014 at 1:04 AM, Jeff Newmiller jdnew...@dcn.davis.ca.us
wrote:
Log out and log back in again. For many people this may happen in the
normal course of daily use.
R on Windows uses the Olsen timezone database, a copy of which is stored with R
in the Program Files directory (e.g. R/R-3.1.1/share/zoneinfo). You could
update the file yourself if you can find a corrected version, or download an
updated version of R.
Thanks Jeff Newmiller, it is working now.
I would like to know, the same kind of configuration can be done in Linux
and Solaris platform.
Instead of R is mapping to operating system(/usr/share/. /usr/share/lib/)
zoneinfo directory.
Thanks
Vasanth
On Fri, Nov 7, 2014 at 1:18 AM, Jeff Newmiller
Working now... after what action?
AFAIK on *NIX systems R uses the OS installation of the Olsen database, so on a
fresh login R should pick up any OS update you have installed.
---
Jeff NewmillerThe
In LINUX, I don't want R to use the operating system zoneinfo(Olsen
database) instead of that I like to point different path which has the
latest zoneinfo(latest Olsen database).
On Fri, Nov 7, 2014 at 2:34 AM, Jeff Newmiller jdnew...@dcn.davis.ca.us
wrote:
Working now... after what action?
?timezones
You probably need to recompile R.
---
Jeff NewmillerThe . . Go Live...
DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...
Loops are not slow, but your code did a lot of unneeded operations in each
loop.
E.g, you computed
D$id==i D$t==t
for each row of D. That involves 2*nrow(D) equality tests for each of the
nrow(D)
rows, i.e., it is quadratic in N*T.
Then you did a data.frame replacement operation
D[k,]$y
From: silvano silv...@uel.br
I fitted a genetic model using kinship2 and I compared it with MTDFReml program
output.
The residual variance of both are very close but the genetic variance are very
differents.
The output are:
MTDFReml:
genetic variance = 1.24015
residual variance = 5.93424
R
Hello,
a day or two ago I submitted VGAM 0.9-5 to CRAN, which has
myriads of changes to family functions---their names, their
arguments, and their order thereof. Especially regarding
family functions for discrete and continuous distributions.
In a nutshell, I found lots of inconsistencies while
Hello, all!
So, as stated in the title, the Lme4 package used to output p-values for the
fixed effects. What happened?!
Literally 2 weeks ago, I ran code, got output with no errors, and had
p-values listed for my fixed effects.
Now, running THE SAME CODE with THE SAME DATASET (nothing at all
Hi,
I noticed your email on the R support page. I am currently using SPSS 21 for
mac, and need to install a macro plug-in which requires me o install R 2.14 (I
already have R 3.0). When i download the R 2.14 package from the CRAN website,
and try to install it on my Macintosh HD, it gives me
Hi,
I noticed your email on the R support page. I am currently using SPSS 21 for
mac, and need to install a macro plug-in which requires me o install R 2.14 (I
already have R 3.0). When i download the R 2.14 package from the CRAN website,
and try to install it on my Macintosh HD, it gives me
On 06.11.2014 20:21, Deepansh Dalela wrote:
Hi,
I noticed your email on the R support page. I am currently using SPSS 21 for
mac, and need to install a macro plug-in which requires me o install R 2.14 (I
already have R 3.0). When i download the R 2.14 package from the CRAN website,
and try
Daniel Mello dmello2 at ucmerced.edu writes:
Hello, all!
So, as stated in the title, the Lme4 package used
to output p-values for the
fixed effects. What happened?!
Literally 2 weeks ago, I ran code, got output with no errors, and had
p-values listed for my fixed effects.
Now,
On 06.11.2014 20:20, Dalela, Deepansh wrote:
Hi,
I noticed your email on the R support page. I am currently using SPSS 21 for
mac, and need to install a macro plug-in which requires me o install R 2.14 (I
already have R 3.0). When i download the R 2.14 package from the CRAN website,
and
SNIP
On Thu, Nov 6, 2014 at 2:05 PM, Matteo Richiardi matteo.richia...@gmail.com
wrote:
SNIP
Final question: in your code you have mean(M[t-1L,]): what is the 'L'
for? I removed it at apparently the code produces the same output...
SNIP
The constant 1L is stored as an integer; the
Hola,
Los ejes los controlas con scales. Es una lista, y uno de sus
componentes puede ser y. Es a su vez otra lista en la que puedes
configurar el eje con detalle. En help(xyplot) está descrito en un
apartado que empieza así: scales: Generally a list determining how
the x- and y-axes (tick marks
Me uno a los agradecidos!
El 6 de noviembre de 2014, 10:35, Lorena Tudela Marco
lorenatudelama...@gmail.com escribió:
Muchas gracias por compartir! ;)
Un abrazo, y buen día.
Lore
El 4 de noviembre de 2014, 22:33, Igor Sosa Mayor
joseleopoldo1...@gmail.com escribió:
Estimada Lorena Tudela
Voy a escribir algo por como yo lo aprendí (quizás mal redactado para un
estadśitico), lo observado es lo que se mide a campo y lo predicho es lo
que surge a partir del modelo estadístico. Luego aprendí como realizarlo
con R, pero no siempre es posible en forma fácil,
Funciono perfecto, muchas gracias Oscar.
Debo decir, que lei esa parte de la ayuda de xyplot, pero jamas se me
habria ocurrido que servia para controlar lo que necesitaba.
Slds y gracias,
Eric.
On Thu 06 Nov 2014 05:49:06 CLST, Oscar Perpiñan wrote:
Hola,
Los ejes los controlas con
Lorena
Ni idea:
Parte 1, nlme puede calcular con predict..., tendría que leer las
especificaciones de nlme, porque dentro de las opciones puede ser que ...
Parte 2, tantas librerías juntas, ..., puede ser que la última que carga
sobrepone una función a las anteriores.
Yo no puedo opinar sobre
Hola Lorena.
Trabajé hace tiempo con binomiales negativas y esta página me sirvió
bastante: http://www.ats.ucla.edu/stat/r/dae/nbreg.htm
Ahí tienes una manera de sacar los predicted values para una negbinom. No
se si es exactamente lo que buscas, pero te puede servir como punto de
partida para
Hola Gerardo, antes hice algo parecido para una presentacion a publico
general, y nadie se dio cuenta, pero ahora tengo que enviar ese grafico
a una revista cientifica.
Gracias.
Slds, eric.
On Thu 06 Nov 2014 14:20:39 CLST, Gerardo Gold-Bouchot wrote:
Puedes sumarle 1 a todos tus
Es un truco muy com�n en varias disciplinas, particularmente en Ecolog�a,
as� que no deber�a haber problema. De hecho en la Biblia de Ecolog�a
Cuantitativa, de Legendre y Legendre, se recomienda sumar un n�mero que sea
del mismo orden de magnitud del valor de la media op mediana de la
variable,
39 matches
Mail list logo