similar to: nls weights warning message

All the usual MASS data sets are nicely sitting in D:\Program Files\rw1022\library\MASS\data on my laptop and yet: > library(MASS) > data(wtloss) Warning message: Data set `wtloss' not found in: data(wtloss) what could be wrong? I'm running 1.2.2 on Win98. Murray Jorgensen -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list --

I'm running R 1.2.2 under windows 98 on a Pentium 133 laptop. I can't seem to retrieve the package MASS data sets: > library(MASS) > data(wtloss) Warning message: Data set `wtloss' not found in: data(wtloss) > data(abbey) Warning message: Data set `abbey' not found in: data(abbey) And yet all the .rda files for the MASS datasets are in D:\Program

Hello dear R-helpers, I'm working with R-2.15.2 on Windows 7 OS. I'm stucked with a merge of two data frames by characters. In each data frame I got two different list of names, that is my main-key to be merged. To figure out what I'm saying, I build up a modified "?merge" example, with errors by purpose: # Data for authors: authors <- data.frame( surname =

Hi, In the third edition of the MASS book, Chapter 8, Page 252 ~ 253. I tried applying the same technique on using a negative exponential model (as suggested in Exercise 1), and trying to get a confidence interval. What I have changed is follows: expn <- function(b0, b1, th, x) { temp <- exp(-x/th) model.func <- b0 + b1 * temp # Differentiate temp derive.negexp <- (x

Hi, I am trying to figure out how to lable the second y-axis after the following codes: plot(x, y, xlab="time", ylab="pay1" ) par(new=TRUE) plot(x,y2, ann=FALSE, xaxt="n", yaxt="n", pch=7 ) axis(4) Then, I want to label the second y axis "pay2". I tried "title(ylab="pay2")", but it put this lable on the first

Ben Bolstad experienced something odd while running R CMD check over a package of ours using r-devel (yesterday's fresh). He tracked down the problem to a vignette calling 'library(MASS)'. Please kindly ignore if this is a known (and temporary) issue. Laurent ----- Forwarded message from Ben Bolstad <bolstad@stat.berkeley.edu> ----- Return-Path:

Dear R users, I am trying to print heatmaps in a loop (with a pause). Idea is to visualize changing correlations over time and for testing I wrote this simple (reproducible) code below. My problem is that levelplot() does not produce any output when I run the code (though heatmap does). Ideally I would like to use levelplot() as it produces a neat index on the side indicating the color and the

Dear all, I am using the R 2.4.0 environment on Windows XP SP2 machine and trying to use the lattice package version 0.14-9 which you have kindly written to share with the R community. I have a question concerning the stripplot which I'd be very grateful if you can kindly advise me on: I used the dataset called trydata to plot a graph using stripplot (dataset attached), and the scripts

Hi, I'm running Aaron's tunequeue scripts for my radio site. (tunequeue.sf.net) I have ices2 to read Ogg Vorbis files and feed the icecast2 server with those files. Could someone please tell me how can I disable ID3 updating? I want all clients to only show the name of the radio, but don't get any of the ID tag information stored in the Ogg files. Where can I set that? Sorry for

I knew it wasn't ID3, I just didn't know the name for it. This way you knew what I meant. :) Thanks for the reply. Then I'll modify my Ogg files. I have jlGui Java player embedded into my web page, and it keeps reconnecting every time a song has ended. Unfortunately by reconnecting it looses the first 1-2 seconds of the songs. According to what you've said there's no way

It's very possible that the Polycom IP600 will work with this. As it is just an implementation of a SIP standard for subscribing to the state of other extensions. As for the feature improvements you might see them from me, but not very likely. It is easier for me to train my customers to hit *8 (I will probably just program a pickup button for them) than it is for me to figure out what I

The snom 200 and 220 have five programmable buttons. Each button has a led that can be used to indecate if an extension is idle, in use, or ringing. A button pannel for the 220 is also comming out soon that will have 20'ish programmable buttons on board. This is a killer app for any company that has receptionists handle calls, and pretty usefull for everyone else. As a matter of fact,

Hello! I have a data frame with some NAs. test<-data.frame(a=c(1,2,NA),b=c(10,NA,20)) I need to sum up values in 2 variables. However: test$a+test$b procudes NAs in rows that have NAs. How could I sum up columns while ignoring NAs (the way the function sum(..., na.rm=T) works? Thank you! -- Dimitri Liakhovitski marketfusionanalytics.com

hi, we want to replace our nt4 pdc/bdc with samba and we use windows nt4 terminal servers with citrix in our environment. as far as i know, we will get a problem here. (extensiver reading of the mailing list led me to this conclusion) regarding: http://ma.ph-freiburg.de/tng/tng-users/2001-05/msg00224.html >Terminal Server has its own User Manager (which has an extra button -

Hi all, I'm having some difficulty passing arguments into lm() from within a function, and I was hoping someone wiser in the ways of R could tell me what I'm doing wrong. I have the following: lmwrap <- function(...) { wts <- somefunction() print(wts) # This works, wts has the values I expect fit <- lm(weights=wts,...) return(fit) } If I call my function lmwrap, I get

Dear Community, I have a data set with two columns, bird number and mass. Individual birds were captured 1-13 times and weighed each time. I would like to remove those individuals that were captured only once, so that I can assess mass variability per bird. I¹ve tried many approaches with no success. Can anyone recommend a way to remove individuals that were captured only once? Thanks, Ray

Hello, I've been working with "nnet" and now I'd like to use the weigths, from the fitted model, to iterpret some of variables impornatce. I used the following command: mts <- nnet(y=Y,x=X,size =4, rang = 0.1, decay = 5e-4, maxit = 5000,linout=TRUE) X is (m x n) Y is (m x 1) And then I get the coeficients by: Wts<-coef(mts) b->h1 i1->h1

Hi, I am doing double for loops to calculate SDs with some weights and wondering if I can get rid of the outer for loop as well. I made a simple examples which is essentially what I am doing. Thanks for your help! -Young #------------------------------------------------------ # wtd.var is Hmisc package # you can replace the 3 lines inside for loop as # sdx[i,] =

Hi. I am running SAMBA 2.0.7 on HP-UX 11.00 and HP-UX 10.20. I have tried adding the MultipleUsersOnConnection registry key onto my Microsoft Windows NT Server 4 Terminal Server Edition servers, running with Service Pack 6. After I make this change in the registry and reboot my WTS servers all my samba connections from the servers are still being made with one process instead of being split

I could use some help understanding how nls parses the formula argument to a model.frame and estimates the model. I am trying to utilize the functionality of the nls formula argument to modify garchFit() to handle other variables in the mean equation besides just an arma(u,v) specification. My nonlinear model is y<-nls(t~a*sin(w*2*pi/365*id+p)+b*id+int,data=t1,