search for: hadassa

...nt in each row) and I would like to convert this matrix to a binary one according to the minor allele frequency. For non-geneticists: I want to have a binary matrix for which in each row the 0 stands for the less frequent character and 1 for the more frequent character. Thanks for any suggestions. Hadassa -- Hadassa Brunschwig PhD Student Department of Statistics The Hebrew University of Jerusalem http://www.stat.huji.ac.il

...have someone already posted this question but I still have not figured how I could do this: Say I have var.1 var.2 var.3 which are large matrices. I would like to do the following loop: for(i in 1:3){ # Pick var.i and do something with it } How can I tell R to do this? Thanks for any comments. Hadassa [[alternative HTML version deleted]]

...size = theta Theta, however, is not an integer. So how can size (which is the number of successes) equal theta? I know there is a function rnegbin which does what I want. I would still like to raise the above question (probably a more statistical one than R). Thanks a lot for any comments. Hadassa [[alternative HTML version deleted]]

...with the variable PKG_LIBS=-L/directory_of_tools -l/tools? 2. Does this file need to be located in the same directory as main.c? 3. When I then run the command R CMD SHLIB it does not make use of Makevars. So what do I need to add? 4. Does tools.c need to be a shared object as well? Thanks a lot! Hadassa -- Hadassa Brunschwig PhD Student Department of Statistics The Hebrew University of Jerusalem http://www.stat.huji.ac.il

...<- seq(0,0.9,by=0.1) redshade <- rgb(red=1,green=a, blue=a) However, I don't really know how to plot things using hexadecimal colors. I used a function which tries to find the closest color to the rgb shades but it didn't work very well. Any comments would be appreciated. Thanks, Hadassa [[alternative HTML version deleted]]

...,1,1,1,1,1,1,1,1,1,1,1,1),C0=5,st90=4,C0.pop=5,st90.pop=4,tau.cpop=0.2,tau.stpop=1) mcmcA <- bugs(dataA,inits,parameters,modelA,debug=T,n.chains=7,bugs.directory="c:/Program Files/WinBUGS14",working.directory="C:/Documents and Settings/Daikon/Roche") Thanks so much... Hadassa -- Hadassa Brunschwig Birmannsgasse 10A CH-4055 Basel Switzerland Phone: +41 78 797 6065 Email: h.brunschwig at utoronto.ca

...meters of the extreme value distribution are equal for each stress level. So in order to perform the likelihood ratio test I need the likelood of the model under HA: location parameters are not equal (i.e. each stress level has its slope and intercept). How can I do this? Thanks for any help! Hadassa

...on might be a basic one. I did browse the help and archives but I still haven't understood how to do the following: I run R locally and would like to read in data from a server for which there is a username and password. That is, how do I open a connection to a server with a password? Thanks, Hadassa -- Hadassa Brunschwig PhD Student Department of Statistics The Hebrew University of Jerusalem http://www.stat.huji.ac.il [[alternative HTML version deleted]]

...shell during the run. I would like to run this program from R and as far as I have seen system() does not handle manual options (?). Is there thus a way to use these C programs anyway through R (i.e. without touching the C programs and fixing the problem in R)? Thanks a lot for any comments. -- Hadassa Brunschwig PhD Student Department of Statistics The Hebrew University of Jerusalem http://www.stat.huji.ac.il

...way to control the strip size in xyplot for example using the strip.default function. Or do I need to redraw the strips by myself? I would like to make the strip slimmer and I have reduced the font size using cex. However, I don't know how to change the size of the strip. Thanks for any tips. Hadassa [[alternative HTML version deleted]]

Hello everyone I have data grouped by subjects and i used the trellis plot to show the curve for each subject. I have another vector "a" which gives me the highest points the curves for each subjects can achieve ( the curve is exponential ). Now i want to draw a horizontal line, whose values are saved in "a", on each plot for each subject. Is there any possibility to do that?

Thanks to everyone who replied to my first problem. Here is the second one: Is the function "lme" (Mixed Model) also in a foreign library(shouldnt, as this is a basic statistical function)? I am trying to use lme in R and it replies that it couldnt find the function lme...weird... Thanks for replying!

Hello everyone I've got another question on Trellis plots. I plotted my fitted model with plot(augPred(fitted model)). According to the description of this method, one can use any optional argument which is passed down to the xyplot function for Trellis plots. Additionally to the trivial plot i would like to insert a horizontal line. I tried to do it with the argument panel= function (...)

...dimensions, the more rows in the array. So for example the array could look like the following: , ,1 1 2 3 4 , , 2 5 6 7 8 9 10 11 12 , , 3 13 14 15 16 17 18 19 20 21 22 23 24 Has anybody a comment on how to build this in the most efficient manner (I know i could do looping)? Thanks! Hadassa

...d Z respectively. The output of lme gives another estimate called "range" which I assume belongs to the parameter estimate needed for the Gaussian correlation. Are my assumptions correct? And where can I get the variance for the W(t) from? Thanks for any answers...and happy new year... Hadassa

...168 0.016 0.003 -0.008 -0.018 -0.008 0.250 I(age^2) 0.146 -0.010 -0.010 0.009 -0.007 0.005 -0.156 -0.917 Standardized Within-Group Residuals: Min Q1 Med Q3 Max -4.22819684 -0.45649041 -0.02996715 0.43379089 2.70363057 Thanks for any help. Hadassa

Hi, I know that with readlines you can start executing a code upon user input. What I would like to do is to have the ability to abort a code based on a new user input. E.g. when the user hits enter the code stops running. That is, I would periodically check for a user input and stop the code when the correct input was entered. Below is an example of a function which would run an endless loop

...to where the 100 of plot2 is such that the it looks evenly distributed. But that means that an interval between tick mark 1 and 2 of plot1 (which is 5) is as long (graphically) as the interval between tick mark 1 and 2 of plot2 (which is 10). Any ideas on how to do this? Thanks a lot!! Cheers, Hadassa [[alternative HTML version deleted]]

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems

Im trying to import data from an excel sheet or a sas file to R...im not succeeding. Apparently the function read.xport for reading a SAS file doesnt exist. What do i have to type in EXACTLY to read from an excel sheet(i guess i would be using read.table?)? Thanks in advance for an answer Dassy