search for: maggy

Dear R-Users, I want to ask a question for a colleague of mine. He wants to put a character vector and a numeric vector into one matrix and still have the old character and numeric type for the respective columns. Unfortunately, I am just starting using R and I could not help him. Is there an easy and straightforward way to do this in R? Maybe a little example facilitates understanding our

Hi, I had a hard time in learning nlm in R and appreciate any help. I encounted the following error message from time to time when I tried different starting parameter values (three parameter values in this case) in nlm(f=SS.fun,p=c(0.1/40,0.1,2),hessian = FALSE,N.measure=object,h=20) Error in f(x, ...) : only 0's may mix with negative subscripts Basically I know the three parameter

Dear List, With regard to the question I previously raised, here is the result I obtained right now, brglm() does help, but there are two situations: 1) Classifiers with extremely high AIC (over 200), no perfect separation, coefficients converge. in this case, using brglm() does help! It stabilize the AIC, and the classification power is better. Code and output: (need to install package:

Dear R-users, I use glm() to do logistic regression and use stepAIC() to do stepwise model selection. The common AIC value comes out is about 100, a good fit is as low as around 70. But for some model, the AIC went to extreme values like 1000. When I check the P-values, All the independent variables (about 30 of them) included in the equation are very significant, which is impossible, because we

my dataset looked like this in the beginning: >Daten V1 V2 V3 1 Dosis Gewicht Geschlecht 2 0 6.62 m 3 0 6.65 m 4 0 5.78 m 5 0 5.63 m I need box plots for V2 with all combination of V1 and V3, so I deleted the first row, and tried this: boxplot(Daten$V2[Daten$V3=="m"]) but it does not work and I

Hi, I run the following tuning function for svm. It's very strange that every time i run this function, the best.parameters give different values. [A] >svm.tune <- tune(svm, train.x, train.y, validation.x=train.x, validation.y=train.y, ranges = list(gamma = 2^(-1:2), cost = 2^(-3:2))) # where train.x and train.y are matrix

Hello, I am now running a multiple linear regression program, but I do not know the difference between the command step and stepAIC. Thanks. Maggie [[alternative HTML version deleted]]

Hi, I am writing a program to plot the correlation coefficients of 45 elements using the commend "plot". Error message "Error in plot.new() : figure margins too large" appeared. I think it may because there are too many elements. But I really need to find their correlation coefficients. What can I do? Thank you. Maggie [[alternative HTML version deleted]]

Hello, I can't link a c code with Mathlib according to introduction of R manual "Writing R Extensions", page 60. It is written : "It is possible to build Mathlib, the R set of mathematical functions documented in 'Rmath.h?, as a standalone library 'libRmath? under Unix and Windows. (This includes the functions documented in Section 5.7 [Numerical analysis

Hello, everybody, I created a R package which includes C code. But I load this package, and carry out the R function in it. It shows C function is not in load table as follows. Would you tell me what is the problem? Where do I make mistake? Maggie [Previously saved workspace restored] > library(var) Attaching package 'var': The following object(s) are masked _by_

I have the following code that aligns the two graphs. Problem is that in .pdf it gives me it x-axis (0-100) is broken down into 0-20, 20-40..and so on. I wonder if there is for it to display the x-axis (and y-axis) in more detail than that. I'd appreciate your input -- pdf(file="VECTOR & ICA ALIGNMENT.pdf", height=5, width=5)

Hello, I am trying to run a virtual machine with qemu-kvm and libvirt as a user. In order to start the virtual machine I setup a tap device manually, add it to the bridge and activate it via ifconfig. But I cannot tell libvirt to not try to add the definded tap device to the bridge by itself. It always fails telling me it cannot add tap0 to br0. Thats the xml snippet for defining the network:

I can''t get past the setup step in Capistrano. I know the paths are correct and I''m using zsh locally (bash on remote), so the shell shouldn''t be a problem. It never asks for the SSH password, even. Any help you could give would be greatly appreciated. Thanks. Kenneth -- => the blog from beyond <= => www.eyeheartzombies.com <=

The plastics industry says polycarbonate bottles are safe. http://www.bisphenol-a.org/about/faq.html#g I'm sure Maggie and here friends would say ALL plastic bottles are very dangerous. This lady seems to be at a reasonable middle ground. http://trusted.md/blog/vreni_gurd/2007/03/29/plastic_water_bottles Polycarbonate plastics the kind of bottle you bought contains BPA. "In 2006 Europe

my data looks like this: PM10 Ref UZ JZ WT RH FT WR 1 10.973195 4.338874 nein Winter Dienstag ja nein West 2 6.381684 2.250446 nein Sommer Sonntag nein ja Süd 3 62.586512 66.304869 ja Sommer Sonntag nein nein Ost 4 5.590101 8.526152 ja Sommer Donnerstag nein nein Nord 5 30.925054 16.073091 nein Winter Sonntag nein nein Ost 6

Resend the mail. ???: Zhan, Cindy (???) ????: 2015?8?28? 16:19 ???: nut-upsdev at lists.alioth.debian.org<mailto:nut-upsdev at lists.alioth.debian.org> ??: Xie, Maggie(???) ??: Request for ew device<0x2B2D/0xFFFF> listed into NUT Dear NUT Developers: We really appreciate for your great job and accomplishment in NUT. It?s interesting and fabulous. It would be great if

I want to make a stacked bar plot with one bar for two variables from my data "chir", the two variables have about 100 values like no, yes and na. I want to show how many no, yes and na they both have together with the stacked bar. I tried to melt these to variables first like this: melt1=melt(data_chir, measure.vars=c("N1_re", "N2_re"), var="zpd") but it

Hi, I have been able to generate association rules for Market Basket Analysis using the following codes: **************************************************************************** ******************************************* library("arules") rules <- read.csv("write1.csv",na.strings=c(".", "NA", "", "?"),header=TRUE)

I've been working on adding code coverage support to my ELLCC project (http://ellcc.org) and I thought it would be nice to fun the index page of the HTML output of llvm-cov into something a little prettier and more useful. The initial result is here: http://ellcc.org/coverage/ If there is any interest, I'll submit the llvm-cov changes. -Rich

Hello, I have some data covering contaminant concentrations in fish over a time period of ~35 years. Each year, multiple samples of fish were taken (with varying sample sizes each year). Ultimately, I want an estimation of the variance between years, and the variance within years + random effects. I used a linear mixed model to estimate these variances, but after reading a number of different