similar to: extra panel arguments to plot.nmGroupedData {nlme}

Dear R useRs, I have a problem with nls.lm function of minpackl.lm package. I need to fit the Van Genuchten Model to a set of data of Theta and hydraulic conductivity with nls.lm function of minpack.lm package. For the first fit, the parameter estimates keep changing even after 1000 iterations (Th) and I have a following error message for fit of hydraulic conductivity (k); Reason for

Hi, I am using the lme package to fit mixed effects models to a set of data. I am having a difficult time understanding the *meaning* of the numDF (degrees of freedom in the numerator), denDF (DF in the denomenator), as well as the Intercept term in the output. For example: I have a groupedData object called 'Soil', and am fitting an lme model as follows: ## fit a simple model #

Hi, I am curious about how to interpret the table produced by anova(ols(...)), from the Design package. I have a multiple linear regression model, with some interaction, defined by: ols(formula = log(ksat * 60 * 60) ~ log(sar) * pol(activity, 3) + log(conc) * pol(sand, 3), data = sm.clean, x = TRUE, y = TRUE) n Model L.R. d.f. R2 Sigma 1834 1203

I'm trying to write a function that has a dataframe as an argument. Within the function, I call groupedData to create a grouped-data object from the dataframe argument. However, the groupedData function doesn't seem to see the data. I'm guessing it's getting the argument name rather than its value, but I'm not positive. My system: R 1.3.1 on WinNT 4.0 (Pentium II), nlme

Dear all, I'm using the nlme package. I would like to introduce the plot of my groupedData object in my latex report using Sweave but it doesn't work. I have no error message in R (Sweave(my.file_Rnw)) and the plot is generated in R, but when I want to see the pdf (I'm using TexnicCenter) of my report, there is an error message and no pdf file is generated. It seems that the plot

Hi, I am trying to install the nlme package, but during installation I get the message <snip> pdSymm text html latex example plot.ACF text html latex example missing link(s): xyplot plot.Variogram text html latex example missing link(s): xyplot plot.augPred

Dear List, This may be related to this email thread initiated by Ben Bolker last month: https://stat.ethz.ch/pipermail/r-devel/2011-July/061630.html In answering this Question on StackOverflow http://stackoverflow.com/q/7195628/429846 I noticed that `methods()` was not listing some S3 methods for `plot()` provided by the mgcv package. At the time I wanted to check the development version of R as

Dear Group: I have been trying to write an R function within which the function groupedData() would be used. The following is a sample program: ####################### sid<-rep(1:6,times=2) time<-c(1:12) trt<-rep(letters[1:3],times=4) tem<-data.frame(sid,time,trt) test1.fun<- function(dat=dat) { groupedData(time ~ trt|sid, data=dat) } test1.fun(dat=tem)

Hello There, I tried your advice and i don't think it works either, See the logs, " /tmp/nut# ./drivers/usbhid-ups -DDDDDD -d1 -s test -x port=auto -x vendorid=... -x productid=... -x subdriver=... Network UPS Tools - Generic HID driver 0.52 (2.8.0-2557-g81fca30b2) USB communication driver (libusb 1.0) 0.46 ?? 0.000000??? [D3] main_arg: var='port' val='auto' ??

hello, I am very new to R. My current data set is a mix of values and categories. It is a geoscience data set, with values per rock sample. Case in point, each sample belongs to a lithology class, and each sample has several physical property measurements (density, porosity...). I want to be able to plot these physical properties for all samples in each lithology class. this is how i'm doing

<facepalm ;)> the ellispis "something=..." were supposed to be replaced by your device ID numbers, etc. On Sat, Sep 23, 2023, 18:26 FatGear <fatgear1 at free.fr> wrote: > Hello There, > > I tried your advice and i don't think it works either, > > See the logs, > > " > > /tmp/nut# ./drivers/usbhid-ups -DDDDDD -d1 -s test -x port=auto -x

For an axis label I want to include the Greek letter mu within the string. I've not found the proper way of including that expression within the string. What I want is "Conductivity (uS/cm)" with the 'u' replaced by mu. When I try "Conductivity (" expression(paste(mu)) "S/cm)" I get an error. If I don't separate Conductivity and S/cm with

hello mailing list! i still consider myself an R beginner, so please bear with me if my questions seems strange. i'm in the field of biology, and have done consecutive hydraulic conductivity measurements in three parallels ("Sample"), resulting in three sets of conductivity values ("PLC" for percent loss of conductivity, relative to 100%) at multiple pressures

I am having difficulty understanding how I would constrain a data set by filtering out 'records' based on certain criteria. Using SQL I could query using 'select * from my.data where LithClass in ('sand', 'clay')' or some such. Using subset, there seem to be ghosts left behind (that is, all of the LithClass *.Labels* remain after subset) > dput(tcc)

Dear all, I have tried to calculate a GLMM fit with lmer (lme4) and glmmPQL (MASS), I also used glm for comparison. I am getting very different results from different functions, and I suspect that the problem is with our dataset rather than the functions, but I would appreciate help in deciding whether my suspicions are right. If indeed we are attempting the wrong type of analysis, some

The data is of two forms, ie one numeric and another ordinal ie 1.5 and <.5. The X Value is a factor. I'm trying to create a lattice conditioned plot with the following characteristics: 1) The plot is conditioned using the form (Conductivity~Day|Valve) 2) The plot should use a barplot for the ordinal (<.5) and dots for the numeric (1.5) 3) A line should be created specifying a

Hello, I have a time series dataset in long format (15 minute data over approximately 2 years). The data set is composed of time series of multiple sensors; much of the time series are duplicate dates and times, but of different lengths. I can plot the data converting the datetime to as.numeric, but they do not plot if I keep the data as POSIXlt. I know xyplot would plot the data if it was a

Nilesh M. Vaghela Research Scholar c/o Prof. N. V. Sastry, Department of Chemistry, Sardar Patel University Vallabhvidyanagar-388 120 Gujarat, India --- On Thu, 29/7/10, Kurt Hornik <Kurt.Hornik@wu.ac.at> wrote: From: Kurt Hornik <Kurt.Hornik@wu.ac.at> Subject: Re: R statical-script for analysis (India) To: "Nilesh Vaghela" <nilesh_spu25@yahoo.in> Cc:

On Tue, 8 Sep 2015, Rob Groner wrote: > I executed lsusb to verify the USB device is there, and it is. I tried > the shutdown command again with debug enabled, but it didn't seem to > reveal much more: > ----------------------------------------------------------------------------------------------- > rtd at linux-5048:~> sudo /usr/local/ups/sbin/upsdrvctl -DDDDD shutdown

Dear all, this is probably more a statistics question than an R question but probably there is somebody who can help me nevertheless. I'm running a regression with four predictors (a, b, c, d) and all their interaction effects using lm. Based on theory I assume that a influences y positively. In my output (see below) I see, however, a negative regression coefficient for a. But several of the