This likely isn't a research-level question, but it is at least a question of interest to this researcher. I'm happy with an answer that sends me somewhere (preferably online) to read about a well-known (to somebody) solution, if there is such a thing. First the question, then the motivation.



> 
> Let $X\_1,X\_2,\dots,X\_N$ be independent $N(\mu,\sigma^2)$ random variables 
>  ($\mu$ and $\sigma^2$ are unknown), and let $Y\_1,Y\_2,\dots,Y\_N$ be defined by 
>  $\{X\_1,\dots,X\_n\}=\{Y\_1,\dots,Y\_N\}$ and $Y\_1 < Y\_2 < \cdots < Y\_N$. 
>  For $k < N$, how can one estimate the parameters $\mu,\sigma^2$ from $Y\_1,Y\_2,\dots,Y\_k$?
> 
> 
> 


The application that I need this for is a strategy to buy light bulbs. I just purchased a stunning chandelier with 27 tiny light bulbs. The bulbs themselves are strange, and I'll have to special-order them. After the first several burn out, I'd like to be able to estimate how many will burn out in, say, the next 6 months. Although the bulbs are independent, the lifetimes of the first 3 to burn out are not.


This problem must come up all the time: medical researchers don't know how long their patients will last, they only know how long some of them didn't. 



How to estimate a time distribution
As requested, I'm making my comment an answer.


While this does not answer the problem as you've formalized it, I would argue that the correct method to consider this problem is actually Weibull analysis--generally failure times are not distributed normally. Essentially, plot the failure times on Weibull paper; the slope of the line you get will give you the information you want.


As Kevin O'Bryant has pointed out, [this site](http://www.weibull.com/AccelTestWeb/calculating_the_parameters_of_the_weibull_distribution.htm) does Weibull computations.


For a beautiful theorem explaining the relevance of the Weibull distribution, see the [Fisher-Tippett-Gnedenko theorem](http://en.wikipedia.org/wiki/Fisher-Tippet-Gnedenko_theorem). The model to keep in mind is the following: objects generally have more than one possible point of failure, and failure occurs when the *first* point of failure breaks. So in an ideal case one wants to compute the minimum of a collection of IID's. The Fisher-Tippett-Gnedenko theorem gives that the extreme value distributions, of which the Weibull distribution is an instance, are universal for this situation, just as the central limit theorem gives that the normal distribution is universal for the *mean* of a family of IID's.


