search for: nrcan

...que Canadian Forest Service / Service canadien des forêts Natural Resources Canada / Ressources naturelles Canada Government of Canada / Gouvernement du Canada 506 West Burnside Road / 506 chemin Burnside ouest Victoria, BC V8Z 1M5 Phone (250) 363-3342 Facs: (250) 363-0775 Email: Alex.Song@nrcan.gc.ca <mailto:Alex.Song@nrcan.gc.ca> [[alternative HTML version deleted]]

...eck if the association shape for a specific copula is more appropriate for your application. If this is not possible fit models assumig different copulas and check which one provides a better fit to your data. I hope it helps. Best, Dimitris Quoting "Chaouch, Aziz" <achaouch at NRCan.gc.ca>: > Hi Dimitris, > > I'm not sure to understand your suggestion. How would you build that > contour plot for a particular copula and on what is computed the > Kendall's tau? > > Thanks, > > Aziz > > -----Original Message----- > From: Dimitr...

Hi, I'm dealing with wind data and I'd like to model their distribution in order to simulate data to fill-in missing values. Wind direction are typically following a vonmises distribution and wind speeds follow a weibull distribution. I'd like to build a joint distribution of directions and speeds as a VonMises-Weibull bivariate distribution. First is this a stupid question? I'm

...y PC clock to the samba server. I still get a lot of messages from the Studio informing me the source code files and workspaces have changed outside the editor or Studio and asking me if I wish to re-read the files. Any advice? Jim -- Love is like the ozone layer, |Jim Craven craven@cg.NRCan.gc.ca| You don't miss it until it's gone. | Geological Survey of Canada | ---Frank Drebbin, Police Squad| 615 Booth St., Rm 211, Ottawa | --------------------------------------| K1A 0E9 (613) 996-9935---------|

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

...right? Does this correspond to SAS's R2? Many thanks, Yan Ressources Naturelles Canada Service Canadien des Forêts - Centre de Foresterie des Laurentides 1055, rue du PEPS CP 10380, Succ. Ste-Foy Québec, QC, G1V 4C7 Tel. : +001 418 649-6859 Fax : +001 418 648-5849 email : Yan.Boulanger@nrcan.gc.ca [[alternative HTML version deleted]]

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

This may seem to be an odd request, but here it goes: I would like to be able to recursively load pxelinux. Simple reason being that many useful mini distro's now use it (PXES, Recovery is Possible, Fedora Anaconda, my own homebrew setup program, unattended.sf.net, etc), and I would like to be able have drop in upgradabily for those distro's. The idea being I can create a master menu and

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

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)

...e protocol version 1 seems to behave as expected, while version 2 appears to look at the password file despite "PasswordAuthentication no" in the config file ? Or am I just missing something important ? I'm not a regular subscriber, I'd appreciate comments to rosenberger at pgc.nrcan.gc.ca or andreas at arescon.com -- ___________________________________________________________________ a. rosenberger arescon ltd. andreas at arescon.com 9706 First St. www.arescon.com Sidney, B.C. Canada V8L 3C7 ____________...

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