similar to: assessing the fit of a LME model

Dear List, I'm getting weird and unexpected behaviour using time and %%, or %/%. It's likely I'm not appreciating the nuances of floating point arithmetic. Or it could be a bug. I'm running > R.version _ platform x86_64-redhat-linux-gnu kernel 2.6.33.8-149.fc13.x86_64 arch x86_64 os linux-gnu system x86_64, linux-gnu

There is a bug in Matrix package, please check it, thanks! Matlab result: x = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 >> lu(x) ans = 21.0000 22.0000 23.0000 24.0000 25.0000 0.0476 0.9524 1.9048 2.8571 3.8095 0.7619

I cannot understand the following error printed out when I try to get the extrema of my time series. I would appreciate some suggestion as I really cannot interpret the error. I might not be using a proper set of parameters in calling such functions. I am learning by doing ... > aa.peak <- wavCWTPeaks (aa.tree) Error in `row.names<-.data.frame`(`*tmp*`, value = c("1",

There seems to be 2 functions call ecdf... http://lib.stat.cmu.edu/S/Harrell/help/Hmisc/html/ecdf.html http://127.0.0.1:11885/library/stats/html/ecdf.html How do I get the one ecdf {Hmisc} to run instead of the ecdf {stats} A pointer in the right direction would be greatly appreciated. Tried to instal Hmisc but got this message, so I assume I have it > utils:::menuInstallPkgs() Warning:

I'm doing the time-series data analysis using FFT. After fft(), I used spectrum() to calculate the power spectrum, which is used to define the main frequency. Now I want to get the phase information according to this frequency. For example, input 12 time-point data,

Dear list; I just created a utility function that replicates what I have done in the past with Excel or OO.org by putting a formula of the form =sum($A1:A$1) in an upper-corner of a section and then doing a "fill" procedure by dragging the lower-rt corner down and to the right. When divided by the grand sum of the entries this function then calculates a 2D-discrete-ECDF. I

Am I trying to perform multiple linear regressions on each 'VARIABLE2'. I figured out that there are different ways, using the following code: (data is given at the end of this message) reg <- lapply(split(TRY, VARIABLE2), function(X){lm(X2 ~ X3, data=X)}) lapply(reg, summary) Which produces the following: $`1` Call: lm(formula = X2 ~ X3, data = X) Residuals: Min

All - I think I'm being dense, but for the life of me, I can't figure out why I get error message with the code below. I have data that looks like param level perc.surv asin.tran DO 3 0.6864407 0.9764544 DO 3 0.1250000 0.3613671 DO 3 0.8738739 1.2077299 DO 4 0.4615385 0.7468986 DO 4 0.5691057 0.8547258 DO 4 0.8504673 1.1737517 DO 6 0.5494505 0.8349297 DO

I have a dataframe that I'm trying to reshape, and need advice. My data: > klam.merge[200:225,] stream lulc x sumlength pct.lgth 200 1223030419685 92 0.25000000 9.89 2.52780586 201 1223030419686 23 0.00274154 4.73 0.05796068 202 1223030419686 41 0.75009917 4.73 15.85833341 203 1223030419686 42 2.65000000 4.73 56.02536998 204

Eric ¿Que dataset utilizo? Por curiosidad probé su código, pero me salen errores. Copio y pego todo como me sale en la consola. > # https://raw.githubusercontent.com/datasets/covid-19/master/data/time-series-19-covid-combined.csv > datos <- read.csv("C:/Users/HP/Downloads/time-series-19-covid-combined.csv") > > library(data.table) > library(ggplot2) >

Muchas gracias a todos por su ayuda. Hace un rato antes de abrir el correo, he encontrado dos soluciones a mi problema de como calcular el porcentaje de variación. La primera es usando el paquete dplyr: https://stackoverflow.com/questions/48196552/calculate-percentage-change-in-r-using-dplyr/48196871 La segunda usando la función PercChange del paquete DataCombine

Dear R-helpers, I have a data.frame (bcpe.lat.m) containing 13 countries, ages 0-50yrs per month, and the corresponding mu&sigma (see below). * I would like to limit the age range to include all 12 months for the 1st 5 years and only whole years for all ages thereafter for each of the countries present in the data frame. * I would like to create separate data.frames according

Hello, i have following problem and I hope you can help me a little bit My dataframe looks like: df a m d typ value 1950 1 1 5 -4.1 1950 1 2 9 2.7 1950 1 3 3 -1.3 1950 1 4 5 -1.9 1950 1 5 2 0.2 1950 1 6 8 0.5 1951 1 1 4 1.3 .... It consists by daily observations from 1950- 2009. Now, I get with.... for (i

In the new version 0.90.0, rbind won't take a vector and a matrix from me, but works OK if I coerce the vector to a matrix. The following was run on a new session (i.e., no prior work), and has been duplicated on 2 machines: an SGI running Irix 6.5, and an Intel box running Red Hat Linux 6.0. Either case works fine in version 0.65.1. > t3 <- c(.5, .5) > t4 <-

Hi I have been looking around for a couple of days without being able to find a solution. I am trying to plot the following table using the package ggplot2 Month Mean SEM cry3m 3 -0.05227273 0.1607422 cry6m 6 0.65826087 0.6638734 cry1y 12 0.38333333 0.3062907 cry2y 24 1.04811111 0.4921863 cry3y

Gracias Ruben, muy interesante. Algo así estuve buscando, por ciudades y a nivel internacional, pero no lo encontré. Me bastaría con el número de fallecimientos por día, de cuantas más ciudades del mundo mejor. Si alguno supiera de algo así, le agradecería que me lo comunicase. Un saludo, Manuel El mar., 24 mar. 2020 a las 21:04, Rubén Fernández Casal (< rubenfcasal en gmail.com>) escribió:

In the new version 0.90.0, rbind won't take a vector and a matrix from me, but works OK if I coerce the vector to a matrix. The following was run on a new session (i.e., no prior work), and has been duplicated on 2 machines: an SGI running Irix 6.5, and an Intel box running Red Hat Linux 6.0. Either case works fine in version 0.65.1. > t3 <- c(.5, .5) > t4 <-

Full_Name: Madeleine Yeh Version: 1.9.1 OS: AIX 5.2 Submission from: (NULL) (151.121.225.1) After compiling R-1.9.1 on AIX 5.2 using the IBM cc compiler, I ran the checks. One of them failed. Here is the output from running the check solo. root@svweb:/fsapps/test/build/R/1.9.1/R-1.9.1/tests/Examples: ># ../../bin/R --vanilla < stats-Ex.R R : Copyright 2004, The R

Dear R users, Please can you help me with a relatively straightforward problem that I am struggling with? I am simply trying to plot a baseline survivor and hazard function for a simple data set of lung cancer survival where `futime' is follow up time in months and status is 1=dead and 0=alive. Using the survival package: lung.wbs <- survreg( Surv(futime, status)~ 1, data=lung,

Dear all, I have used following code but everytime I encounter a problem of not having coefficients for all the variables in the predictor set. # code rm(list=ls()) library(caret) # generating response and design matrix X<-matrix(rnorm(50*100),nrow=50) y<-rnorm(50*1) # Applying caret package con<-trainControl(method="cv",number=10) data<-NULL data<- train(X,y,