search for: dyang

Dear all; Using read.ssd, I tried to import a SAS dataset from a network drive; an error occurred: > stemattr <- read.ssd("g:/olmn2/dyang/siteprod/datasasv8", "stemattr") SAS failed. SAS program at C:\DOCUME~1\dyang\LOCALS~1\Temp\file7785.sas a log and other error products should be in the vicinity Warning messages: 1: sas not found 2: ls not found 3: SAS return code was -1 in: read.ssd("g:/olmn2/dyang/sitepr...

R-helpers; I tried to get a superscripted 3 in the following xyplot example but failed: >data(whiteside) >xyplot(Gas ~ Temp | Insul, whiteside, panel = function(x, y, ...) { panel.xyplot(x, y, ...) panel.lmline(x, y, ...) }, xlab = "Average external temperature (deg. C)", ylab = paste(paste("Gas consumption (1000", expression(ft^3),")"), aspect

Dear All; While using EMclust() in the mclust package, I frequently encountered a program error. A message window popped up with the message " Rgui.exe has generated errors and will be closed by Windows. You will need to restart the program. An error log is be created." > version _ platform i386-pc-mingw32 arch i386 os mingw32

Dear R list members; I have a 10 x 423 data frame which consisting of response, time, subject, site, plot and covariates (continueous and categorical) measured at the plot level. When the data frame was converted into a groupedData object, a warning appeared > A <- groupedData(ht ~ time | Subject, data = tt, outer = ~ site * plot, + labels=list(y = "Height", x =

Dear all; I experienced a memory allocation error while running gnls: > wBlplfe.wc <- update(wBlplfe.w, corr=corAR1()) Error: cannot allocate vector of size 37950 Kb In addition: Warning message: Reached total allocation of 256Mb: see help(memory.size) Increased memory limit to 512, 640 M, the allocation error persisted: > memory.limit(512) NULL > wBlplfe.wc <-

R-Helpers; This seems simple to set graphic parameter, but I tried plot a graph with 10 intervals in x axis with par(xaxp = c(0, 150, 10) but to no avail: par(xaxp =c (0, 150, 10)) plot(age, y18, type="n", ylab="Height (m)", xlab=" age (Yrs)", font = 2) lines(age, y6, lty = 4, lwd = 2) lnes(age, y10, lty = 1, lwd =2) lines(age, y14, lty = 2, lwd =2)

Dear R-help; Following the installation of R 1.5.0 and library packages to my Win 2K machine, I tested the new R system: > library(nlme) Loading required package: nls > library(lattice) Loading required package: grid >data(Soybean) >plot(Soybean) A program error appears: "Rgui.exe has generated errors and will be closed by Windows. You will need to restart the program. An

Dear all, I have a data frame with one factor and four numeric variables and wish to obtain the var-cor matrix separately by factor. I tried by() and sapply() but getting nowhere. I understand this can be done by subsetting the dataframe, but there should have some sleek ways of doing it. Here is a simulated dataframe; s <- rep(c("A","B","C"), c(25,22,18)) d

Hello, I run R on Windows 2000. When I load lattice library and plot some lme class objects, the graphic window displays a dark-green color. Graphics look ok on the screen, but lines and points become invisible whem I send the plot to a printer because of the background color. Are there ways to reset the background color? Both par(bg=) and windows() do not work. Any ideas? Richard

Hi R-help; In using nlme, I often encounter the following error message: Error in nlme.formula(model = response ~ ......, data = ... : Step halving factor reduced below minimum in PNLS step Resetting the pnlsTol default value does not solve the problem. Any suggestions or workarouds for solving the problem? TIA, Richard

Dear R-help; I created a nlme model (Xnlme) with 3 fixed effects B1, B2, and B3. The random effects and available covariates in the groupedData were extracted with raned(): >Xnlme.ranef <- ranef(Xnlme, aug = T) When the fixed effect was plotted against covarites, an error occurred: >plot(Xnlme.ranef, form = B1 ~ ps + dr + el + sl + as) an error occurred: "Error in

Many thanks to Peter Dalgaard BSA and Uwe Uwe Ligges for their prompt responses. It turns out xaxp can not be set by par(xaxp = c(x1, x2, n)) as indicatedd in Help for 'par' document. Both pointed me to the right direction by setting axes = F or xaxt = "n" in plot() and using axis(1, seq=c(0,150,15)) to produce an axis with desired intervals. Richard My original question:

Dear All; Searching on the achieve, many questions on optim() have been asked, but I haven't seen the following. The question began with my original inquiry on "Optimization failed in fitting mixture 3-parameter Weibul l distribution using fitdistr()" which I posted on Jul. 28, Prof. Ripley kindly advised me to look into options of optim() for the answer. Following his advice and

Dear All; Attempting to reproduce Figure 4.15 of MEMSS (p. 171) in R using plot(Wafer, outer = ~ Wafer) yields an error: > plot(Wafer, outer = ~ Wafer) Error in order(na.last, decreasing, ...) : Argument lengths differ The plot() produces the figure without problems for all versions of Splus (4.5 to 6.2) on Windows. Noticing that the plot() differs slightly between

Dear R-helpers; While using gnls() to fit a function > Gbht0t.gnls <- gnls(h2 ~ Rht(b0, b1, b2, h1,t1, t2), data=gbht10, + params=list(b0 + b1 + b2 ~ Sisp -1), start=c(strssb0,strssb1,strssb2)) I encountered an error: "Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid variable type " Rht is a defined function to be

Dear all; This an FAQ. I tried to access lme source script so I can step into it to debug the problems resulting from a lme() call. I used getAnywhere("lme") or nlme:::lme, both produced only the function definition and "UseMethod("lme"). Any idea how to list the source code? TIA, Richard Yang [[alternative HTML version deleted]]

R-helpers; I used Proc Model in SAS to fit the following model to data: proc model data = dbsmv; a = a1*F**2; b = b1*F + b2*T + b3*F*T; tph2 = tph1 *((1 - exp(-a*age2)) / (1 - exp(-a*age)))**-b; fit tph2; and yielded the following estimated parameters after iterations: a1 = -0.15943, a2 = -1.8177, b1 = -0.01911, b2

Dear All; I have used Winedt editor in conjunction with R on the Win2K platform for more than a year without any problem. Recently I purchased a P4 machine with Windows XP. Following the installation of R1.62 and WinEdt 5.3 to separate directories, I copied R-WinEdt to the Plugins subdirectory under the Winedt directory, clicked on the "install" file and edited the Rprofile with two

Dear all; I wish to create a graphic object combing an xyplot() and an mtext(). My code looks like following, gmv <- { trellis.device("windows", bg="white", width = 7, height = 7) xyplot(Mvol ~ Age | Nl * Th , data = Hft1, prepanel = function(x, y) prepanel.loess(x, y, span = 1), xlab =list(label = "Age (Years)", font = 2),

Dear all; Installing the GLMMGibbs package to my Solaris Unix box, I got an compiling error: ars.c:497:10: missing terminating " character ars.c: In function `dump_arse': ars.c:498: error: parse error before "mylesj" ..... The compiling error was reported to the list on Jul 3, 2003. According to Prof. Brain Ripley this is a known problem with the package and gcc 3.3,